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RELATIVE MOTION OF ORBITING PARTICLES 
UNDER THE INFLUENCE OF 
PERTURBING FORCES 

ByJ.Bo Eades, Jr. * 

SUMMARY 

This report describes the mathematical developments carried out for this 
investigation. In addition to describing and discussing the solutions which were 
acquired, there are compendia of data presented herein which summarize the 
equations and describe them as representative trace geometries. 

In this analysis the relative motion problems have been referred to two 
particular frames of reference; one which is inertially aligned, and one which 
is (local) horizon oriented. In addition to obtaining the classical initial values 
solutions, there are results which describe cases having applied specific forces 
serving as forcing functions. These forces are designated as having components 
parallel to the triads of the reference frames ; hence the analytical results ob- 
tained are for motions referred to both frames and influenced by both forcing 
function systems. Also, in order to provide a complete state representation the 
speed components, as well as the displacements, have been described. These 
coordinates are traced on representative planes analogous to the displacement 
geometries. By this procedure a complete description of a relative motion can 
be developed; and, as a consequence range-rate as well as range information is 
obtained. 


♦Senior Analyst and Manager, Seabrook Office, Analytical Mechanics Associates, Inc. 
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I. INTRODUCTION 


Two bodies in motion, and in close proximity to one another during a 
space flight follow trajectories which are relatively independent of one another. 
Contrary to this, when physical connectors are in place, or "attractive forces" 
of sizeable magnitude are present, then these trajectories may be coupled so 
that the ensueing motions are dependent on one another. 

The investigation reported here has considered several aspects of the 
relative motion problem. In particular, the influence of initial state (position 
and velocity), and the consequence of selected, applied force systems have been 
examined. The applied forces assumed in this study were designed to have specific 
orientations; the consequence and purpose for this will become evident as the 
results are brought forth and discussed. As a prelude to what follows herein, 
it is noted that these individual forces, which are presumed to be acting on one 
of the bodies, have components of fixed magnitude parallel to the reference 
triads. By this mechanism it is demonstrated that analytical solutions may 
be obtained for a "properly described" equation of the motion. Subsequently 
the relative motion descriptions which evolve will exhibit an influence from 
these particular forces. Also, in relation to this, the reader may recognize 
that a variety of simulated physical situations could be described in terms of 
these same actions. In this regard it is apparent that the results from this 
investigation represent a significant addendum to the present collection of in- 
formation for the relative motion of orbiting particles. 

As an aid to clarifying the physics of the present situation the sketches 
shown below indicate the two particular frames of reference to be utilized. In 
both drawings the motion commences with the reference triads parallel to one 
another. However, as time progresses the base particle (P) moves to a new 
location, as does the relatively moving particle (Q). Throughout this time the 
two masses may remain separated, as suggested by the relative displacement 
vector (r). However, the main difference, as noted on these sketches, is that 
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in one instance the displacement is referred to a moving frame of reference 
(x, y, z) - see Sketch (A) - while in the second (Sketch (B)) it is referred to an 
inertially aligned triad. 


Q 



Incidentally, it is in reference to these triads that the applied forces are 
assumed to have their fixed magnitude components. And, it is these forces 
which are studied to learn how the relative state of motion is affected. In 
substance, then, it is apparent that the position and velocity, at any time, re- 
ferred to either triad, is a consequence of three input quantities. One, the 
initial state of relative motion (this case is hereafter called the "initial values 
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problem"); and secondly, the two separate applied force systems (this will, here- 
after, be referred to as the "zero initial values problem"). 

. In the discussions and descriptions which follow, here, there will emerge 
some definite patterns regarding the relative motion state. Most noticeable are 
the reactions produced by the applied force systems and by the initial state. In 
particular the consequences of these inputs on motions referred to the two triads 
of reference, are primarily seen in the in-plane displacements and velocities. 

As a general statement these in-plane relative motion components (displacements 
and speeds), when referred to a rotating frame of reference, trace out figures 
which are counter to the orbital motion. On the other hand, corresponding situa- 
tions when referi’ed to the inert ially oriented frame of reference, trace out 
figures which are in the direction of the orbital motion. Secondly, the initial 
values problems are represented by different geometries in these two frames 
of reference. In particular, for general inputs the geometries which result 
show a divergence - as the time in motion progresses. However, the traces 
representing these various cases are not alike when described in terms of 
the various coordinates. The figures in a rotating frame are generally cycloidal 
and/or parabolic in character. On the other hand, the motions expressed in 
inertial (relative motion) coordinates are nominally spiral-like. This would 
seem to imply that the relative state of motion is not bounded, except as it is 
related to the central attracting mass. Fortunately this is not the case; there 
are means available which will provide a "closed" relative motion path. The- 
conditions by which this may be brought about are described in the report under 
f, non-secular" cases, 

. As a word description, illustrative of this last condition, consider the 
following situation. Suppose a sub-satellite is to be 'launched" from an orbiting 
vehicle; and, also, suppose it is desired that this sub-satellite remain in the 
immediate vicinity of its parent. It can be shown that when the parent body is 
In a circular orbit, one means of attaining the desired relative motion is to 
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impulsively launch the sub-satellite orthogonal to the parents velocity vector. 

Of course, the geometric shape of the subsequent traces will appear differently 
on the planes for the two reference frames. For example, the in-plane trace 
i (referred to the rotating frame of reference) would appear as an ellipse*. On 
the inert ially oriented plane of reference the corresponding figure would trace 
as a limaeon*. (Necessarily displacement traces alone do not fully describe ; 
the relative state of motion — one needs the velocity diagram, or hodograph, 
as well. These are readily obtained and easily plotted; once the full state of 
motion is described the investigator has a complete graphical representation of 
the problem r s actions). 

One may wonder about the reasons behind the two reference frames. This 
is easily reconciled when it is recognized that some spacecraft may be Earth- 
oriented while others could be star— oriented; thus the two coordinate systems. 

In addition there is, occasionally, some question regarding the need for both 
the displacement and the hodograph traces of a given motion. Generally this is 
most easily explained by noting -that the state of motion (absolute or relative) 
is comprised from both coordinates. Also, there may be a distinct need for 
range and range-rate information, especially in conjunction with maneuvering 
operations. Needless to say this information can be readily obtained from these 
relative motion coordinates; and, in addition, the graphical display of both can 
enhance the understanding of how contemplated flight operations might progress. 

A few words regarding the out-of -plane motions would be worthwhile at 
this point. As a general observation it can be verified that this aspect of the 
relative motion, be it numerically defined, analytically described or graphically 
represented, is more complicated and less easily defined than the in-plane case(s). 
From an analytical point of view the main difficulty comes about from the over- 
all dependence of the representative equations on the full initial state and/or all 
the force vector components appearing in the input. As is demonstrated in the 

♦These situations are described in the body of the report. 
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report, it is indeed difficult to estimate (say) the geometry of these motion 
traces as an a priori piece of information. With some study, and a few example 
situations, it is reasonable to expect that the reader could acquire some ,f feel” 
for these problems; and he could develop a sense of prognostication with regard 
to the motion traces, but only in a general manner. 

For those readers who are not too well acquainted with the ideas and 
applications of a relative motion study it would be worth the time and space 
needed to mention a few examples. Probably the first example to come to mind 
would be that of intercept and rendezvous between two orbiting vehicles. Whether 
the situation should consider a one active -one passive pair, or two active 
satellites, is of little consequence. When the two vehicles are cooperating, the 
rendezvous problem is straight forward and easily implemented (mathematically). 
In the event that one body is '‘uncooperative” (an "evade” situation) the rendez- 
vous is more difficult to achieve, but remains a problem in relative motion - one 
which fits the overall concepts of this investigation. With the option, here, of 
adding "thrust” to the relative moving vehicle it is evident that (say) the rendez- 
vous can be accomplished under this action rather than by impulsive means. 

The use of sub -sate llites, experiment packages, and other orbiting 
particles come directly under the relative motion problem concept when one 
seeks to track these bodies and/or predict where they might be at some subse- 
quent time. 

It should be remembered that in general the studies which could come 
under this classification need not be considered simply as "external” motions. 
That is, the equations and procedures developed for this investigation can be 
utilized to study the motion of particles inside a spacecraft as well as outside 
of the vehicle. By properly interpreting the physics of many problem situations, 
and adapting these results to the conditions at hand, the ideas set down here 
have a wealth of applications possibilities. One somewhat trivial but interesting 
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example which comes to mind is that of a particle free to move about inside 
a spacecraft. Recognizing that the most significant force which would be acting 
(there) is that due to gravity gradient, then it is evident that the particle cannot 
remain motionless — in a relative sense — unless it is located at the dynamical 
center for the spacecraft -particle system. Actually the particle is most likely 
to be, and remain, in motion as a consequence of the gradient force. Necessarily, 
the closer the particle is to this center the smaller the acceleration acting on the 
body. It can be shown that the trajectories for such a particle, referred to a 
moving frame of reference, would be regular hyperbolae. The particle would 
tend to move toward* the dynamical center and in the gradient direction simultan- 
eously. 

Before leaving this section and moving to the mathematical developments, 
it would be well to mention some pertinent facts regarding the following sections. 
First, in the writing of many formulae in this report the matrix notations which 
appear may be somewhat unusual to the reader. In this regard it is suggested 
that those who may be interested in the details of developments, leading to and 
represented in these various expressions, should consult the Appendices as a 
first reading assignment. These will introduce the reader to the notation and 
make the reading of the mathematics much easier and more understandable. 

Also, the developments herein make ample use of matrix notation; and, the 
solutions, per se, are carried out by the method of Laplace transforms. Those 
who are unfamiliar with either or both of these notational operations may find 
some difficulty in following the mathematics. 

Finally, in conjunction with discussions of the various types of motions, 
and those inputs which influence them, the reader will find compendia of re- 
sults sections. These data presented collections of results which are typical 
of the various motion classifications, as described. The purpose in presenting 


*This suggests an idealized situation wherein the spacecraft flies a circular orbit. 
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II; KINEMATICS AND PROBLEM DEFINITION 


Kinematics (An Inertial Frame of Reference). In this section certain kinematical 
relations are developed to describe a relative motion situation as it would be 
viewed in an "inertially oriented’ 1 frame of reference. 




P(0) 

P 



M 


Fig. II. 1 - Geometry of a re- 
lative motion. The vectors P, 

Q locate the two particles (of 
interest) with respect to M 
(origin). The relative position 
vector (r) locates "Q” with res- 
pect to "P". 


It is assumed that the (~)°-frame 
represents the inertial reference system. 

On the other hand the (~) f -frame undergoes 
angular displacements (<p) for the regular, 
orbital motion (particle P). Let an adjacent 
particle (Q) have its motion defined relative 
to P. Then, in this discussion, all motions 
are described (kinematically) in the inertial 
frame. 

' Accordingly, particle ”Q” moves 
relative to ”P T ', as both move about M (the 
attracting center). Here r^ locates r, Q fr 
with respect to ”P”, thus its position 
vectors is : 

Q = F + V 

which can be described in either frame of 
reference. 

Assuming that r^ is measured in the 
local inertial frame, and since ”P” (origin) 
moves along its own orbit, then the location 
of ”Q” in the inertial frame is given by, 



is a general designation for the relative position vector. 

| Preceding pagesblank 



Qj^H +T(<p ) P. 


01 . 2 ) 


(Herein R is used to denote the inertially defined relative position vector; 
r will be reserved as the description for r in a rotating frame of reference, 
one attached to "P". Also, (~) will be attached to P and Q to denote their 
inertially defined values. The transformation matrix, T(<p ), transfers the 
vector P (<p), back to the inertial frame, (~)°; see Appendix A for a discussion 
and development of the transforms). 


Kinematic Definition of Velocity. From Eq. (II. 2) above, it follows that the 



Fig. n.2. This figure is a sketch 
showing the inertial frame, 
attached to P, used to locate *Q , T 
with respect to ,r P M . 

i 

[ 


velocity of Q, defined in inertial space, is 
(generally) given as: 

Q j = R + (T_ ) P + (T_) P . (H. 3 a )+ 


But, since Tf P ri is assumed to move on a 
circular path, and to have only planar ro- 
tation (|p| is not changing), then it follows 
that, 

• • 

Q r “R + (T-) P, (II. 3b) 



-sin <p 

-cos (p 

0 

• • 

where (T_)=*p 

+cos ip 

-sin <p 

0 


0 

0 

0 

m 



b(t.) 2 

o 

II 

O 

1 

1 

o 


(see Eq„ (A. 6a), 
Appendix A) 


Here (T_)^is a "reduced matrix", as de 
fined, in the Appendix. 


*(If |P | would be allowed to change, then this full expression would be used. 
Also; here, and elsewhere, the notation (T_) means T(<p“); a(T + ) will be used 
(subsequently) to denote the transform, T(<p + )). 
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Malang use of this description for T, then a more convenient form of 
Eq. (II. 3b) can be had; substituting Eq. (A.6d) leads directly to, 



Q 


I 


= R + <p (BTJ P, 


wherein B is the constant matrix defined in Eq, (A. 6c), Appendix A. 

This last expression shows that the (kinematic) velocity of ’'Q" is com- 
posed of the relative (but inertially described) velocity (R ), plus the velocity of 
T, P" transferred back to the inertial reference system. 


Kinematic Acceleration. The "acceleration" for "Q" is readily evaluated from 

Eq. (II. 3c), by differentiation. That is, as a general statement, 

0 0 0 

^ = R + ^(BTJ P + (p J&TJ P +<p <p [BT- ] F ; 

here the deleted quantities are a consequence of the invariance of terms 
(<p, |p|, and B). Consequently, the applicable kinematic description is, 

«• *• 

Qj =R +<p [BT_] P, 

** 

or, recalling that BB = - I , 

2 

• * ^ 

Qj-R-tp (I 2 TJ P. 

One should note that the (kinematic) acceleration for "Q" is composed- of the 

relative value (R) - for the displacement of "Q" with respect to T, P M - plus a 

• 2 

centripetal term, proportional to <p . 

A reorientation of Eq. (11.4a) shows the relative acceleration (R), ex- 
pressed in the inertial frame of reference, written as: 

- 1 * 9 - 

r = q j + i a 2 Tjp. 


*The symbol denotes equations of importance. 

**See Appendix A. 


(II.3c)* ** 


(II. 4a) 


(II. 4b) 
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Here, the quantity (Q^) is symbolic of the "specific accelerations" arising from 
the forces acting on the particle (Q); while the second term is a centrifugal 
force component associated with the position of "P"l 

Concluding Remarks . The expressions developed above describe a means for 
representing a particle r s relative motion, as this is expressed in an inertial 
frame of reference. However, this case is restrictive in that it pertains to a 
regularly circulating base particle (P). Of course, the analytical operations 
introduced above are fundamentally geometric (or kinematic) in kind. In order 
to apply these relations to a particular case it would be necessary to introduce 
the "dynamical" (Newton mechanics) aspects of the situation — that is, to des- 
cribe the "forces" which are acting on the test particle (Q). 


Kinematics, for a Rotating Frame of Reference. In the next developments, to 
follow, ’hnotions" will be referred to a moving frame of reference — one which 
has a leg parallel to the radius vector locating the base particle (P). In this 
rotating-frame vectors P and r will be defined relative to this moving triad; 
hence, in order to refer these vectors back to the inertial axis system, it 

x° is required that the transform T(<p“) be 
applied to both. (See the sketch below 
for a geometric representation of this case). 


Here r^ ~ r (x 1 , y\ z r ) and P - 
|p| e f ; hence, the position vector 

Q p r + P, in general) 
r 

is written as, 



^)P(O) 


M 


Fig. II. 3. Sketch depicting the 
moving frame of reference. 


or 


Qj = T (<p ) r + T(<p ) P, 


Qj ~ T (<p )(r + P). 


♦Here r describes the relative position vector; the lower case letter is used 
t o distinguish from the previous development. 


(II.5)* 
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As before, the transformation matrix, (T_), is applied to rotate the vectors 
r, P back to the (e^)° triad. 

It must be remembered that these vectors are expressed in the rotating 
frame; consequently, they relate to the moving triad. 


^1. 6a) 


(IL 6b) 

recall that (T_) means T(<p~). This is the transformation used to rotate vectors 
back to the (~)°-frame of reference. 

An examination of Eq. (II. 6b) indicates that the "local 11 velocity (r), and 
the associated rotational component(s), are added vector ally (but transformed 
also! ). It is worth noting that the rotational component is proportional to the 
orbit angular rate (<p), and to the vector position of T Q’ f , which is (r + P). 

Necessarily, to transform this part of the velocity back to the inertial frame 
the operator (BT_), which is a geometric entity, is utilized. 

In the next paragraphs the. kinematic acceleration is defined. 

Kinematic Acceleration . Carrying the differentiation one step further leads 
directly to the acceleration definition. From this operation the quantity is 
obtained as: 


Kinematic Velocity Definition. Differentiating Eq. (II. 5) to obtain the velocity 
of ”Q n - as it would be referred to the inertial frame - one obtains the general 
expression: 

0 

• ■ . • y 

Qj = (T_)[f + P] + (T_) [r.yp], 

wherein (P) vanishes due to the assumed Invariance of |p| . In addition, 
since (T_) = <p B(T_) (see Eq. (A. 6a), Appendix A) the velocity is most con- 
veniently expressed as: 

• • 

Q = (T_) [ r ] + (p- (BT_) [ r + P] ; 
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0 0 

Q = <Tj[r] + (TJ [r]+#T [BTJ [r + P] +<p + BTJ [r + p] 

0 

+ <p (BT_) [r +/1 , 

wherein the vanishing components are due to those invariants previously 
mentioned. 

Making substitutions, and collecting terms, one finds that: 

** •• * p n 

Q. = (TJ r + 2 <p (BTJ r - <p (I TJ [r + P] . 

I a 

In this expression the three components (shown) are the familiar quantities: 

(a) (r), the local acceleration, 

« «— 

(b) (2<p r ), the Coriolis acceleration, 

and 

(c) (<p [r + P] ), the centrifugal term. 

Equation (II. 7) defines the kinematic acceleration; this is a geometric 
description; however, in solving the relative motion problem it will be necessary 
to rearrange the expression in order to acquire an equation for r. That is, 

(from Eq. (II. 7)), one writes: 

f« 44 . • _ 

r = (T_ ) Q - 2 <p (T_ BT_) r + <p (T_ I 2 T_) [ r + P] , 
wherein 

<T S 1 = T(0“) T = T(<p + ); (see Eq. (A. 5), Appendix A). 

Evaluating the matrix products in Eq. (II. 9a), it is seen that (T + BT_) 

= B 2 and (T^I^T^) = I^; consequently, 

•« •« * _ 

^ f = f (T+) Qj } - 2<p (B 2 ) r + <p fl 2 ) [r + P ] . 


(n.7) 


(II. 8) 


(II. 9a) 


(II. 9b) 
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In this last expression all quantities, except Q , are referred to the local rotating 

«• i 

triad. The parameter (Q ) is a symbolic representation for the specific accelera- 
tions due to outside (or external) forces. E for a central field attraction case 
Q would be replaced by the gravitational force, and/or the other forces which 
would influence the motion of "Q Tr . These particulars will be discussed in more 
detail subsequently. 


Motion of the Reference Particle (P). In order to examine the motion of M P M , due 
to a central field of force, one should consider (first) the dynamical equation: 

7. (GM)m p 


which assumes an inertial frame of reference. 

Now, considering the kinematics of this motion, return to Eq. (II. 9b); 
but,there,let r => 0, thence "P", and 

_ # •• •• ■ _ 

q“p =(T + )Q 1 -^> (B 2 )? +<P(I 2 >[? + p 3}= s ' 

{(T + )P t + p 2 d 2 )p}; 

therefore : 

<r + ) p = - <p 2 a 2 ) p . 

(Note: here the quantity (T+) P represents the transforming of P^ to the 
displaced position, and to the reference frame located by <p). 

Now, using Eq. (II. 10a) to describe P — which applies solely to a 
simple field of attraction — then Eq. (II. 10b) leads to 


GM *2 



A , L . . .n -j- _l j ■ " i --j._ p r- * 

*It must be evident that (p represents a (constant) angular speed for the 
reference particle ,, P n as it moves on its circular orbit. 


(II. 10a) 


(II. 10b) 


(II. 10c)* 
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* 

This last result is the well known expression for the angular rate ( 0 ) on a 
circular orbit (|p| ) due to the gravitational attraction (GM). This result 
states that the circular velocity, squared, is defined by the quotient GM/ |p |. 

The Relatively Moving Particle (Q), Having determined a motion for "P", 
the next situation to consider is that for Tt Q”. Of course, this will quickly lead 
to a description of the relative motion as it develops under the influence of 
whatever force system(s) assumed. In this regard suppose (for the moment) 
that only a gravitational force is present; then, 

1 IQjI 3 1 

where m is understood to represent the mass of the particle located at "Q"; 
and where Q is its position vector (relative to M). Recalling that 
Qj S (T_) [r+P], see Eq. (II. 5), then it follows that 

i- cr-)Cr- + P]. 

jr + P | 

(This models the specific force, or acceleration, acquired from a mass 
attraction, M). 

In attempting to evaluate this non-linear expression it is evident that 
an analytic solution cannot be had; however, a linearization (based on the 
smallness of |r | with respect to | P | ) will allow one to find an approximate 
analytical result. (It can be said that this formulation, acquired above, de- 
fines the "gravity gradient" effect for Q ! ) 

An esqpansion for the denominator term, above, is outlined below: 

Since 

|r + P| - J (r + P) • (r + P) , 

then 

*The transformation matrix, (T_J, is not included in the denominator since it 
does not influence the magnitude of the vector (sum). 


(11.11a) 


(IL 11b)* 
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and 


r + P| 3 - 


|p|3r i + (j£[) 2 +2 

|P[ 


r ° P 13/2 

p 2 J 


i? + pf 3 - ip r 3 fx -| (-J|L ) a -s ]. 


1*1 


(II. lie) 


Q, 


Now, making use of Eq. (ILllc), in Eq. (II. lib), it follows that; 

- r • P 3 - 2 


GM 

IpI 3 


cr_, { [(? + p> - 3P ] - .3? - f fr + p) + ' ' • ' ) ; 


p| 


p' 2 


H. O.T, 


and, to first order in (r), the approximate expression is: 

► ¥-^‘ 3 (rj [ (?+2) - 32 ^]- 


(It should be noted that as r “* 0 (Q => P), then 


Q => - 


GM 



(T J P . 


(II. lid) 


This is equivalent to Eq. (II. 10a). As before, this specifies a motion for "P" 
as it would be described by Newtonian mechanics. A consequence of this last 
result is that the two remaining terms in Eq. (II. lid) are "additions" to that 
influence attributed to the attraction of M). 


Of the two remaining "force" components, in the equation above, the 
one quantity; 



GM 

|P| S 


(TJ r = 


• 2 - 
- <P <T_) r , 


(II. lie) 


has a same form as before. Therefore it behaves as a "restoring force"; i.e. , 
a force tending to return ’Q"' to the position of "P" during the motion. 
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However, the remaining quantity, 

A Q = + 3 M. <TJ P 1= 3(^ 2 (I.) P [? • P ] , (II. Ilf)* 

1 IPI 3 L IPI 2J 

has a different character (as compared to the others); this is a force which acts 
opposite to the "attraction” of gravity. Consequently, since it is radially directed 
(ocP), it behaves like a "reduction” in gravity and has the tendency to drive the 
particle away from the origin (at P). 

These two partial results define the basic relative motion of TI Q”. It is 
apparent that one part of the motion is restorative while the other part tends to 
the opposite effect. This combination leads to an interesting physical interper- 
tation of the motion, for Q, in the vicinity of P. Subsequently, a discussion of 
this will be made. 

Equations of Motion . Having obtained a kinematic expression for the relative 
motion acceleration, Eq. (II. 9b), and having developed a (linearized) des- 
cription of those effects attributed to gravitational attraction, Eq. (II. lid) - 
then, a combination of effects* leads to a descriptive estimation of the relative 
motion. That is: following Eq. (II. 9b), 

« * t» * _ 

' ?=(T + )Q r -2p <B 2 ) r +< 0 I 2 (? + P) f 

wherein, now, 

(T+) Q T ~ - — 3 <T+T-) [(? + P) - 3P . 

1 |P| 3 L |P| 2 J 

Taking into account that, 

Q O 

CT+TJ = I , andGM/ |p| =<p (see Eq. (11.10c)), 

4 4 

then r can be written as : 

** O ^ ^ • 

r = <p [(I, - I )<r+P) + 3P (r • P)] -apB 2 r . (II. 12a)** 

♦The description herein Includes a force system aligned with the triad of reference. 

A — 

** P is a unit vector corresponding to P . 
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In this expression P is a unit vector pointing toward tr P’ T from the center 
of attraction (M). 

To be more precise with the format of Eq. (II, 12a) it should be written 
as shown below. (Recall that in this equation r = r (x, y, z), where the scalar 
components are described in the moving frame of reference). As a consequence, 
Eq. (II. 12a) is expressed as: 



wherein P=|p| e t =|P| J , and (r • P) = x. Finally, then, 

X X 



(II. 12b) 


Eqs. (11.12) describe a (linearized) form of the equations for the (Q) 
particle's relative motion. This motion is with respect to an orbiting particle 
moving on a circular path about M. Of course, the equations are written in 
terms of Rotating" Cartesian coordinates; also they are approximations for 
the case wherein the test particle (Q) is acted on solely by the gravitational 
attraction, of M. 
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In a "real life" situation the elementary mass (Q) would be acted on by 
forces other than gravity, also. In order to represent these additive effects, 
let any such specific accelerations be represented (symbolically) by a vector, 

F. Of course, to be compatible with Eq. (11.12), it is necessary that F be 
described in the moving frame of reference. Consequently, a more complete 
description of Q’s motion is given by the equation: 

{>» ? =<p [(I 2 -y(r + P) + 3I 3 P (r- P)] -2<PB 2 ? (II. 12c) 

3 

In Eq. (II. 12c) the F. would account for all of the perturbative forces which 

j 

could influence the relative motion of the test particle. For the present these 
forces are not described (other than symbolically); any particular descriptions 
will be left for the reader to describe, as needed. 

Discussion on the Linearized Equations . Most likely the clearest descriptions 
of the motion for Q, as given in Eqs. (11.12) can be had by examining Eq. 

(II. 12b), wherein the additive perturbations () F.) may or may not be included. 

However, when the perturbing quantities are included, they should be described 
in a manner such that they do not unduely complicate the solution procedure. 

Haying gone to the trouble of ’linearizing" the gravity gradient term, one should 
not model the perturbations in more than a linear format. 

It might be expected that a general solution to these equations would 
evolve as an initial value problem; for such cases the results should be ex- 
pressed in terms of elementary functions. 

From Eq. (II. 12b), it is seen that the in-plane (x, y) scalar differential 
equations are coupled through the Coriolis terms while the out-of-plane ex- 
pression is independent of the other two. This suggests that the one scalar 
solution can be (conceptually) acquired independent of the others; however, the 
couple equations must be handled together and simultaneously. (One approach 
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to a typical solution will be suggested subsequently). Of course, one of the 
simplest solutions to these expressions can be had by modeling the perturbations 
as constants. 


In view of the separability of the in- and out-of-plane differential equations, 
it would be both practical and convenient to recast (say) Eq. (II. 12c) as follows: 


(a) For the in-plane components, write: 




_ Tx "I « 2 F3 O'] 

LvJ ^ LO o- 


01 - 
r 


-2»m 2 )t * 


(I*,)* • 


wherein r B = X L co is the angular rate of the reference particle (P), and 

LyJ 

() F.) . would be composed of the perturbing force (s) in-plane components, 

Li j p-C 

only. 


(b) The out-of-plane expression (a scalar) is given simply as: 

j 

or, as an equivalent mathematical operation, 


e z , - ? = - <p 2 <e z , • + e z , • (Jf. ) . 

j 

Before outlining a solution procedure for the above expressions it is 
felt advisable to develop a set of corresponding, linearized differential equa- 
tions, for the motion of Q, as it would be referred to the inertial frame of 
reference. Consequently, this is carried out in the following paragraphs. 


The Relative Equations of Motion (for the Inertial Frame of Reference). Under 
the assumption that a relative position vector (r ) is known (and described) 
in the moving frame of reference, then this information can be "transferred*' 
back to the inertial frame using the idea set down in Eq. (II. 5); that is, by the 
application of: 


(II.12d) 


(II. 12e) 
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R= T(<p”) r. 


(II. 13a)* 


On the other hand, having available the result set down in Eq. (II. 12c) it 
is expedient to begin there and to develop the desired differential equations of 
motion, accordingly. 

For present purposes rewrite Eq. (11.13) as 

r = Tfto-) _1 R; 


or, as an equivalent expression, 

r = T (<p + ) R . 

From this last statement it is seen that the (kinematic) velocity is 

. 

r = (T + ) R + (T + )R, 

* 

where (T + ) and (T + ) are given as Eqs. (A. 6e) and (A. 3a, or A. 7c), respectively. 
Of course this can be simplified by introducing Eq. (A.6e); this leads to 

• 9 

v = -<p (BT + ) R + (T + ) R . 

Consequently, velocity, referred to the "moving frame", is expressed in terms 
of the inertially defined state vector (s). 

Carrying the differentiation another step (forward), the kinematic des- 
cription of acceleration is found to be : 

.. - • .. 

T = <I 2 T+) R - 2<p (BT + ) R + (T + ) R , 

after suitable substitutions have been made. 

Now, when Eq. (II. 14c) is solved for R (the inertial-frame's des- 
cription of the acceleration), the result is: 

♦Recall that R is the representation for r r in the inertial frame, while r is 
its designation in the "local” (rotating) frame. 


(II, 13b) 


(II. 14a) 


(11.14b) 


(n.l4c) 


22 



(tI.14d) 


i R = (T+f 1 r + i p [(T+)" 1 ] R + 2<p [(T+f BT+] R . 

This expression is "symbolic" of the motion; to account for the force systems 
which have (already) been assumed, it would be necessary to introduce Eq, 
(II.12c), for r, into the result above. As an outline on how to proceed (from 
here), the following notes are included, with explanations: 

-I - - - 

(a) When evaluating (T + ) r (=T_r), note that (I^-I^P vanishes 

since Poce , and (L-L)e =0. Also, in the expression for r it is seen that 
. x 2 3 x 

both r and r are present - these are obtained (for present purposes) from 
Eqs. (II. 13b, and II. 14b). Thus, with (TJ replacing (T+) 

(TJ r = i 2 [(T_)(I 2 -I 3 )(T + R) + 3(T_) I 3 P (T+R • P)] 

-2 <p (T.) B 2 [-tp (BT+) R + T+R ] +<TJ £ F , 

} 

= iP 2 {[(T_)(I 2 -I 3 )(T + )]r +3[(T_)P(T+)]r • p} 

+ 2 <p 2 (T.) B 2 B <T+)R - 2(p (T_B 2 T+) R + (T_)£ f . 

j 

Here, it can be shown that (T.B BT + ) = - (T_I T+) = - I and (T_B T + ) = B • 
so 

(T+) R]- 2<^ 2 I 2 R - 2<p B 2 R + <T_) Y Fj • 

j 

(b) Also, in Eq. (II.l4d) there are transformation products which 
must be defined. In this regard it is noted that (T_I T + ) = I and (T_BT + ) 

= B^, as before. 

Introducing all of these results into Eq. (II.14d), it is found that: 


(T_) r - <p 2 [(T_) 


3 0 0 
0 0 0 
0 0-1 
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~ .2 r 

R =9 {(TJ 


3 0 0 
0 0 0 
0 0-1 


(T + ) - 2I 2 } R -JSC UBjET' H (T_) £ F +<p 2 (Ig) R 




or, 


£-i 2 {vj 


3 0 0 

0 0 0 
0 0-1 




Before leaving these expressions, it should be noted that the matrix 


I can be replaced by (T_) I (T + ); hence a contraction of Eq. (II.14e) 
2 2 

follows, as: 



*2 0 0 
0-10 
0 0-1 


(T + )} R + [fTJ If . ] . 
j 


For those cases where the displacements are known (implicitly or ex- 
plicitly), then the linear solution for the motion is (in concept) capable of 
being determined. Of course here it should also be noted that the term M (T_) 
^E\ T ' is a symbolic representation of the perturbative force (s) as it would be 
expressed in the inertial frame of reference. * 


Summary . To this point in these analytical developments, expressions have 
been obtained, with a linearization introduced, which (in theory) can be 
solved for the relative motion of 'Q” with respect to ,r P M . As described, 
these results are for both the moving frame of reference and the correspond- 
ing inertial frame. Both representations spring from the classical (Newtonian) 
equations governing the motion; and both contain symbolic descriptions for 
forces other than the gravitational attraction. It must be recognized that the 
appropriate representation for such forces must be made (in the described 
frame of reference) in order for the equations to be valid. 


♦Any specific description of the forces, referred to this reference frame, 
would modify Eq. (II.lAf) accordingly. 


(II. I4e) 


(II.14f) 
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III. PROBLEM SOLUTIONS 


In the next few paragraphs a method of solution, for an initial value 
problem, will be presented and discussed. For simplicitly and conciseness 
an added simplification will be included, for these solutions, one which will re- 
duce the numbers of manipulations which will need to be carried out. 

A Solution for the Relative Motion. In the next several paragraphs a solution 
for the relative motion of TI Q" with respect to "P" will be described. For 
present purposes this solution is acquired from Eqs„ (II. 12c); these equations 
describe the relative motion (vector r ) in reference to the moving frame of 
reference. Also, since the solution is to be obtained as an "initial value" 
problem the method of Laplace transforms is a natural and convenient mathe- 
matical tool to use. 

Before entering into the mathematics of the solution it is worthwhile 
to study, somewhat critically, the governing differential equation, Eq. (II. 12c). 
In this regard, looking at 

; = ^ 2 [(I 2 -I 3 H; + P)+3I 3 P(?. P)]-2^B 2 ; + I 3 ^F. , 

j 

o 

It is evident that the variables (r, P, r and F) are all defined in the local 
(moving) frame of reference, (~) f . The remaining quantities (including <p) 
are constants; generally, matrices. 

Recalling that P (the position vector for "P") = |P | e f , then it is 
apparent that: 

(a) (i 2 -l 3 )P = 0 ; 

A 

(b) I P - J (see Appendix B, for the definition of J ); 

o X X 

A 

and, (c) r * P = x (the radial relative motion coordinate). 

Introducing other special matrices, described in Appendix B, and as required, 
it is easy to show that Eq. (II. 15a) can be recast as: 


(II. 15a) 
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(II. 15b) 



- -2 

v - ip 


H 3 J x ~ J 3 D r - *p (B 2 ) r + I 3 



j 


{This is the general expression which will be manipulated for a solution, 
r — r(t). ) 


The first operation to be undertaken in this solution is to transform 
Eq. (II. 15b), by direct Laplace transforms*, wherein the only time -dependent 
variables are those noted above. Now, after transforming and collecting 
terms, the resulting expression is found to be: 

[s 2 I 3 + B 2 s - <p 2 (3^-Jg)] ?(s) = [sl 3 + B 2 ]? 0 + V o + 


Y 


wherein s is the transformation variable, and the vectors r t r^ represent 
an initial state of the relative motion. (The special matrixes, J^, in- 
cluded here, are defined in Appendix B. Also, in these operations the 
perturbation acceleration term(s), F, are assumed to be constant (s))» 

In order to acquire a solution to Eq. (II. 16a), the inverse transform 
of the above algebraic equation is needed. However, it must first be solved 
for r(s) explicitly. 


Defining the coefficient of r(s) as a matrix, A(s), then its inverse 
(A _1 ) is needed since a symbolic form of the solution is expressed as: 


: (s) = a' 1 [sl 3 + 2<p b 2 ] ? 0 + A' 1 [y J 


-1 


l 1 /- 


(Note that the matrix, A, has a same form as that described in Appendix B; 
consequently its inverse is of the form shown there. That is; since 


A ~ aJ^+ bjg cJ^ eBg 1 


then 



^ 1 ^ 2 ' eB 2 
ab + e 2 


+ 



♦The transforms used here are the more familiar ones (see most standard 
mathematical texts); the so-called r, p -multiplied" transforms (see texts by 
L. Pipes) are not used in this work. 
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(II. 16a) 


(II. 16b) 

(II. 17a) 
( .17b) 



where the coefficients (a, b, c, e) are identified from Eq. (II. 16a). In this case, 
when identifying the coefficients in A(s), recall that I = J + J + J ; then, 

O 1 Z O 


_ 2 _ -2 . 2 
a - s - 3<p , b = s , 


2*2 • 
c=s +9 , and e = 2 <ps). 


aM7c) 


Introducing Eq. (II, 17c) into Eq. (II. 16b) and collecting terms on powers 
of s, an intermediate result is: 


: (s>=[ 


S % + ^ (I 3 B 2 “ B 2^ + ^ 4I 2 “ 3J 2 “ J 2 B 2 


[ 


2 2 * 2 
s (s + 9 ) 

s 2 I g - s(2<a B 2 ) - s<p 2 j g 


] 


2 , 2 .2 

s (s +9 ) 


]? + [ 


s \~ s (29 B 2 )-39 2 J 2 

3 2 *~2 

s (s +9 ) 




(II. 18a) 


wherein F is understood to mean £ F, (sum of the perturbing specific forces). 

j 

A more meaningful collection of terms, from the above expression, in 
preparation for taking the inverse transform, is as follows: (here, the rationale 
applied is obvious): 


; < s)= (^Ti)Vo 


Vo * (4I 2- 3J 2 )r o- 2 ’ 5B 2 r o + I 3 F 

-r — — t 


s +9 


, 2 * 2 , 

(s +9 ) 


[s (s +9 2 )] 


-{ 


6^ 3 (J 2 B 2 ) +3^ 2 J 2 ? o + 2(5 B 2 F 

— — _ 

[s (s +9 )] 


}- 


3<0 2 J 2 F 

r 2 / 2 . * 2 i 
Ls (s +9 )J 


(II. 18b) 


(In this egression obvious cancellations, etc. have been accounted for; also, 
prior to the next step, the reader may wish to consult Appendix B regarding the 
special matrixes used herein). 


Recognizing that the solution (here) is to be expressed in terms of ele- 
mentary functions; and, knowing the inverse transformation equivalents for 
these, it is prudent to "reduce" some of the algebraic functions above. Using 
partial fractions, the following equivalents are established: 
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( 1 ) 


( 2 ) 


and, (3) 


1 

J_ /I 


s 

'i . 

, 2 '2, 
s (s +<P ) 

• 2 Is 
<P 


2 *2 
s + <p 

/ • 

1 

, -i- f 

i 

1 


2 . 2 . *2 
s (s + p ) 

•2 \ 
<P 

2 

s 

" s 2 + 

a 2 

<P 

1 

- i r 

1 

_Lf 

1 _ 

3.2 *2 
s (s +<p ) 

•2 L 

<p 

2 

s 

= 2 \ 
<P 

s 


2 , 
s + ip 


*)]■ 


(IL 19)* 


Relative Motion Displacements. Making use of the expressions above, transforming 
Eq. (II. 18b) and clearing terms, a solution for r(t) may be expressed as: 

r =, 3J F 

— + I — + — — I - 1 (J B ) r +3J r 

0 


T Z, od X •• • 

;«-[(« -3J )? -2B 2 + I 3 — + — 

<P <p (ft 


+ 2B. 






3J 2 F 2 r - r o 

t + (I 0 "4I rt +3J ) r +2B (I c 

2 L v 3 2 2 o 2* 3 

<P 


2B F 


+ 


3J 2^ "*~2 J COS ^ t) + L 6 ^ J 2 B 2^ r o + ^ I 3 +3J 2 ) ~^ + ^2 J sin (II. 20a) 

<p <p<p 


A more uniform representation of this result can be had by making use 
of the matrix equivalents mentioned in Appendix B. Thus, after some manipu- 
lation, Eq* (TI. 20a) can be recast as: 

?(t) = [(« + J > r ' 2B — + (J + 4J 2 + J 3 ) -^] 

(p <p 

r - r o FI* T3 FI • ,2 

- [ 6J 2 B 2 r o + 3J 2 T + 2B 2 71 J "Li J 2 72 J 


<P 


<P 


<P 


J* “ 

+ [(Jg - 3J X ) r o + 2B 2 T 2 - ( J ! +4j 2 +J 3 ) Ti J 008 (it) 

(equation continued on next page) ^ ^ 

*Some inverse transforms are noted here; (a/s) (constant); ^ (l/s ) - t 
(the independent variable); s£“^(2/s ) = t ; ^ l<p / (s + <p )] ~ sin (<p t); 
i -1 [s/(s 2 + <p 2 )] = cos (<p t); and, s£ -1 [sr(s) - r(0)] = r (t)* 


28 


+ [ 6J 2 Vo + (J 1 +4J 2 +J 3 ) “ +2B 2 -k ] Sin W' 

<P <P 

Eq. (II. 20b) describes the relative motion of f, Q" with respect to T, P”, 
in ’local coordinates” (x, y, z), for the case where a central attraction and 
fixed perturbative forces (V F.) are considered. 

An expansion of the various matrices, etc. is instructive at this point 
to illustrate the composition of the various terms involved. For convenience 
in the resulting expansions, the specific force terms have been separated from 
those involving the initial state vectors. Thus: 
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0 
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*Some of the quantities included here are in preparation for the non-dimensional 
zation to be introduced subsequently. 


(II. 20b)* 
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Here, the state vectors: 

o 

fji ^ „ ■ * * 

i(t) * (x, y, z) , r o s (x o , y Q> z o ) , r 0 H < x 0 > V z 0 > ■ are < each > 

column vectors, as described. Incidentally, this result (above) may be checked 
with that given as Eqs. (A. 29b), in Reference Cl]** 


The Relative Motion Velocity. Equations ([1.20) provide analytic expressions for 
approximating to the relative motion displacements of ’Q" with respect to "P". 
However, the corresponding relative velocity is needed if one is to complete a 
description for the state of motion. 


The simplest, and most direct, means for acquiring this latter informa- 
tion would be to differentiate the displacements given above. When this procedure 
has been carried out, on (say) Eq. (11.20), it is found that: 


= - (3J [2B rt r +1 
0 l 2 L 2 o 2 
<P 



+ 2B 

2 







-{<V 3 V ? o + 2B 2 • - (J 1 +4J 2 +J 3> 7£ } sin 

<f> ; JP 

+ K Vo + (J 1 +4J 2 +J 3> ~ + 2B 2 Til ° 0S ^ 

<P <P 


For clarity, and continuity, an expanded form of Eq. (H.20d) is presented 
below with terms involving F separated from those described in terms of the 
initial state. (This is a companion expression to that given as Eq. (II. 20c)). 


^Reference [l]: "Relative Motion of Orbiting Satellites’’, byJ.B. Eades, Jr. 
NASA CR 112113, dated July 1972. 
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r(t) / x v z \T - T 

Here, — ^ =(—,—,—) » r - (x , y , z ) , etc. By vlrture of the fact 

<P <P <P <P 

that this result is akin to that in Eq. (II, 20c), then these scalar expressions 
also may be checked against corresponding ones in Reference [l] . 


(IL 20e) 


Non-dimensionalization . To this point in the development the state equations 
have been expressed in terms of dimensional quantities, see Eqs. (11,20). For 
many operations, and situations, it would be useful to have these expressed 
in a dimensionless form. 

For purposes of compatibility with past work, and for the convenience 
of notation, etc. , non-dimensionalization will be carried out as follows: 
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(a) All distances, such as r(t) (= r(x, y, z), and the initial value 
r(0», are normalized using the reference particle*s radius, |p| . 

(b) The velocity for , 'Q M } and the associated initial value, will be non- 
dimens ionalized by division with |p[ <p (the speed of the "P M particle). 

(c) To retain compatibility, the perturbing forces are normalized by 

o g 

the quantity |p| <p (the centrifugal force acting at the position of "P"), 

(d) Also, since Ip is a fixed parameter (here), the time derivatives 
will be replaced by position derivatives; this is accomplished by writing, 



A 

dt 


A 

d ip ° 


In explanation of this procedure the following variables (and associated 
initial values) are introduced: 


'i ' 


X 


'V' 


X 


'V 


‘ F ‘ 

X 

V 

1 



A. f 

1 

• 



1 

T? 

Ipl 

y 

> 

V 

|p|<p 

y 

» 

T n 

1 

to 

r 

y 

_c _ 


z 


£' 

m — 


z 


. T c. 


F 

L zJ 


d ^ 1 d 

where — = (~) = — — ; here A (p), symbolically represents the dimensionless 
^ . ■ <P d 

relative position vector ; corresponding to this, A *(<£), will be used to describe 
the associated dimensionless velocity vector. 


( 11 . 21 ) 


State Equations in Dimensionless Forms. Making use of the definitions introduced 
above one can easily recast the solution equations into appropriate dimensionless 
forms. For instance, Eqs. (II. 20b, II.20d) maybe written as: 

* «P> = [(«! + J 2 ) * 0 - 2B 2 I - q + Pl + 4J 2 + J 3 ) f ] 


-[ 3 J 2 ( 2 B 2 A 'o + V + 2 B a f ]* , -[| 

(equation continued on next page) 
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+ [ (J 3 - 3J l )J o + 2 V 'o - (J X + 4J 2 +J 3 )? ] C ° S * 

+ [ 6J 2 B 2- ' k o + (J 1 + 4J 2 + W + 2B 2 ? ] SiQ ^ ’ 

and 

X’ W --[8J 2 (2B 2 I o + I 2 I' o ) + 2B 2 f] 

- 3 [v> - [ (J 3 - 3J l U ‘o + 2 V ’o ~ (j l + 4J 2 + V Stn ^ 

+ [ 6J 2 B 2 Jl "o + (J 1 + 4J 2 + J 3 ) *o + 2 V] C0S * * 

Due to the simplier notation introduced here, these expressions will be 
the ones manipulated hereafter (in these developments). Needless to say, the 
dimensional form(s) of these quantities may be recaptured by an application of 
the transforms shown in Eqs. <11.21). 

Prior to describing the state of motion, in terms of the inertial relative 
motion variables (R, R), a more compact and revealing form of the position ex- 
pression, above, will be developed. This particular representation combines 
the terms in Eq. (II. 22a), making a geometric description of the relative motion 
traces easier to "see". This development is outlined in the following paragraphs. 

Rearranging the State Equations . In this rearrangement of the state equations 
(see Eqs. (11.22), above) it is expedient to commence with the dimensionless dis- 
placement expression (A); and, to acquire the corresponding "velocity” ex- 
pression by direct mathematical means. Also, since the out-of -plane displace- 
ments are uncoupled from the in-plane ones, then it is simplier to work with the 
corresponding separated expressions. Consequently, Eq. (FI. 22a) are separated 
as seen below; see Appendix R for descriptions of the special matrices used here: 


(II. 22a) 


(II. 22b) 
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A «P)F v+v> =[(«! + J 2 )i 0 - 2B 2 I; + (J x + « 2 > f] 

- [ 3J 2 (2B 2 l o +I y + 2 V ] ” - [ I j 2 f ] ” 2 

+ [' 3J /o + 2 Vo ' (J 1 + 4 V] C0S V 

j + [ 6J 2 B 2 i o + (J 1 + 4J 2 )J o + 2B 2 ? ] Si “ * 

j + jjg + U q - f) cos <p +A ^ sin (p J j- . 

The first terms to be considered (here) are those involving the circular 
trigonometric functions. It is to be demonstrated that these can be related to 
the transform matrices, T (<p ) and T(cp + ) — but when they have been reduced to 
a second order form. Thus for the defined matrices (T) 2 , write; 


T (<p~) = I 2 cos (p + B 2 sin <p + J g = (T _> 2 + J g , 


and 


T(<p + ) = I 2 cos ip - B 2 sin cp + J g ■ (T + > 2 + J g . 

Now, in Eq. (11,23) define coefficients (K , K and/or K f ) as shown below: 

c s s 

I 2 K c cos (p = {“3J 1 j4* o + 2B 2 A' q - +4J 2 ) T } cos (p , 

and 

r 2 K s Sin<P 5 { 6J 2 Vo + <J 1 +4J 2 U ’ o + 2B 2 Sto 

A manipulation of the "sine" expression is needed here; in particular the 
,f B " operator is removed (see Appendix B for the equivalent operations used 
below). For this purpose the following operations are described: 

Since 

I 2 K s Sin H 6 (B 2 J Vo + < J l +4J 2 ) <- B 2 B 2 )I o + 2B 2 f } Sln<P 
then the right side is manipulated so that it becomes: 


(11.23) 


(11.24) 


(II. 25a) 
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(II. 25b) 


{b 2 ( 6Jl I o + 2 T) - [4J 2 + J x ] B 2 (B 2 I' o ) } Sin p 

= B 2 K J o - [ 4J 1 + J 2 ] Vo +2? } Sin<P; 

or, this is written as, 

B 2 K ' sin V - B 2 {6J l A o - (4J 1 + J g ) B^ > 0 + 2 T } sin <p. 
Next, let it be assumed that: 


IK cosc? + B K 1 sin^D = (Ti A +(T + > A , 

2 c ^ 2 s 2 a ^ 2 s’ 

(where A and A are two matrices to be defined below). The present problem 
s a 

Is to determine the quantities (A , A ). By expanding the matrices ,f T ”, using 

a s £ 

Eq. (11.24) above, 

I K cos (p + B K r sin<p = I (A +A ) cos <p + B (A -A ) sine? ; 

2 c ^ 2 s 2 a s ^ 2 a s 

then, by matching coefficients, it is found that: 


K 

- A + A 


K 

+ K ! = 2A 

c 

a s | 

j 

\ , or, 

c 

S i 

K* 

> 

> 


K 

- K f = 2A 

s 

a s * 


c 

S 1 


Substituting from Eqs. (II. 25b), it can be shown that 

A a " I {Vo + [ J 2 - 2J 1 ] V o + [ J 1 ~ 2J 2 J ? } > 


A s = f f 3 V o + K + J 2 ] B 2 K - [ J 1 + 2J 2 ] ? } • 

It is quite apparent that these two matrices have a marked similarity; 
they differ by a constant multiplier, and by a change in the sign of the matrix, 


♦Subscripts have been added to the matrix designation "A", indicating that they 

are acquired by ,f addition ,r and by "subtraction” of the constant matrices, K 

and K f . c 

s 


(n.25c)* 


(II. 26a) 
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j (Note that the matrix operator altering the sign of while retaining the 


sign of J is 


-1 

0 


"( 1 - 2 ) 

0 

0 

1 


0 

■ 

1 . 


= (J 2 -V = < I 2- 2J l> i 


hence 


Therefore, 


A s =3 ( J 2- J l )A a-> 


(n. 26 c) 


(I 2 K c ) cos <p + (B 2 K^) Simp = 

= [^ T - > 2 +3<T+ ) 2 < J 2 “ J lO A a * 

This expression, Eq. (II. 27a), describes the contribution tojfc (p) provided 
by the trigonometric terms. Hence, this partial solution is expressed as: 

l A * Htrig " [ (T ->2 +3 < T + >2 (J 2 - J l > ] V 
or, since (J - J ) is its own transpose; 

a 1 

[ AA Htrig = [ I 2 +3 ( J 2- J l ) ] fr ->2 A a> 

where A is defined above, in Eq. (II. 26a). 
a 

Next, consider the constant and first order secular terms In Eq. (11.23). 
First, however, a manipulation of the secular term is needed. For this, write 

9 t 6J 2 Vo + 3 V'o + 2B 2 f ] = B 2 [ 6J 1^ 0 " 3J 1 (B 2 ^o ) + 2f ] * • 

"■ 3 

Now, when the constant term in [i Jv((p)l is multiplied by M — J ”, 

£ X 


01.27a) 


01. 27b) 
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I j iK +j 2>i 0 - 2 <v"’o> + (j i +4j 2 )f Hv"o - w'v + fv ; 

and, therefore, using these expressions, the two quantities can be combined as 
follows : 

(4J X + J 2 « 0 - 2(B 2 J< o ) + B 2 f [(4J X + J 2 )i 0 - 2B 2 i-J<p} + (J 1+ 4J 2 )T-2B 2 ^ 

= 02 - I Vl* 1 > [< 4J 1 + " 2 Vo ] + K + 4J 2> “ ^2 J ? ' 

Adding to this expression the remaining secular term, then the last part 
of the in-plane partial solution, for I h (<p), can be joined together as: 

A *tP) = P 2 -f^ B 2 V[ (4J l +J 2 )£ o ~ 2 V'o] + { (J l + 4J 2 ) -’ 5 [ 2B 2 + | J 2 !fl ]} f - 

Finally, the zth component of this (positional) relative motion solution is 
seen to be: 

J 3 J((f>) = ^) cos< P + ^o sin ^} ■ 

Now, collecting the various partial solutions, then, Eq. (II, 23) can be 
symbolically represented by: 

A te» = [l 2 + 3 «J - J x )] (T.) 2 A a + [l 2 - f <PB 2 J x ] K o + * r T + J 3 {f + U 0 - T) cos 0 


+A 


1 sin <p \ ; 
o J 


wherein, 


A a " K 3 Vo + [ J 2 - 2J 1 ] Vo + [ J 1 - 2J 2 ] ~ T ) * 


K o S(4J 1 +J 2 U o' 2 V’o> 


(II. 27c) 


(II, 27d) 


(11.28) 
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and 


V^V^K+fV 5 ]- 

Note that the terms A and K , above, are constants — they depend 

a o 

mainly on the initial values; however, the last term, (the coefficient of T), 
is secular in its composition. It is important to keep this in mind since the 
relative velocity will be determined (next) by differentiation. For this, 
bears a derivative. 

The Relative Velocity Equation. Rather than follow a procedure analogous to 
that above, for the description of relative velocity, a more direct approach is 
taken. That is, by differentiation, Eq. (11.28) yields the corresponding velocity 
equation. 

( The result of differentiating the displacement can be written, generally, as: 

A'(V) ■[l 2 + 3(I 2 -J 1 )]fT.)' 2 A a -|B 2 J 1 K o + V T f + J 3 {t' o eos*- . 

• O 

Of course, as shown in Appendix A, (T_) =pB (T_) ; and hence, (T_)* = 

2 2 2 2 

B (T ) . Also, from the definition of '3/ it is evident that: 

2 -'2 T 

n =-[ 2B 2 + 3J 2 <p]. 

Therefore, 

V(<P) = [ I 2 +3(J 2 ' VX (T - ) 2 A a'2 (B 2 J 1 )K o + *T ?+J 3 {* o COS< ^" <V ? )sln<P } ’ (n - 29) 

with as shown above. 

r 

Summary . Eqs. (11.28, n.29) describe the relative motion of a test particle , 'Q" 
with respect to a reference particle ”P M in terms of variables which are defined 
in the local (rotating) frame of reference. Here, the initial state parameters 
(4 , A f ) Q are dimensionless quantities as defined previously. Since this problem 
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has been solved as an initial value situation, then the coefficients (A , K ) are 

a o 

expressed in terms of the initial state parameters. Also, included here is a 
representation for a "fixed value" perturbing force system (i.e., r, the dimen- 
sionless specific force vector). 

The expressions noted above are merely rearrangements of the (similar) 
equations developed earlier (see Eqs. 11.22). The purpose in presenting this last 
set is that they provide a simpler and more direct means of representing the re- 
lative motion state, geometrically. 

In subsequent paragraphs a description of these geometries will be given; 
however, prior to this, a parallel development for the inertial description of these 
motions is to be obtained. 


NOTE: Before moving on to these next developments it would be useful 
to write the state equations (II . 28, 11.29) in their expanded scalar format. For 
this, proceed as follows: 

(a) First, from Eq. (n. 28), the in-plane position variables are de- 
veloped from : 

i 2 i «p) =[v 3 < j 2 - j i>] f^vEv! * b 2 j iX + V- 


In matrix form this can be written as: 


i 2 -fc (<p) = 


-cos (p sinp 
+2sin<p 2cos < p 


+3 0 

0 0 




+2 

0 


k' + 
o 


1 

0 



1 

o 

i—4 

J 

r 

o 

1 


0 2 

-1 

| 

f 

r 

o 

tH 


o 

to 

9 

rH 

! i 

1 

0 1 


.- 2 

*i 

>+ \ 

f 

l 

0 4 

— — 

+ (p 

1 

o 

<M 

L, 


0 






° % 


>T 


(II. 30a) 


(II. 30b) 
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This equation, when separated into its scalar form, gives: 

£ = (_3 £ ^ cos (p + ^ sirup - 27?^ cos <p- cos <p - 2T^sincp) + (4^ + 27?^) + (T ^ + 2cpr^) 

and 

7?^) = (6^ s in<p + 2^ cos^+477^sin<p+2r^sm<^-4r^cos0) - (6£ o <p- 7^ + 24 ^ + 3^7j£) 

-[ 2 ^- (4 -!^ t „]- 

Finally, the zth-equation is seen to be: 


CfiP)“r c (1 -cos<p) + C 0 c °sp + C^sin^. 

Next, the velocity equation (s) will be expanded, and presented. 

Beginning with Eq„ (IL29), and separating the in- and out-of -plane com- 
ponents, then the in-plane resultant is: 


i/»=[ 


I +3(J 
2 2 


-v] 


B rt (T )«A -;{BJ)K + ^T 
2 v “ ; 2 a 2 2 1 ; o T 


this can be set down in matrix form as : 


* 

sin <p cos<p 


" 

3 0 

A + 

m -* 1 

0 2 

A r + 

1 0 

-1 

T 


0 0 

to 2cosp -2 sin<p_ 


0 0 
U J 

o 

_1 0 _ 

0 

0 -2 
w tm 

J 

2 

- 1 °. 


f 

’4 0 ' 


o 

to 

1 

_ f 

“o 2 " 


o 

o 




I 0 + 




- 


- 

l 

iH 

O 

1 


-2 0 

1 l 

-2 0 


J 

© 

1 

4 


(It is recognized that the parameters A , Ji f , and f (here) are necessarily 

o o 

composed as the in-plane vectors only). 


Next, when these matrix expressions are expanded into a scalar format 
one has: 


(H.30c) 


(II o 30d) 


(II. 31a) 


(II. 31b) 
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ai.3ic) 


£ r (P) = (3£ 0 sin^+^cos<p + 2T)^sm<p+T^sln<^) -2 t^(cos^3-1), 
and 

V'(<P) = (6£ cosp -2£ , o Bin<p+iri' o COs<f)+2T£COS<{>+4T r jSin<p)-(6£ o +3ri' o )-(2T£+3<pT' 
Lastly, the zth velocity equation easily yields: 

^ t (<p) = £^cosp“(C o -r^)sin^ . 

The companion dimensional expressions, to those shown above, are 
readily deduced from Eqs. (II. 20c, and II.20e) presented earlier. Of course, one 
major difference between these two resultants is the reformation leading to the 
geometric interpretation provided by the latter expressions as compared to the 
earlier ones. 

The development procedures for a relative motion, referred to a local 
frame of reference, are completed with this summary section. Following this, 
a similar description of this relative motion, referred to an inert ially oriented 
frame, will be developed. This appears in the next few paragraphs. 

The Relative Motion, in Inertial Coordinates . In determining the desired state 
equations, here, those which have been described (above) will be modified — 
i.e. transformed — according to results obtained earlier in the section(s) on 
kinematics. Accordingly, using Eq. dl. 13a), the relative position vector, re- 
ferred to the inertial frame, is related to the corresponding vector, described 
in the local rotating frame, by: 

R =.T(<p") r. 

From this single relationship it appears that the transformation could 
be carried out quite simply (in concept). However, it must be remembered 
that the vector (r) has been determined as an initial value problem; therefore, 
it will be necessary to transform the initial state vectors, also, in this develop- 
ment. This task will be undertaken first, below. 


). 

(II. 31d) 


(II. 3le) 
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(Initial Values) . According to Eqs. (n . 13b), II .14a, and A. 7), the variables 
of interest here are described, and related, as follows: 


r = T«0 + )R, 
r = T(<p + )R + TftO + )R, 

T(<p + ) = I cosco-B^sin <p + J^. 

Now, at the initial position (t = 0) these expressions are specialized to read: 

r(0) = ? o = T(0) R o , 

* 

and ?(0 )^r o = T(0)R o + T(0)R o . 

Herein, T(0) = I + J 3 = ; also, by definition, T(0) = These results 

lead directly to: 


Since 

then 

wherein 


and 


? o = I 3^o‘^ B 2 I 3^o' 


When these expressions are separated into in-plane and out -of-plane vectors, 
one obtains: 


(in-plane) 

and, 

(out-of -plane) 


Lr = I R , I r = I R - <pB 0 R ; 
2 o 2o 2o 2 o 2o 


Jr - J n R , J_r ~ J„R • 
3 o 3 o 3o 3 o 


The corresponding dimensionless forms of these results are noted to be: 


(in-plane) 

and, 

(out-of-plane) 


r Ji. « L ft f I * Lft» - B ft ; 
2*o 2 o 2 o 2 o 2 o 


J £ = J ft , = J 0 R* . 

3 o 3 o' 3 o 3 o 


♦Here ft and ft 1 are the dimensionless (state) variables, referred to Inertial 
axes; ft s ft (H ,H, Z) and ft 1 s (H», H', Z'), analogous to Eqs. (11.21). 


(II. 32a) 


(D. 32b) 


(11.32c) 


(II,32d) 


<R.32e)* 
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(Equation Transformations). For the transformation from "local-coordinates " 
to "inertial" ones, the relationship to be applied is given in Eq. (II, 13a); namely, 

wherein, generally, 

T(<p~) = I cos<p + B sincp + J „ 

a Z o 

For the actual use of Eq. (II. 33a), above, the expression for r(<p), 
appearing as Eq. (11.28), will be employed. Consequently, in dimensionless 
format: 

m •r.mn * «p> - t 2 #>-) [i 2 +s«i 2 - j i>] t 2 *d A a+ T 2 «,-j [i 2 - § *> b 2 jJk o 

+ T^ftT) + J3 |t + (ft^-T) cos (p + sin pj- . 

Before proceeding, a closer look at the lead term on the right of the 
equals sign is in order. 

Expand the matrices shown there, as follows: 

(T ->2 [ J 2 + 3(J 2 - J l>] (T ->2 = < T -> 2 y TJ 2 + 3(T ->2 (J 2 - J l )(T ->2 = < T ->2* < T ->2 
+ 3 (TJ 2 (T+) 2 (Jj-Jj). 

since the matrix (J g - J ) is its own transpose. Thus, in view of the multiplica- 
tion property noted in Eq. (A. 7e); 

( t _) 2 [i 2 + 3(J 2 - J x )] (T_) 2 = T 2 (2o~) + 3(J 2 -J x ) 

Next, making the required substitutions, from Eq. (II.32e), it is easy 
to show that* : 


♦Subscripts (^)j are added to indicate that these matrices (now) infer a repre- 
sentation in the inertial frame. 


([1. 33a) 


(II. 33b) 


(II. 34a) 
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|[i 2 »oMJ 2 -2Ji )B 2 R^<Jr 2J 2 )r] 

(2J -JJK -2B^ ? , 

“ 1 z o z o 


and 

* = V<V 4 V-^ 2B 2 + f v]- 

Introducing Eq. (II. 34a) into Eq. (II* 33b), and noting that the appropriate 
forms for A , K , # are those in Eqs. (II. 34b), then 


&(<P) = [' 


T 2 (2<p-) + 3(J 2 -J x )] 


A + T 
a 2 
I 




K + T (<p”) ^ r 

Oj 2 ' T 



- T) cos (D + ft 1 sin 
o 



Eq. (II. 35) is the dimensionless, relative motion position vector referred to the 
inertial frame of reference. Note that T is the same vector as described before. 


The Inert ially Described Relative Velocity . Following the procedure used for 
the development of Eq. (11.29), then the velocity vector, here, referred to an 
Inertial frame of reference, will be obtained by differentiation. (It should be 
apparent that this is equivalent to the differentiation of Eq. (H. 13a), and the 
subsequent substitution for r and r (or their equivalent dimensionless form (s)). 


Taking the direct mathematical approach, that of acquiring the velocity 
(ft*), from the position vector (ft), then (symbolically); 

*’ #»=2 [ T 2 '< 2 ^>] A a f [V^ViK 

I I I 

+ T r (</T)* T + T faf) p£ r ]f + J„ (8* cos<p -(ft -T)sin<p}-. 

2 T 2 L T J 3 l o o J 

For this expression it should be recognized that (see Appendix A), 

T'(<pl = BgTg (<£»”), and T^(2<p") = . 


(II. 34b) 


(11.35) 


(H. 36) 


(II. 37a) 
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(II .37b) 


Also, 'I' , being the secular coefficient of T, bears a derivative; namely. 

When these various quantities are inserted into Eq. (II .36), and some 
manipulations are carried out, it is found that (one form of) the velocity vector 
is: 

- B 2 {[ 1 2^>] + [ T 2 (2<P_) ' 3 < J 2- J X ) ] A a I ' T 2 (,p ' ) [f + 2 V>]} 

- 3T 2 (^ _ ) [jg^] f + J cos ( P ~ (® G ~ ?) sin . (II . 38) 

In this expression the vector ft = ft(p.) is the relative position vector given 

by Eq. (11.35); also the coefficients (A , K ) are those quantities noted as 

a l °I 

Eqs. (II. 34b). 

Summary . The expressions for ft and ft 1 , given by Eqs. (II .35, and 11.38), re- 
present a linearized solution for the stated relative motion problem when it is re- 
ferred to an Inert tally oriented frame of reference. These results describe the 
same problem as that in Eqs. (11.28 and 11.29); the difference between these two 
solutions is the frame of reference implied by each. As in the earlier solution, the 
external force system, here, has the same fixed components as before; that is, each 
component is parallel to one of the axes for the ’’local rotating frame of reference”. 

There is a second exception for these results. Here, the parameter co- 
efficients (A , K ) have been altered (necessarily) to describe the proper set of 
a o 

initial values , those required for the inertial frame’s representation. 

In order to see the more definitive representation for the inertial frame's 
results, the state equations are rewritten in the expanded format shown below. 

That is, from Eq. (II .35) write: 
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(a) for the in-plane displacement solution; 
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and; (b), for the out-of-plane displacement; 
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In place of the relative velocity (Eq. (11.38)), a more convenient ex- 
pression, for manipulation, is: 

B'(p)=B 2 {[2T 2 (2<p-)]A ai + T 2 (^-)[(l 2 -|pB 2 J l )K 0i+ 'l' T r]} + T 2 ( ( p-)[^f 

* I (B 2 J 1 ) K Oj] + J 3 h C0S * ' ( ®o ~ ' T) SiD V) ■ 

Written in matrix notation this becomes; 

(a) for the in-plane solution: 
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(continued on next page) 


(II. 39a) 


(II. 39b) 


(II. 40a) 
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and, (b), the out-of -plane solution: 

cos <p 

Equations (11.39 and 11.40) are a more descriptive set of expressions for 
this solution to the relative motion problem. Here, as previously, the displace- 
ments and speeds depend on the initial state quantities (S . R 1 ) and on' the applied 

oo 

(fixed) specific force (T). It is seen, and expected, that these force components 
are the same as those used earlier since these equations describe the same re- 
lative motion problem but as it would be seen in the inertial frame of reference. 

With one aim of this investigation being that of describing a geometry for 
these motion traces, that task will be undertaken in the following sections. For 
the illustrations to follow, both solution descriptions - i.e. , the one represented 
in the "local rotating" frame of reference, and the other referred to an inertial 
triad - will be examined. Subsequently, the geometric figures which these describe 
on the various coordinate planes will be described graphically and analytically. 
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IV. GRAPHICAL DESCRIPTIONS 


Geometric Representations . In the following descriptions certain relative motion 
traces will be examined. Those which are selected have been chosen so that com- 
parisons may be made between representations found for the two reference frames. 

The first geometries discussed will be those for in -plane displacements. 
For convenience and conciseness all discussions and graphical representations 
are cast in terms of dimensionless variables. 


Only a limited number of problem situations are described, in detail, below; 
however, a compendium of these results is included herein for reference purposes. 
The details for constructing the various trace geometries is contained in another 
report*. The first case study is that for: 


(a) In-Plane Displacements; Rotating Triad . The general expression to 
be investigated here was previously given in Eq. (n .30a); that result was: 


r 2 * < = = [h +3 (J 2 - J X> ] T 2 «*"> A a + k " I* Vl l \ + * r 


(11.41) 


From this equation one sees the displacements expressed as three vectors ; 

one, each, to describe contributions from the parameters A , K and T, re- 

a o 

spectively**. In order to describe a geometry which these produce, and to do so 
in a systematic fashion, each vector statement will be examined separately; and, 
then, in combination with the others. 

The first component vector (that due to A ) is noted below. It is recalled 

£L 

that A s A ( A , A f , T); however, at present it is more convenient to let 
St si o o 

A a be expressed as a general two component vector; i.e. , A^= A^ (A^, A ^), 

Then this contribution to A (<p) is described as: 


* "Construction of Relative Motion Traces" (an Interim Report), by J. B. Eades, 
Jr. , AMA Report No. 73-39, August 1973. 

**Thls composition is "mathematical" in its make-up; a more "realistic" con- 
struction will be obtained in the following section. There the problem is recast, 
to separate initial values from the effects of an "applied force". 
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( 11 . 42 ) 


r - i _ 

' 64 ' 


-cosp sinp 


'V 

Lm #»J, " 


= 




4 -JJ. 

67} 

1 

2sin<p 2cos<p 


• V 


where these two scalars (A.), in terms of the initial values, are: 


and 


A = 34 + 27}' + T , 
1 O O 4 


A 2 = _ 2T r} * 


It is evident that the expression above describes a (2:1) ellipse on the 
(4* 7?)-plane. This ellipse is centered at the coordinate origin and described by 
the parametric equation 


( V>‘ 




= 1 . 


A general sketch of the trace is shown on Fig. II. 4 below*. 


The second vector contribution, arising from the K q term, involves 

only initial value quantities; i.e,, K S K (A , a 1 )° However, for convenience 

0 0*00 

this partial solution is expressed as: 





, wherein 


K =2(2£ + f?*) 

1 o o 

K 2 SW o- 2 ^ 


The parametric equation for this trace is readily determined as: 




The equation describes a line on the (4, V) plane which is parallel to the 77 -axis. 
Note that the line originates from 677 (0) = K n and that the expression is linear 
in the independent variable <p (-pt). 


*This is a familiar result for a relative motion described on this plane; see 
Reference [l] . Also, note the direction of motion about this trace geometry. 


( 11 . 43 ) 


( 11 . 44 ) 


(II. 45a) 


(II. 45b) ; 
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V 

FIG. II. 4. Sketch showing an in-plane trace produced by the coefficient, A a . The figure is described 
in the local, rotating frame of reference. Note that the trace’s initial point (0), and 
selected angle (<p) positions as indicated. (Refer to Eqs. (II .42) for a mathematical des- 
cription of this figure). 



The secular influence for the in-plane trace is eliminated when vanishes* 

This condition could be achieved w r ith either of the starting conditions: 4 = ? 7 , = 0, 

o o 

ox* 73 * --24 . (Such constraints on the initial state parameters may occur in either 
l o o 

an ’’intentional” or ’’chance” situation). 


The third component vector is expressed in terms of the applied force (T), 


where T = T (T^ , T , T^). 

The trace for this part, of the solution is obtained from: 



1 (2<p) 

tm 4 - 1 <p 2 


r i 

T 

l r} -> 


(II. 46,a) 


A descriptive equation for the geometry is: 


6 77 - 4 T 1 
3 r}i 





which describes a parabola on the (4, V) plane; but, with its apex not at the 
origin. The figure opens in the negative 77 -direct ion*. 


(II. 46b) 


The full in-plane trace for this problem is constructed by adding the three 
partial solutions above. It is apparent that the general figure would be an ’’ellipse” 
which moves away from the origin at a ’’linear -plus -parabolic rate”. 

Instead of examining these various figures in detail, now, a physically 
more realistic study of this problem will be undertaken. This considers, as 
separate entities, those influences which arise as a consequence of the initial . 
values (X Q t and those due to the applied force system (T). The individual 

effects are described below. 


The influence of initial values (A q , Jt^) on the in-plane trace is obtained 

from A and K , collectively. In matrix format this partial solution is: 
a 0 


*This description assumes all r-components are positive valued. 
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-eos<p sinV 3 


(11.47) 


A [ I 2 jt >)]x S 


Afi 
A r) 


2sin<p 2 cos <p 


i r . n 


2 J 


1 0 
3$2 , 


L 2 


K_ 




wherein, now, 


and 


Ar^ + ^' K fV A 2 = «o ; 


Kj =2(2i o + T]' o )- 




as before. 


The geometric figure describing the equation above is a f *wandering 
ellipse"; one whose center is continually moving from the coordinate origin*. 
(This is typical of initial value problems; see Reference Cl] ). A quadric equa- 
tion for the trace is: 


"2 T ~2 

+ A 2 


[(A^-ty + f Kjp] 2 


MW) 


= 1 , 


Note that when the secular term is eliminated (i. e. , K s 0) the figure becomes 
a closed locus. It is the (2:1) ellipse with its center not at the coordinate origin. 
An appropriate cartesian equation for this figure is: 




«<? + «J 2 > 


<A 1 ? j-K 2 > 

< + «; 2 » 


=i. 


(See Fig. II. 5 for the traces). 

The effects of an applied specific force (T) produces an in-plane trace all 

its own. For this study the initial state values are zeroed (A - X 1 =0); con- 

o o 

sequently, this partial solution is acquired from terms in A and from the f 

£L 

solution noted earlier. The equation for this case is found to be:. 

*A more definitive name for this figure would be cycloid. 

**Due to the linearization which has been imposed for these solutions, the closed 
figures have a periodicity which matches that of the circular (base) orbit. It 
must be remembered that this is an approximation to the more realistic case. 


(II. 48a) 


(IX. 48b)** 
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-.002 0 .002 .006 


Sketch (b) 



FIG. II. 5. Graphs to illustrate the influence of initial values on the in -plane 
relative motion displacements, as referred to a rotating frame of reference. 
Sketch (a) describes the more general case (see Eqs, (II. 48a)), wherein both 
and r£ are present. Sketch (b) illustrates the consequence of removing 
the secular influence (i.e. , 77 ^ =0). The two ellipses (sketch (b)) describe 
the influence of u (curve I) and a ^ 11 (curve II) on the motion traces. 
Note : All curves here originate from the origin {A = 0). 
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A ftP ) J < 
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A Xf 
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-cos(P sin<p 
2sin<p 2cosp 
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-2 r 
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L n J 
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(29) 

(-2<0> (4- |<P 2 ) 
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i 

1 


This expression suggests: (I), a (2:1) ellipse; couplied with (2), the para- 
bolic figure found for the r)- terms. (Note that this geometry is not a closed 
figure; and that there are multiple secular Influences present). 


It is seen here that when T = 0, the degree of divergence diminishes, 

7) 

This constraint leads to a cycloid for the motion's trace; or: 


and 


A 2 4 = T £ (1 - cos <p), 


At? 

— — * “ (<p - sin 9). 


On the other hand, when T h vanishes the figure is akin to both the cycloid and a 

€ 

parabola. The describing equations for this case are: 


and 


^ 2 ^ 

— — - T (<p - sin <p), 
2 V 




(A general quadric for this particular in-plane trace is: 

,2 


<V>‘ 


r / + 4T n 

€ n 


(A + 9 2 ) 2 

+ — 2 2 [(1 -cosp) 2 + {<p- sin(P)^] ) . 

4 <T fi +41^ ) 


(b) In-Plane Displacements; Referred to an Inertial Triad . In this section 
the motion’s displacement traces are referred to an inertially-oriented triad. 

The general analytical expression pertinent to this part of the study has 
appeared previously as Eq. (II .35); its in- plane part is: 


^Geometrically this may be thought of as a 'Vandering ellipse", but, one which 
has a time dependent coefficient, as shown on the right. 


(11.49) 


(II. 50a) 


(II. 50b) 


(n. 50c)* 
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(II. 51) 


* 


I 2 B(«P) = [3(J 2 -J 1 ) + T 2 (2 V3 -)]A ai+ T 2 ^-){[l 2 -| • 

It should be noted that here the displacement vector may be separated into three 

component vectors, also; one for each of the coefficients: A , K and 7. Each 

a i °I 

of the components is described below so that the reader can become acquainted 
with how each adds to the overall geometry. For reference purposes the co- 
efficient vectors (A , K and 7) are repeated below. Thus: 
a i °i 

A +(J_ -2JJB R 1 + (J -2J )t! , 

a l 212 o '2 V 2 o ' 1 2'J’ 


and 


K = (2J - J J d - 2B& 1 , 
Oj ' X 2 o 2 o 


t . r c ). 


On close examination one finds that the A and K quantities are not 

a o 

precisely those noted above. This is because the initial state parameters have 
been transformed into a form compatible with the present reference triad. How 
ever, the force vector (7) (and its constant scalar components) are unchanged 
since these are not "rotated" to the new frame of reference. 


The first of the vector components is expressed in terms of A ( = A 

a, a 

(A^, A g ); with: 1 

A nlfs +2H 1 +T t 1, and A 0 =”-[h + -27 1. ) 

l2lo o £ J 2 2 L o o 77 J 


This particular partial solution is given by: 


*[' 


X “ [ S€T 2 - J X> +T 2<W*">] A a^* 


and, may be set down as the matrix equation: 
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A 2. 


(II. 52a) 


(II. 52b) 


(11.53) 
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The present displacement trace, on the (H, H) plane, is seen to be com- 
posed of: (1), a constant vector , 3{J -J )A ; plus (2), a circular locus of 

2 x a 

/ I 2 2 \ ^ 

radius ^ A j , The full figure is a circle of double orbit frequency {2tp) 

with its center located at (H, H) = (-3A , 3A ). Analytically, the figure can be 

C Xu 

described by: 

(\ S +3A 1 ) 2 + (6 1 H -3A 2 ) 2 = (A x 2 +A 2 2 ). 

This particular geometry is analogous to the ellipse found for the (£, 77)-plane, a 
trace which evolved as a consequence of the coefficient, A . 

The second contribution is acquired from the constant, K . The corres- 

°I 

ponding partial solution is: 


6 [y^>lr 


or, as a matrix equation: 
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6H * 


COS (ft 

-sin<p 

3m 

— X 

r . n i 

smip 0 

} 

r i 
K 1 

6H 

j 

„ i 

sin<p 

cos^ 

. 

2 

-cos (p 0 

i 

- K 2. 


The scalars K. (above) are defined from: 

K = K (K K ), 

Oj Oj 1 2 

where, in terms of initial values: 

K s 2 (H +H»), and K = - (H + 2H»). 
1 ' o o' 2 x o o' 

An analytical parametric equation for this trace is: 

(« 2 B,2 + <*2 H ) a - K l 2 + (? K l- K 2) 2 - 


Geometrically the figure is a spiral originating at a locus away from the 
coordinate origin. 


(11.54) 


(11.55) 


(II. 56) 
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On the (4 , 7?)-plane the trace due to K q was a line parallel to the 4 _ax ^ s - 
Here, one could show that the present trace is merely a transformation of the 
uniformly moving point (on the ( 4 , 7?)-plane) into the spiral (on the (H, H)-plane). 


Once more it is seen that a secular influence developes from one com- 
ponent of K^, namely K^. When this quantity vanishes the resulting geometry 
is a closed figure. 


The applied force system, T (= t(t 
trace solutions as follows: 


4’ 




contributes to the partial 


) 'V- 


This expression is recognized to be (simply) a transformation of the earlier so- 
lution (see Eq. ( n.46a)). In a matrix format the expression is written as: 


6l. 
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o> 

hi 


[an 

to 

II 


-sin<p 

cosp 
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which can be expanded to; 

6 H^p) -Ti. [cos<p + 2<psin<p ] 
3 £ 

and 

6 H(^>) [sin^3-2^cos^] 

3 £ 


Ci 

+ T [2pcos^>- (4 - ^~) sin<p] , 
T} 2 


+ 7 [2<psinp+ (4 - ) cos<p] . 

V 2 


Apparently the specific force (T) leads to a divergent (or secular) condi- 
tion which cannot be circumvented without removing the force, per se. 

A quadric expression for the trace is: 1 

2 

(5 3 E) 2 + (6 3 H) 2 = [T ? + 2 V 3T r) ] 2 + [2<pT^~ (4-^|~) T^] 2 . 


This expression describes the spiral-like influence typical to these 
figures. Actually each force component, alone, leads directly to a spiral 
(form) as its contribution to the overall geometry. 


(II. 57a) 


(II. 57b) 


(II.57c) 


(n.58) 
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The general traces produced here, as in-plane displacement geometries, 

are a composite of the figures generated from the partial solutions above. Thus 

the curve found on the (H, H) plane would be a form of spiral generally, but one 

incorporating the circle described through A . 

I 

These vector descriptions provide a convenient, mathematical means to 

represent this solution; however, this is not a scheme which is always consistent 

with the physics of a problem. For instance, to show the influence of initial 

values on this solution, one should include influences from both A and K .To 

a i °i 

clarify this situation, and that arising from the specific force (f), the next dis- 
cussions are included. 

The effect of initial state (ft , ft 1 ), alone, on this solution is obtained by 

o o 

setting r=0 in the analytical results. For this case a partial solution would 
appear as: 

A [l/OP ) \ = [3(J 2 - J 1 ) + T 2 (2‘P -)] v /T 2 «p-)[l 2 - | 

wherein [A 1 = A (A , A ), but with: 

L ajd a r 1 2 

A + 2H* 1 , A s “ [*H +H»1. 

12Lo oJ 2 2 L o oJ 

As before, it is assumed that 

K =K (K , K ), (see Eqs. (II. 55)). 

Oj Oj l 2 


In matrix form Eq. (II ,59a) is: 


[y^>]i s 


r AE 


Ah 


-3 +cos 2<p 
sin 2 <p 


-s in 2 (p 
3 +cos 2 <p 
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cos<p- 
sin (p - 


^ sin<p 
2 




cos <p 


-sin<P 
cos <p 



(II. 59a) 


(II. 59b) 
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This expression suggest a geometry composed of: (1), a circle, dis- 
placed from the origin, due to A ; plus (2), a spiral, from the K term. Here, 

a i °I 

again, the figure's divergence develops as a consequence of the scalar, K^. 
Therefore, if that term is removed the trace becomes a closed curve with a 
period matching that of the circular base orbit. 


The non-secular trace geometry, defined for = 0, is described by: 


I 1 / Hr 


AK(ip) 

= 1 


“ 2 
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-3 + cos 2 (p -sin2cp 
sin2cp 3 +cos 2 (0 J 


-j r 
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H +H» 
o o J 
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r 3(0 , „ . 

cos<p- -- sm(0 -sm <p 

3(0 

[ sin^-*~cos^ cos(0j 


0 

K 


Generally speaking, this equation describes a limacon, 

A special case of the limacon arises when HJ s 0. For this condition, 
added to K “ 0, a symmetric cardioid results; one which is symmetric about 
the H-axis. 

The in-plane displacements associated with the applied force (T), are 

described in terms of A and the transformed vector, T Up “) r. The matrix 

aj 2 T 

expression for this partial solution is: 

a[i 2 %»] 2 - 


AH((0) 

_ 1 
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-3 +cos 2(0 

-sin 2(0 

AH(<P)_ 

~ 2 

2 

sin 2(0 

3 +cos 2(0 j 


-2 T 
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3(0 

cos(p + 2(0sin(0 2(0cos(0- (4 - — sin(0 

3(0 r 


T i 


|_ sin<p-2(0cos<p 2(0sin<p+ (4 - — — ) cos(0 J I T I 

2 * TJ J 

Here the first matrix describes a displaced circle (of double orbit fre- 
— 2 2 

T ^ + 4T ). The second term represents a spiral-like 


(II.59c) 


( 11 . 60 ) 
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figure ; thus, the combination is a spiral-like trace typical of inertial, in-plane 
displacements. 

This is the last description for in-plane traces referred to the two 
reference triads. The next section describes some of the companion hodographs, 
for the situations noted above. Since these discussions are restricted to the in- 
plane cases, all considerations of out— of— plane traces are deferred until later, 

(c) In-Plane Hodograph; Rotating Frame. The hodographs, or velocity 
diagrams, for this relative motion problem are obtained from Eq. (11.29). The 
in-plane portion of that expression is: 

X 2 ^ , ^) = [l 2 -3(J 2 -J 1 )] B 2 T 2 <^)A a -f flWW, 

where the constant parameters A , K and T have been set down in Eqs. (II .43) 

ct O 

and (II. 45a). 

The geometric figures to be described here are obtained by a procedure 
analogous to that used for the displacements. That is, selected situations are 
described to show how the various parameters affect the motion and its traces. 

A first contribution to be examined here is that due to the quantity A . 
Recalling that this term is acquired from the trigonometric parts of the dis- 
placement solution, then making use of the same definitions as before one finds 
the partial solution: 



! 

05 

iTff 


sin<P 

cos <p 

V 

(II. 61b)‘ 

J1 

& V' . 
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2 eos<P 

-2 sinp 

- A 2- 



Here, also, 


A a A a (A l’ A 2 )- 
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It is seen that this expression describes a (2:1) ellipse centered at the 
coordinate origin. A cartesian equation for that figure is: 

(\ V ) 2 (\ v ') 2 

(This is an expected result also; its equivalent has been obtained and discussed 
in Ref. [l]). 


(11.62) 


The second component vector equation is, from Eq. (II .47): 


r - -i 3 

*64*' 


0 o ’ 


2(24 + T) 1 ) 

n 

L I 2 Jl ' (¥,) J2 S -i (B 2 J l> K o = 

.by* . 

2 

-I 0 

l 2 J 


V -H' 

'o 0 J 


(II. 63a) 


This plots as a single point on the (4 7} 1 ) hodograph plane. Obviously, 

this expression defines the motion rate associated with Eq. (II. 45a). 

Next, the influence of T on these hodographs is found from: 




64- 


0 2 
-2 -3<p 


T 


Here the trace is a line parallel to the t^'-axis; it represents a point 
moving, continually, in a specified (±) direction. An equation for this line is 
easily found as: 

6 n' 

In summary, these three influences produce a trace which is composed 
from an ellipse moving away from the origin (in an rj 1 -direct ion), and which 
has an offset central point. The offset is due to the combined influence of K 

o 

and T. 


(II. 64a) 


(II. 64b) 
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As noted earlier these component expressions, taken independently, do 
not define a physically real situation, hi order to clarify this and to illustrate a 
more realistic problem the equations are reformulated as shown below. 


The contribution of a prescribed initial state to the relative motion hodo- 
graph can be expressed as: 


a [v*’Hi s 




sin^> cos<p 

3£ +27?* 
0 0 

+ 

0 

Ar}' 

1 

2 cos<p -2 s in <p 

. - 


_-3(24 o ^;). 


Examining this result the figure is seen to be a (2:1) ellipse not 
centered at the coordinate origin. If the secular factor is removed (i.e. s 0) 
the hodograph trace reduces to the (2:1) ellipse, 


( A ^’) 2 ( A ^ 1 ) 2 

(Necessarily, the quantities describing this ellipse may be varied according to 
how is made to go to zero). 

Finally, to determine the influence of r alone, on the hodograph, 
consider the matrix equation: 


a[i/>)] 2 - 


> 
i 


P* 

2 


sin<p cos <p 
2 cos <p -2sin<p 



* - 


1 

CM 

o 

1 


T h 


4 

+ 



4 


-2T 





J 

V 


-2 -3<p 


T 


» d 




. . 


Here, the first matrix is obtained from A , while the second comes 

a 

directly from terms in (^^.T). The geometry associated with the first part of 
this expression is, again, a (2:1) ellipse while that defined from the second is 

i a line locus paralleling the 7 )’ -axis. Consequently, the composite geometry would 
appear as a "wandering ellipse", or cycloidal curve . An analytical expression 
for the trace is : 


(11.65) 


( 11 . 66 ) 


(II. 67a) 
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(II. 67b) 


| V A 2 n ' ) ~ r i (A 2^' ) = <4T V " T 4 2) sin ^>- 

Physically, one should recognize that each of the T scalars affects the 
hodograph differently. For instance, suppose that T^ s 0; then, 


and 


A z t'~ zr r 1 = - ar^cosip, 


A rt 7?'+3<pT = 4 t sin<p. 
2 ?? 7? 


These equations describe a regular cycloid on the hodograph plane. 


Next, if T = 0, the hodograph's trace is found to be the ellipse: 
V 


<V>‘ 


<A 0 ??'+2r t ) 




r i 


(2 r^y 


(Sketches of these curves are found on Fig. II. 6 below). Note that the : 
secular character of the hodograph vanishes when is zero. 

This completes a discussion on this hodograph plane's geometry. In 
the next paragraphs a parallel study, for the inertial triad's hodograph, will 
be presented. 

(d) In-Plane Hodograph; Inertial Frame. A general expression for the 
relative velocity referred to an inertial frame of reference has appeared earlier 
in Eq. (II. 40a). 


Since only the in-plane geometry is of interest here, then the applicable 
equation is: 

Js'ifpy, 


I 2 B'(«p) = 


H'fP) 


= B2{[2 T 2^-)] A a i +T 2^')[( I 2-|’ 5B 2 J l) K o I + ’ I 'r ? ]} 

+ t 2 <p")[*; ? -| (Vi)k 0 ]• 


(11.67c) 


(II.67d) 


( 11 . 68 ) 
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Sketch (a) 


Sketch (b). 


Graphs illustrating the effect of r on the relative motion hodographs« 
Sketch (a) describes an influence due to Tn\ while sketch (b) shows 
the consequence of T^, applied separately (see Eqs„ (n . 67))« Note : 
These curves originate from the coordinate origin (A 1 =0). 











It should be remembered that A and K are primarily written in the initial, 


1 


I 


inertial state of the motion (ft , ft’ ), though A does include the applied force. 

o o 

(See Eqs. (II. 34b), through (11.37) for descriptions of these and other quantities 
of interest in this expression). 

As a first part to this study consider those contributions made by the 

three parameters, A , K and T. 

a I °I 

The vector A influences the velocity trace in the following manner: 
a I 


6[i 2 R»]i " 


<5H 1 
6h* 


= 2B. 


[VP')] 


J X 


= 2 


-sin 2^5 -cos 2 <p 

V 

cos 2 <p -sin2<p 

• A 2. 


Here, A s ~[s +2H r +r t 1 , and A 0 s ~|"h + H* - 2r 1. 

12Lo o £ J 2 2 L o o 7 ) J 

This geometry is recognized as a circle centered at the coordinate 
origin, but having a motion frequency which is twice that of the base orbit. 
(The quadric equation for this figure is: 

(6 1 S') 2 + (6 i H') 2 =A 1 2 +A 2 2 ). 

By comparison, the corresponding trace on the (£*, r ) ? ) plane was an 
ellipse, traced over once during each orbit. 


The second partial solution for this hodograph is due to the constant, 
defim 

[^j ={[ B 2 T 2<’ 3 ‘>][ I 2-f’’ B 2 J l]-| T 2^">[Vl]}V 


K . With K defined as K (K , K ) then: 

Oj Oj o r 1 2 


(II. 69a) 


(II. 6 9b) 


or, 
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* x 

— (sinp + 3pcos<p) 
2 

-cos (p 


K 1 

_5H'. 

2 

1 

— (-cos <p + 3<p sin (p) 

" 2 

-sin <p 


K 

l 2 J 


— .1 i— t r =o/W xllt \ _ 

WllCiClXl, IV, III 1 - 1 'n / «*»*'“• 

1 O O 


K = 


- (H + 23* ), 
'“o ” o 


The secular terms in this expression produce a trace which is the involute 
of a circle; however, the trigonometric terms alone describe a circle. When the 
components of Eq. (II . 70a) are joined they lead to the following quadric equation: 

( 6 s-> +( B h») =( — kJ +k 2 . 


It is noted here that when vanishes (a condition removing the secular 
influence for the displacement geometry) the resulting figure is again a circle: 

(6 2 S ’ )2 + (6 2 H ’ )2 =K 2 2 • 


This figure is centered at the origin of coordinates and has a radius = K . 

Lt 

For the converse situation, that of having K vanish, the figure reduces 

Ld 

to a spiral originating away from the coordinate origin. 


The last quantity encountered in the relative velocity vector is the specific 
force (T). The particular solution describing its influence is: 


6H» 

V/’Ha" 6h , 


or, in an expanded format: 


6H* 


[sinp + 2(0cos <p] — 2) cos^+psinpj 


■ 

6H r 

J 

3 

[-cosp+2<psinp] \jr~ 2 — 2) sinp-pcos pj 




An examination of this result quickly shows divergence is present here. 
(The scalar equations here can be cast into the equivalent form: 


(II. 70a) 


(II . 70b) 


(II.70c) 


(II. 71a) 
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(II. 71b) 


(6gH ') 2 + (SgH') 2 = ( ) 2 - 6 (r^ r^)<f>+ [{2r^ ) 2 + 6 (t^ - 5 (r^) 2 ]<a 2 




This expression is composed of the force constants plus various ordered secular 
terms). 


Due to the explicit dependence of these equations on <p it is obvious that 
the resulting geometry is a spiral, also. 

Recalling that A a and K lead to a circle and a spiral, respectively, 
I °i 

on the hodograph plane, it is obvious that the general figure will be spiral-like 
too. 


The next task is that of determining the consequence of a prescribed initial 
state Ry, and applied specific force (T), on the trace geometry, when these 
are independently applied. The two situations are examined below. 


The consequence of initial values on the hodograph figure is described 

by specializing Eqs. (II .49). Grouping the state parameters as: G = (H +H r ) 

1 o o 

and G„ = (H + * T ), then for present definitions ; 

2 o o 



[“'1 


1 — ■■ 
> 
a 

K 

1 


sin 2(p 
-cos 2 <p 


cos 2 <p 
sin 2(0 


G x + H 
G_ 


» 

o 



+ 


r i 

— (sin<p + 3<pcos<p) 

-cos<p 

“■(-cos<p+3(psinc0) 

. ^ 

-sintp 



Obviously the first matrix represents the (now familiar) circle of double fre- 
quency, while the latter describes a spiral. Once more the resulting trace is 
largely influenced by the spirals divergent characteristics. 


(II. 72a) 


Since secular effects are suppressed when = 0 (the condition intro- 
duced for the displacement case), the corresponding hodograph geometry becomes 
the closed figure described as: 
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(A, H ') 2 + (A, H’) 2 = (H' ) 2 +2(G„ +s ' ) [G (1 - C0S<P) - H' simp ] . 

1 1 O ct O la O 


This parametric equation defines a limacon. 


In addition to the constraint noted above, when H ! =0 also, an asymmetry 

o 

for the limacon is removed. Finally if G s H* (only) the classic cardioid is found 

la O 

as a hodograph trace. This is noted from a specialization of Eq. (II . 72a); i. e. , 
when the conditions indicated above are included, the equations become: 

A^H*(<p) = -H^ cos 2<p+2H^cos<p, 

and 

A- H‘(<n) = - H 1 sin2(p + 2H» sinP . 

1 v o o 

(These expressions may be combined to yield the quadric : 

(A 1 H') 2 + (A 1 H') 2 = ( 2 Hy 2 [|-cos¥>]. ) 


Next, the influence of (T) alone on the hodograph* s trace is examined. 

Since components of T appear in A , the proper partial solution here is: 

a I 


a [ i 2 «»] 2 " 


AH* 

Ah 1 


-sin2p + sin<p + 2pcos<p 2cos2<p+ 


(¥-0 


cos<p+<psinp 


cos 2^-cos<p + 2<psinp 2 sin 29 +(^^ — 2^) sin<p-<pcos9 




Once again the explicit presence of <p t seen here, eliminates the possi- 
bility of acquiring a non-divergent trace. In fact the figure for this hodograph is 
akin to a spiral; however, this one originates at the coordinate origin; and, it 
does have a monotonically increasing radius. 


Summary . Trace geometries for the initial values problem and the zero-initial 
values problem have been examined briefly in the paragraphs above. There 


(II. 72b) 


(II.72c) 


(II, 72d) 


(11.73) 
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displacement and hodograph loci were described as they would appear on re- 
presentative planes In the local rotating and inertially aligned frames of 
reference. 

One limiting constraint imposed on these solutions was that of a constant 
disturbance vector (T). This vector has fixed components parallel to the coordi- 
nate axes of the rotating triad. In order to examine an analogous case, but one 
with fixed components parallel to the inertial triad*s directions, a second problem 
is formulated and solved next. In that solution one finds that only the disturbance 
itself contributes new information about the motion traces. Hence, an examina- 
tion of these trace geometries will be a much simpler task. Therefore, in the 
next section a formulation and solution of the problem will be carried out first; 
then the trace geometries will be described, and discussed briefly. 

Descriptions in the foregoing paragraphs have been brief, and the 
illustrations there were limited in number. Recognizing that interested readers 
would like to see typical graphs for each of these various contributions, a com- 
pendium of results has been prepared. This catalogue of traces, for the various 
in-plane cases, is attached below. What one will find there is a set of typical 
traces; each one obtained for positive valued parameters {components, etc. 
from the vectors A^, K, and T). In most cases the matrix equations (only) 
are noted; however, the corresponding scalar expressions are easily developed 
from these; or, they may be found by looking into the volume describing trace 
constructions* . 

Before leaving this section it should be mentioned that a similar collec- 
tion of trace data (a compendium) is found behind the next section also. There, 
data are presented for the solution involving the specific force (T^). Also, to 
complete the catalogue one will find another data collection for the out-of-plane 
traces. This is found attached to the other descriptive materials. 


♦Interim Report, "Construction of Relative Motion Traces”, by J. B. Eades, Jr., 
AMA Report No. 73-39. 



y o DATA SUMMARY 


A Compendium of Data for In-Plane Trace Geometries . The collected data, 
which follows, illustrates in -plane motion trace geometries typical of those 
found in both frames of reference. These curves describe both the initial- 
values problem and the zero-initial values case where the disturbance is due 
to T(=r(T., T r) , T^)). 

Trace geometries appear in the following order: 

(1) 0 Displacements for the Rotating Frame of Reference 

(2) . Displacements for the Inertially oriented Reference Frame 

(3) , Hodographs referred to the Rotating Frame 

(4) . A Non-Secular Case, in the Moving Frame Notation 

(5) . Hodographs referred to the Inertial Frame 

(6) . A Non-Secular Case, in the Inertially oriented Frame 

For these various representations, the individual mathematical com- 
ponents are shown first (Cases 1, 2, 3), Next, the initial -values problem and 
the zero- initial values case (due to t) are shown. These more realistic cases 
are numbered 4 and 5. The last situation, 6, illustrates the full trace geo- 
metry — that due to both the initial values and the disturbance, T. 

Following the hodograph sections the reader will find a case study illus- 
trating a non-secular condition for each of the reference frames. 


It should be remembered, throughout, that all of these geometries re- 
sult from an arbitrary selection of initial value parameters (£, 7], £) Q and 

V) 0 and disturbance components (T^, r , T^), In all instances these 
quantities were chosen as positive valued constants. 
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In-Plane Displacements, referred to the 
Rotating Frame of Reference 
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In-Plane Displacement Diagram; Local, Rotating Frame of Reference . 
The describing equation: 


; A( ¥ ,)=[l 2+ 3(J 2 -J 1 )]T 2 (<p-)A a + [l 2 -|^ 


Vi 


K + \E r T . 

o T 


Partial Solutions , 

1. The construction due to A (the 'Trigonometric 1 ' Solution Part). 

a 

The expression for this partial solution is: 

6 [h 1 * > ] - [ J 2 + 3 ( J 2 - J l> ] T 2 *> ' ') A a ‘ = 2 [ 2 J 2 * J ! ] T 2 (<P "> A a ; 
wherein 

A a S 2 [/Vo + V 2J 1> B 2 A o +(J l- 2J 2 )f ]’ 

and 

r s r (r^ r^)„ 

2. Construction due to the K^-term (the Constant Initial Values). 
This partial solution is defined by: 





1. This trace geometry is an Ellipse (2:1) 
centered at coordinate origin. Rate of 
motion over this trace is <p . 



wherein 


K s (4J +j )L -2B A 1 . 
o 1 2 o 2 o 


2. This trace is a straight line, parallel to 
the Tj- axis. 


*0 'describes the <p~ 0 locus; 

P is a general point on the trace. 







3. This geometry shows a parabolic trace 
on the (4, Tj ) plane of motion. 



4. This trace is a combination of 1. and 
2. above. 



In- Plane Displacement Diagram (continued) 


5. The Applied Force Solution (with Ji , Ji 1 , s 0) 

o o 

is obtained from: 

A [ I 2 A " «»]" [ I 2 +3<j 2 - j i ) ] 

wherein 

S (A aV 3 l[ J l- 2J 2?- 

% 3 ( J 1 +4J 2 )-[ 2B 2 + !’’ J 2 ]’ 5 * 

and 

f£f(r 5.V r ^ 

6. A Full Solution (i.e., a combination of 4. and 5. above), 
is described from: 

’ r 2 A «P> = [l 2 + 3(J 2 - J i>] A a + [l 2 - f Vl 0 ] V V * 


€ 



5. (The trace geometry here is a des 
cription from X. and 3. above). 

i 



6. This figure is obtained as a combina- 
tion of 1, 2, and 3, or 4 and 5, above. 


*0 is the initial point for the trace; P is a general position on each curve. 



In-Plane Displacements , referred to the 
Inerttatly Oriented Frame of Reference 
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In-Plane Displacement Diagram; Inertial Frame of Reference 
The describing equation: 

I 2 « m - = [3 (J 2 - ) +t 2 (2 <? -) ] A aj+ T 2 «p -) {[vfew* X 

+ ^ T r}. 

Partial Solutions . 

1. The construction due to A (the ’Trigonometric" part of 

a i 

the solution). 

The descriptive expression for this partial solution is: 
6 [i 2 ««p)]=[ 3 (J 2 -J i)+ t 2 (W-)]a 


wherein 


and 


Vi[ I 2V B 2 (J i- 2J 2 )i ; +(J r 2, 2 |f ]- 


TST(T r T i?' V 1 

2. Construction due to the K (the Constant Initial Values quantity). 

°I 


This partial solution is defined as: 

5[l2««p)>T 2 (^)[ I 2-|( B 2 J 1 )p]K 


wherein 


K ■ (2J -JJft -2B ft* . 
o 1 2 o 2 o 



1. The trace is a circle with its center off- 
set from the coordinate origin. Offset 
is due to fixed terms in A a ; the circle is 
a consequence of the transform, T. 



2. This trace is a spiral-like curve; it is a 

consequence of the secular form of the 

expression. A displaced Initial locus is 

due to IK . 

2 o 


*0 is the initial point for this trace; P is a general locus. 




3. A spiral-like trace primarily due to the 

secular nature of ^ . 

T 



4. This spiral trace is a combination of the 
types shown as 1. and 2. above. Origin 
offset relates to initial state. 




1 1 - 



5. This spiral trace is related to the geometries 
described as 1. and 3. above. Note the 
curve originates at the coordinate origin; 
this is due to the null initial state. 



P 


6. (The special case of motion beginning at 
the origin (Particle Ejection Problem) is 
not shown here. This figure is, also, the 
combined solutions for 1. , 2. , and 3. above). 



Hodographs, In-Plane, referred to the 
Rotating Frame of Reference 
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Partial Solutions 

1. Trace geometry due to A (the 'Trigonometric" Solution), 

a 

The equation describing this trace is: 

6 [y ,(<P >]= ['2 +3(J 2 - J 1>] B 2 T 2 <*"> A a’ 

wherein 

A a 5 J[Vo +( V 2 V B 2 i ; + ( V 2, 2 )f ]' 

and 

fsf (VVV' 

2, The Contribution for K q (the Constant Initial Values term). 
This partial solution is given by: 

6 [l 2 I •(*>]= -fVi K 0 . 

wherein 

k 0 s < 4 v j 2 u V 2 V;- 


*0 is an Initial point for the trace; P is a general locus. 


v 



1. This trace is a 2:1 ellipse, centered at 
the coordinate origin. The motion rate 
over the trace is (p. 


4' 



2. The solution produces a single point on 
the hodograph plane. (Note that the 
point lies on the ?? r ~axis). 




4 ’ 



3. This trace appears as a line, parallel 
to 77 f -axis; trace motion is as shown, 
typically. 



4. This trace is a 2rl ellipse, with center 
shifted due to K Q -term. Motion and 
rate, over hodograph are as shown. 
(This geometry is a combination of 1. 
and 2. above). 



Hodograph Construction (continued) 


5. The Applied Force Solution (A , X ^ 88 °) bas its trace 
defined by: 

1 ’ <*>] = +s < j 2 - J l>>2 [ A a] T + *T*’ 

wherein 

[ A a]r S !< J r 2J 2 )f - 

«rJ.«-(2B 2 +3J 2 <0), 

and 

T-T( T V r v , T t ). 

6, A Full Solution (i.e., the combination of 4. and 5. above), 
is defined from: 

I 2 A 1 (P) = B 2 {[i 2 +3(Jl -J 2 )]t 2 «P-) A a - | J x K o } 

+ T. 

T 


*0 is the initial point for the trace; P is a general position on the curve. 


V 



5. The trace is now a meandering (2:1) ellipse; 
the center is shifted and moving in the 
negative ry'-direction. 

(This is a combination of 1. and 3. above). 



6. The geometry is also a combination of 1, 2, 
and 3, above. 



An In -Plane Non-Secular Case for the 
Rotating. Frame of Reference 
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DISPLACEMENT 


Non-Secular Traces ; Local, Rotating Frame of Reference 

To remove the divergent nature of these in-plane displacements it 
is necessary for (2 £ q + 7? T s 0). 

The resultant defining equation is: 

I 2 i (^)-[l 2+ 3(J 2 - Jl )]T 2 (^-)A a (1) ♦k-hVlK". 

wherein 

A «.iki +J o B oA'l. 

a 2Llo 2 2 oJ 
and 

k o )s ^o- 2b 2 3 1 1 ; =, 2 [v 2 vJ' 

00 (These are reduced forms for A and K , needed to eliminate the 
4 * a o 

secular influence). 

The Hodograph corresponding to the non-secular requirement (s) 
noted above is obtained from the expression shown below: 

\h » = B 2 { k ■ + 8^ ■ - J 2 )> 2 ftO A a (1) - f K o a) } , 
wherein 

A K ^ are as defined above, 
a o 



This trace is a 2 :1 ellipse with a shifted 
center (shift due to yc Motion on 

the figure is as shown. 

(A special combination of the displacement 
constructions 1. and 2. (or 4.) above. 
Another geomet^ for this case arises when 



This trace is also a 2:1 ellipse but centered 
at the coordinate origin. 

(This is a special combination of the hodo- 
graph construction 1. and 2. (or 4.) above. 


*0 is the initial locus in this plane; P is a general point on the trace. 



Hodographs, In-Plane, referred to the 
Inerttally Oriented Frame of Reference 
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In- Plane Hodograph Diagrams; Inertial Frame of Reference 


The general hodograph equation is: 

1/ ' «P ) - 2 [b 2 T 2 <2v' -> ] + B 2 T 2 (<p -) { [l 2 - |<PB 2 J x Jk^ 

+ V} +T 2 « 9 '){*V f -|[ B a J i] K o I }- 

Partial Solutions . 

1 . Trace construction due to A a (’Trigometric" Solution) 
The defining equation for this solution is: 


wherein 


6[i 2 B’(«5)]=2[b 2 T 2 (2^-)]a ; 

l 

A a S K I 2 S o + B 2 (J X- 2 J 2 )& ; +(J 1- 2J /] 


and 


t = T r/’ V* 

2 . The construction for K n (the constant initial value partial solution), 

I 

This solution is described by: 

6 [i 2 b- <¥>)] = t 2 «o{b 2 [i 2 - I V (B 2 J x )] - | (B 2 J x )} k 


wherein 


K o " V 2 Vi * 


*0 is the initial locus point, while P is a general point on the trace. 


4 



1 , Trace for this partial solution describes a 
circle with center at the coordinate origin. 
The trace motion is at double orbit 
frequency. 


E 1 



2 . This hodograph trace produces a spiral- 
like curve. The divergence is linked to 
the secular coefficient (B J <P). 

Z _L 





3. This trace produces a spiral; the 
divergence here results from terms 

in both St' and ^ ' . 

T T 



P 


4. The trace appears as a form of the 
archimedian spiral coupled with 
circles. 



In-Plane Hodograph Construction (continued) 




5. This trace is produced as a spiral-like 
geometry. Secular terms appear in both 
and > 3> l . 

T T 



6. This geometry is composed from traces 
in 4. and 5. above; or, a combination of 
1. , 2 . , and 3. 



An In-Plane Non-Secular Case, for the 
Inertial Frame of Reference 
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DISPLACEMENT 


Non-Secular Traces ; Inertial Frame of Reference 

To remove the divergent nature of the in-plane displacement 
trace it is necessary for (E q + H^) s 0. The resultant defining equa- 
tion is : 

1 2 *#» s [ 3 <J 2 "V +t 2 ( 2 ^')] A a < i 1) +T 2 (<Pl {[ T 2 1 

wherein 


r» + b & 1 ”| 

, and K (1) ■ - T 

jj +2B J 1 

L 0 2 0 J 

°I L 

2 0 2 1 0 J 



(A General Construction; 
Special cases not shown). 


(These coefficients are reduced forms of the like expressions noted above). 

O 

The hodograph. corresponding to the non-secular requirement (s) 
indicated above, is defined as; 


S 2 [ B 2 T 2 <2<P ">] \ ( r } + T 2* > " ) { ^['2 * I 

- - (B J )} K * 

2 { 2 rJ o. 


( 1 ) 

>1 


wherein 


A ^ and K ^ are as defined above. 
a i °i 



(A General Construction; Special 
cases not shown). 

These curves are limacons, in general; 
the special cases are noted to be varied 
forms of the limacon (e.g. ,the cardioid). 


*On these traces O is the initial point while P is a general position on the curve. 



VI. PROBLEM SOLUTION 


An Inertially Aligned Force System. In the foregoing sections of this report 
the relative motion problem was solved for the case of a disturbing force 
with components parallel to the local rotating frame of reference. The studs'- 
which follows , below, will differ from the former one in that now the force is 
assumed to have components in the inertial frame of reference. As with the 
former case, each of these components is assumed to be fixed in value. This 
will allow for an analytical solution to be obtained, directly, and will form a 
basis of comparison with the previous solution. 

Return, for the moment, to the differential equation for this problem 
(Eq. (II. 12c)); there the relative motion's acceleration vector is expressed by: 

r = P 2 [( I 2 - I 3 ) (r + P)+ 3I 3 P (?• P)]-2PB 2 ?+S F, 

wherein the summation on F represents all the external forces which may be 
applied to the test particle (Q). Also, here, r is the relative position vector 

A _ 

(for Q with respect to P); and, P - e , is the unit position vector for the 

X 

base particle (P). Other terms are as defined for Eq. (II ,12c). 

The force system of interest now is represented jn the inertial frame, 

hence: 

V 5F Vv F x’ v V- 

where the F scalars are the projections of F onto the inertial triad of 
reference. Recognizing that Eq. (II .74) can be written (for solution) in terms 
of local, rotating coordinates, then it will be necessary to relate F^ to the 
same triad of reference. This is most easily done by means of the transforms 

described in Appendix A. There, Eq. (A. 3a) relates the rotating triad (e 1 .) to 

“* o ^ 

the inertial one [(e^) ] ; accordingly, an inertial vector is transformed back 

to the local frame by this transformation. As a consequence one finds that the 

components of F^ , which are "seen” in the rotating triad's frame, are: 


(II . 74) 


(11.75) 
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where T<0 + ) = T(p + ), since <p is the transfer angle here„ Therefore, Eq. 
(II. 76a) is expanded to: 



" 1 
cos (p sinp 0 


F 

X 

-sin<p cos ip 0 


F 

Y 

1 

rH 

o 

o 

1 


F 

L z J 


Probably a more meaningful statement for T(p + ) is that noted in Eq. (A. 7c), 
Appendix A ; that is : 

T (<P + ) =I 2 cos(0 -Bsin<p + (I -I ). 


(II. 76b) 


(II, 77)* 


A Solution for the Relative Motion. The solution which is found here follows from 
Eq. (II .74) after Eq. (n ,76a), or Eq. (II .76b), has been included. Making use 
of the simplifying interpretations to terms in the differential equation (see the 
statements following Eq. (n .15a)), it is apparent that a solution is to be obtained 
for the expression: 



r 



2 (p B r +T(<p 
z 



(11.78) 


It should be evident that this equation differs from that in the previous case only 
in the last term — the disturbance function F. Here, as noted above, the 
scalars describing F are constants; however, the disturbance term itself is 
a variable through T(p + ). 

Before the formal solution is attempted it should be noted that for this 

case the relative motion displacement (r) is referred to a triad moving uni- 
0 

formly (<p = constant) with the base particle (P). The constant matrices (J fc ,B) 
are employed to properly "locate” the associated vectors in the triad space. 


*It should be recalled that (1^ -I^) - this is one of the special diagonal 
matrices used here. 



The state vectors (r, r) are those of primary interest; however, the initial 
state (r Q> r ) is arbitrarily set so that the full solution is obtained as a com- 
bined "initial-values problem" and a solution for a forcing function (hereafter 
called the "zero-initial-value problem"). 

Again, applying the method of Laplace Transforms, Eq. (11.78) is cast 
into the following transformed format: 


2 - 
s r 


( S ) - t 3 [ S ?o + ? o] = ' *3 ] ■ - 2 <° B {*3 [ S f(S) - ? o]} 


L 2 *2 


V 


2 ° 2 
s -up s +<p 


]V- 


(II. 79a)* 


or, in a more concise algebraic form: 

[l 3 s 2 -< fi 2 (3J X - J 8 ) +2(3b 2 s] ?(s) f A(s) ?(s)}= [l 3 s + 2<^B 2 ] ? 0 + Ig ? 0 

J 


r„s b 2 <p 


M 

L 2 .2 


_ 2 *2 
s + s +(p 


- If !- 

s J I 


(II. 79b) 


Formally, the solution for r(s) is recognized to be: 


. f(s)=[A(s)]- 1 [l 3S+ 2<iB 2 ]; o + [A( S )]- :l [l 3 ? o ] + [A(s)]- 1 [ ^ B2<P 


3 ] f r 


2-2 2*2 s 

s -up s +<p 


(II.79c) 


wherein [a(s)J * is the inverse of the matrix A(s) defined in Eq, (II. 79b). 
(This inverse matrix is set down in Eqs. (II .17)). 


After manipulating and rearranging, and finding appropriate inverse 
transforms, a solution for r(t) is found as: 

r (t) = {l 3 cospt + [4J l + J 2 ](1 - cos<pt) - s[j 2 B 2 ] m - stnpt)} ? 

+ (lg sinpt -3J 2 (tpt - sin^t) -2B (1 -cospt) 

(equation continued on next page) 





*For reference purposes, one should recognize that T(<p*l is transformed as: 

I 2 s _ J 

=^{T(^ ± )}=s^[I 9 cos^^B q sin^)+Jo} = + 2 ■ ^ + — , wherein "s" is 

S S 


the transformed variable. 
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+ { J 3 0- ~ cos ^ t ) + 2 I 2 ^ tsin< ^ t ) “ f b 2 <^t-*tcos*t) 

.... a . p 

- JW +5J 2 j (1 -cospt- ~<pt sin^t) + -| J 2 B 2 (2<pt+<PtcosV3t -3 sin<pt)}-r| • 

Introducing the dimensionless variables shown in Eq. (II .21), and re- 
grouping the above expression in terms of its in- and out-of-plane components, 
a more useful form of this result is: 

i <P) = { [4Jj + J 2 >' 0 - *v ;} - H i 0 - 2B 2 i ; } COS *> + {6[b 2 jJi o 

- b 2 [ j 2 + 4 j Jv J sin ^- 3 Vi { 2 WV ;} fK +5j 2 ] 

-3B 2 J 1 ^- 8i n V) + |[j 1+ 2J 2 ][l 2+ PB 2 ][l 2 co S ^B 2 sin V ] 

+ — ["j + 4J Icos (pi r + J {ji costp+^T’ sin<p +r_ (1 - cos^)} 

2 L 1 2J J I 3 o o I J 


In this expression it must be remembered that <p plays the role of the independent 
variable (<p = (pt), and that all state and disturbance parameters are now in non- 
dimensional form. 


Eq. (11.81) is rather difficult to study and interpret due to its complicated 
form. Thus, it will be presented in a more compact format; the scheme followed 
here is directly analogous to that noted by Eqs. ( 11.24) through (II .27). 

First, the trigonometric terms having initial values as coefficients are 
consolidated, and grouped, using the matrices T(^>^). Secondly, the constants 
and the first ordered secular terms are separately grouped and named. 

(1). For the trigonometric terms, write: 

Vo cos<p + Vs sin<a s M + t 2 (<o+) M 


= r 2 DV^sJ COS(0+B 2 SlnV ’ 


(II. 80) 


(H.81) 


(II. 82) 
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where k and k are representations for the cosine and sine coefficients in 
c s 

Eq. (11.81), respectively. Recognizing that, now, 2 Q_ =k +k , and 2 Q, = 

a c s s 

k - k , then it is easy to show that 
C S -- 

tf a ^kv[v 2j i]Vo}- 

and 

^ = !f 3 VV[v 2J i]v;}* ** < n - 83a) * 

On examining these coefficients it is found that: 

^s =3(J 2" J l )<2 a ; (II,83b) 

here the matrix [j - J ] provides a proper change in signs, while the three (3) 

^ X 

is a needed multiplier. 

Incorporating Eqs. (II .83) into Eq. (II .82) it is easy to show that: 

I 2 k o cos <p + B 2 k g sin<p= {t 2 «0 + 3T 2 («5 + )[j 2 - ^ , 

after some manipulations. 

(2) . Next, when the first term (the constant parameter) in Eq. (II . 81) 

3 

is multiplied by “ J the result assumes a form which is very similar to the 
secular coefficient in the equation. Consequently, these two quantities can be 
combined as follows: 

{ k + j 2 >' o - 2b 2 4 " ;} - sb 2 j i [ 2i 2 

4 2 -k B 2 j i]{k +j 2 ]V 2 V;>}- < n - 85 > 

(3) . Finally, for conciseness and uniformity of notation, the coefficient 
of Tj is replaced by (the equivalent set of terms): 


*When this result is compared to Eqs. (II .26), it is seen that (Z? a , # g ) - (A a ,A g )^ v ; 
i.e., the previous definitions, but now only initial value parameters are contained. 

**When this result is related to the previous one it is noted that the two are identical 
if in A a the T term is removed (see Eq. (II .27)). Hence, here, the term is re- 
placed by [A ]. , signifying this constraint on A . 

^ Vi 
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(II. 86) 


[ 2J l + 5J 2][ T 2« 5+) ^ -hl+f [ J 1 +2J 2 ][ J 2 B 2 +B 2 T 2 ( ’ 5+) J + 4 


When due account is given to all of the results noted above then the 
solution, Eq. ( II . 81), can be replaced by: 

i« = [' 1 +3(, fV] T !4Jl., + [ l f!»Vl] K o t *r f , 


+ J { J cos<p+£ T sin(0 + T (1 - cos<p)f , 
3 ^ o o I J 


wherein 

[ A a]i.v.i[ 3J l 4 'o +(J 2- 2J l ) V;]> 

k ^K +j 2 ] j o- 2 v;' 


and $ is the quantity given by Eq, (II , 86). This solution should be compared 
to Eq. (II ,30a), the solution for A p) when T (the specific force, referred to 
the rotating triad) was the ’’disturbing function”. A comparison of these two 
mathematical statements will show that the only new information provided here 
is associated with the terms. Hence, the initial value problem has remained 
intact (as would be expected) while the ’’zero-initial -value problem” is altered 
accordingly. 


The Relative Velocity Equation. The equation in question, now, is most 
readily obtained by the differentiation of (say) Eq, (II , 87). One must be care- 
ful,noting that: (1), the T (p ± ) terms bear derivatives, as does the coefficient, 
4^. In fact it is only the trigonometric quantities which are involved in this 
operation. Now, after differentiating and clearing, the resulting (velocity) 
equation is: 


(II. 87a) 


( 11 . 88 ) 
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wherein 


•*'«») = {i 2 + 3(J 2 - J i>} [ b 2 T 2 HKl V. - 1 [ B 2 j i ] V *r f r 

+ J 3 {*^0 C0S< ^ + ^ I "*^ o ) sin( p} » 

% S - [ 2J 1 +5J 2 ] B 2 + I [ J 1 +2J 2 ][ J 2 B 2 + ^2 -^ B 2 ) B 2 


iC J l- 2J 2 ] B 2 [ T 2^"> +T 2^ + >] * 


Summary, The descriptive equations for 4{<p) and A ! {<p) above denote a time 
history of the state of relative motion for the particle (Q), relative to the base 
particle (P). Here ”P n is assumed to move on a circular orbit about a simple 
attracting center ( fi ), There are no mutually reacting forces between the 
particles; however, Q is presumed to be acted upon by an external force 
which has constant valued projections referred to the inertial frame of reference. 

The state equations are given in a dimensionless form; nevertheless the 
dimens ionalitj!' of the problem can be recaptured by noting and applying proper 
multipliers (see Eqs. (11.21)). Also, the expressions above define the relative 
motion in the local, rotating frame. Consequently, when the motion is to be 
described in an inertially oriented frame of reference it will be necessary to 
take a proper accounting of the transformation which can accomodate this re- 
suit, (This is the task to be undertaken in the next section; comments are de- 
ferred until that time), 

A study of the expressions thus far developed shows that the separating 
of these results into an initial-value problem and a zero-initial-value problem 
is a natural division. In fact, there is no coupling between these solution types; 
hence, any new information which is to be gained here must arise as a conse- 
quence of the applied force (f^) alone. This reduces the number of situations to 

*See Eqs, (A. 6), Appendix A, for a description of the derivative forms of T^^). 


(II „ 89a)* 


(11.90) 
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be discussed, now, and allows problem solutions to be added to one another 
for a more complete definition and accumulation of effects and simulations. By 
combining the results here (for fp with the previous cases it is possible to 
imply a variety of applied force situations, in addition to the initial values 
effects. Of course, the reliability of these answers cannot be expected to ex- 
tend too far, into physical space around the base particle, or into time, since 
the analytical solutions are a consequence of the linearization imposed on the 
governing differential equation. This is a constraint on the problem; neverthe- 
less it is obvious that these formulae do provide an excellent insight into and 
understanding of the several situations studied herein. In this regard the reader 
can become much better acquainted with these classes of relative motions than 
he would otherwise have through (say) discreet numerical examples. 


In order to provide a better view of these solution equations the primary 
results, Eqs. (II ,87a) and (II .89a) are written below in matrix format. Recall- 
ing that the initial-values solution, now, is the same as before then the principal 
added effects are to be found in terms of the T^. Nonetheless, to avoid confusion 
the full expressions are written here. 


(a). From Eq. (II .87a) the in-plane displacements are given by: 



4<<P) 


-cos<p 

sin <p 

(V) = 

. W> . 

= 

2 sin<p 

2 cos (p 


34 + 2?? 1 

o o 

4 f 

o 



1 0 

r> 

+ 






J 

- “ 2 



44 +27J r 
o o 


I V -24* 
*- o o 


* 

r 3 

- 2 ( 1 - cos (p) ~sinp 

C4 

+ m 

sinp -cos(0 

) 

- ■" 

»— i 


-6sin(2 -5(1- cos <p) 

2 

2(1 + cos (ft) (2 sin<P) 

I 

7h . 


(b). The out-of -plane coordinate is: 

do i (P) s £(<P) = £ coscp+C* s\n<p+T (l-cos<p). 
3 o o z 


(c). From Eq. (II. 89a) the in-plane velocity is described as: 


(II. 87b) 


(11.87c) 
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* 


1 2 Jt '(<p) 55 


n'(<p) 


r - - si n<p 
2 


sinp cos<p 
2 cos <p -2sinp 


34 + 277 t ' 
o o 

.X . 


0 0 

I » 


44 +271' 
*o 

V -24* 

'o o -» 


-3cos<p -2 sin <p 


+ 3 


^-<pCOS<P 

Li Li 

(l-<psin<p) (^cos^O) 



r 


► 

!|J 



T*. 


(d). The out-of-plane component is: 


(II. 89b) 


'(<P) = cos (p+ (T -4 )sin<p, 

o z o 


(II. 89c) 


In the next section the case just studied will be reexamined but, then, 
with the motion referred to an inertially oriented frame of reference. For that 
analysis the state variables in use will be a dimensionless set; these are to be ac- 
quired by the same scheme as was used in the foregoing case study. 


The Relative Motion in Inertial Coordinates . In this section a solution to the 
relative motion problem is developed with the state variables described in 
reference to an inertial frame of reference. The method used to acquire these 
results is identical to that employed earlier; that is, the solution expressed in 
a rotating frame is transformed to the inertial frame. It should be remembered 
that the initial values appearing in the various parameter coefficients are also 
subject to the same transformations. (This operation has been carried out 
previously (see Eqs. (11.32))). The basic operation in this transformation is 
carried out according to: 

R(t) =T(<p’) r(t), (11.91)* 

where 

R(t) = R(X, Y, Z). 

Here R(t) is the inertially described displacement vector while r(t) is the 
vector in a local coordinate representation. The operator, T is the trans- 
form matrix needed to obtain the desired results. 

* It should be noted that R and r, here, are not dimensionless; for symbolic 
representations the physical coordinates are used. Thus time is the independent 
variable, etc. 
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Once this transformation is completed, and the position vector is properly 
expressed, then it is a simple task to develop the relative velocity equation. 
Mathematically this is accomplished by differentiation; symbolically this follows 

from Eq. (11.91). above, as: 

- 1 * 

R(t) = T(<p-) [r +¥> Br] , 

having taken into account the derivative of Eq. <A.6d), Appendix A. These ex- 
pressions represent the transformed state equations in dimensional form; the 
corresponding equations, expressed in dimensionless variables, are: 

«(<P)=T(p-) I«£>), 

and 

Recall, in these f ormulae ,that B and Ji , and their derivatives, (B f , A f ) 

! are written in the previous dimensionless form; the derivative here ls with 
respect to <p (the independent variable for the dimensionless notation). 


The Displacement Equation. A direct application of the first of Eqs. (II . 93) 
leads to the following expressions: 


> I 2^)={ T 2 ( 2 ^‘> +3 [ J 2- J i]}[ A a i ]i.v. +T 2^" ) { 1 2- MvJK^^Vr 

wherein A . , K and <£ maybe deduced and/or are defined in Eqs. 

l ajJi.v. o T 

(II .34b), and (II .86), respectively; and, for the third component: 

j cosp + B f s in + T (1 -cosp)} . 

o 3*o o I J 


For quick reference the coefficient matrices are repeated below: 

| A 1 + fjo-2J, IBoR'I , 

L ajJi.v. 212 o 12 lJ 2 oJ 
K = [*2J -J IS - 2B S' , 

o L 1 2 J o 2 o 

and 

♦The nomenclature here is such that Jl (<p) =*£(£, T), £); -A '(P) 5 I f (£ * 7 )\ £ T ); 
and B(cp) = B(E, h, Z); B f (P) s B’(H *, H ? , Z'). 


(11.92) 


(II. 93)* 


(11.94) 
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(11.95)* 


*T " [ 2J 1 +5J 2 ][ T 2^ + ) - r 2 ] + f [ J ! +2J 2 ]h B 2 +B 2 T 2^ +) 

+ i[ J l- 2J 2 ][ T 2 fP '>- T 2« 9+) ]' 

(The reader should remember that the previous and K q matrices have been 
revised by including a proper set of initial values; the consequence of those mani- 
pulations is noted above as A and K . Also, for this analysis the T terms 

a I °I 

of A are absent; the new matrix contains only initial values (i.v. ) — thus the 
subscript). 

Equations (II . 94) are to be used for the development of in-plane and out- 
of -plane trace geometries. Subsequently there will be some descriptions of 
these presented herein. 


The Inertially Described Relative Velocity . Following from Eqs. (II . 94) one can 
obtain the velocity definition for the present problem situation. Once again the 
analytical expressions may be arranged so that there is a logical separation, 
within the equations, to describe both the Initial - Values Problem and the Zero - 
Initial - Values Problem. By the application of Eq. (II . 93), or by differentiation, 
it can be shown that one form of the dimensionless velocity equation, referred to 
the inertial frame , is: 


(a), the in-plane component: 

V* «° » = B 2 { 2 L T 2 ^ ] [ A a ( ] i. V. + 1 T 2 ^ 1 [ 2J 2 - j l - 3 B 2 V ] K 0] 

+T ^)[ $ rH tT 2^[ $ ;]v 

and, (b), the out-of -plane coordinate: 


r- 


J -j K r coso- ^ sin<P+T sin<p}„ 

3 ■ 3 < o o I J 


*See Eqs. (A. 7), Appendix A, for the general form of the equations (T (fp*). The 
subscript (-J^, attached to these, infers the two-dimensional form of T 
(e.g., T 2 *I 2 T). 


(11.96) 


101 



The coefficient matrices A , K and $ are noted in Eqs. (11.95). The 

a i °i T 

matrix, a derivative of is found to be: 

T T 

- [ 2J 1 + 5J 2. ] B 2 T 2 ^ + I [ J 1 + 2J 2 ] { J 2 B 2 + h - ' P B 2 K T 2 (<P+) I 

+ ^ [ J 1 - 2 J 2 ] B 2 C T 2 ^ _) < n ‘ 97a > 

Summary, Equations ( 11,95) and ( 11.96) represent an analytical solution to the 
relative motion problem for that case when the particle of interest (Q) is acted 
upon by the central-inverse-square gravitational attraction field of force, and by 
an externally applied force. This external force is presumed to have fixed com- 
ponents in an inertial frame of reference. Also the state variables here are de- 
fined in this inertial frame. 

The previously studied problem parallels the present case, except that 

there the force had its fixed components referred to the rotating frame of reference. 

When these two situations are compared, directly, one notes that the initial-value 

problem (i.e, the solution in terms of B and B* ) is found in both. Some re- 

o o 

flection on this will lead to the conclusion that this is a natural consequence, and a 
logical happening, in view of the methodology applied for both situations. The con- 
clusion then, is that here the only new information which is provided is directly 

traced to the effects found for terms in (Remember that t t -T t (T m , T tt , T )). 

I I I a H z 

Once again the ability to separate the initial-values problem from the forcing 
function solution (Zero-Initial -Values Problem) is found to be of advantage. An 
obvious conclusion reached is that now one has the ability to consider the influence 
of initial values disturbances in the rotating frame, and disturbances in the inertial 
frame simultaneously or separately. This enhances the investigator's chances to 
simulate a larger variety of conditions and situations in his analysis. The added 
ability to combine influences and to examine multi-case situations will lead to a 
much better understanding of the overall relative motions problem. Of course 
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the specifies of various cases may not be in evidence, but the intuitive grasp of 
the problem can be significantly increased. 

As an aid to understanding the makeup of these solution equations the 

general formulae, Eqs. (II .95) and (II . 96), are expanded and presented below 

in matrix format. There each coefficient matrix is expressed in terms of 

appropriate initial -values. (Recall that A and K do not involve the f 

a i °i 1 

components. Also, it should be remembered that the initial -values solution 

here is the same as that for the former case (see Eqs. (II .39) and (II .40); 

these are repeated for completeness)). 


(a). The in-plane displacements may be expressed as: 


I 2 r (<P) = 


’s ftp)' 


* 

cos 2 <p-3 

-sin 2<p 

>o + 2h ;>' 

_H(<p)_ 


sin 2 <p 

cos 2<p+ 3 

i (H o +3 v. 


sin<p cos <p 


i/ 

( — 
o 

1 — l 

t 


HI 

2 J 



*\ 

rr 

1 +< 


»/ 

IL 


4-2 co s<p 5 sin p 

13 


-2sin<P 


-5cos<p 



p 7.i 

-2cos2(p sm2<p 


■ 1 4 — 

-(sin <p+— sin 2<p) j (cos 2<p- 3) 


V 

+ 

7 

+ 3 <p 

1 , l . 

> 



-2sm2p — cos2(p 

L 4 J 


cos<p+7(cos2<p+3) — sin2<p 

l 4 4 _i 


. T H. 


(b). The out-of -plane component is: 


j = Z(<p)~Z cosP + Z r sinp+7* (l-cos<p) 
3 o o z 


(c). The in-plane velocity components are, from Eq. (11.96) 
S f (p) -sin2<p -cos2<p 

H'(<p) cos 2<p -sin 2<p 


* 

"S +2H t ’ 


o 0 


H + B> 

J 

L o o. 


(equation continued on next page) 


(II. 97b) 


(II.97c) 
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cos ip -sinp 


1/ 

i 

w 

o 

1 



1 


A 

2 

L 


2(2. + H' ) 
o o 


;S 1 

oj 


f, 

+ ^7 

13 sin 2<P 

“(9 + 11 cos 2 (p) 

+ 

-sin<p 

5 cos<p 

r 

(9-13 cos 2 (p) 

-(11 sin 2<p) m 


cos (p 

5 sin^O 


3 (p 

' cos (p+ — cos 2 <P 
1 

sin2<p ' 

] 

rr„ i 

M 


sin<p+“ sin 2<p 

-cos 2 <p 

1 

T 

L hJ 


(d)» The out-of-plane component is given as: 

J>((p) = Z»to)==Z* cos<p+(T -Z ) sin <p. 

3 o z o 

Following the pattern established earlier, a next task is that of describing 
geometries which represent motions on the various displacement and hodograph 
planes. In the next section this is done; there, graphs and mathematical definitions 
will be displayed. 


(II. 97d) 


(II. 97e) 
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VII. GRAPHICAL DESCRIPTIONS 


Geometric Representations, For the previously described traces, the relative j 
motion figures represented a particular family of general solutions. There the 
complementary function had added to it results acquired by allowing the specific 
(perturbative) force play the role of a forcing function. For that case, however, 
the applied force term was represented by a system of components which were 
aligned with the local horizon frame of reference. | 

In this, the present case, the only new information which is to be acquired 
arises solely from the forcing function — the applied force which is to be considered. 
Necessarily the new solution is altered in view of the fact that the "force” is repre- 
sented, now, by components parallel to the inertially aligned frame of reference. 

(It is recognized that the complementary solution is not altered; thus, the initial 
values problem retains its singular identity for both cases). 

Assuming that the forcing function is representative of the inertially aligned 
disturbing force, then the consequences of this, on the geometrically described 
traces for the motion, will be ascertained. Of course, one must recognize that 
this quantity will affect motion traces in both frames of reference; therefore, a 
full investigational expose will evolve and be discussed in the next few paragraphs. 
Here, then, will be descriptions of the displacement traces and the hodographs, 
for both frames of reference; but, only the consequence of the disturbing force 
will appear. Generally speaking, these descriptions will be presented primarily 
in an analytical format; nevertheless, a few figures are shown to indicate the 
trends which have been found. For a more complete graphical description the 
reader is referred to the compendia of results found at the end of this section. 

It should be mentioned that for the present only in-plane results are dis- 
cussed and described. The out -of -plane motions, traces, etc. are deferred until 
later; at that time all of the out-of-plane motions and geometries will be pre- 
sented, in one section. 



/ 

/ 
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(a). In-Plane Displacements; Rotating Triad . The general analytical 
equation for this problem is: 


V <*> = k +3 < J 2 - l.V. + t 2 - I” B 2 J l] K o + Vi • < n - 98 > 

However, here, only the partial solution dealing with Tj is of consequence; 
therefore interest is focused on: 


a [vH s Vi ; 


(II. 99a) 


wherein 


*r £ K + 5J 2 J [ T 2 <^ + ) - >2 ] + I 2 [ J 1 + 2J 2 ] [ J 2 B 2 + B 2 T 2 <*> + >] 


+ 4[ J l- 2J 2][ T 2« , ‘ ) - T 2 <P+) ] ! 

and 

T = T (T w , T , T ), with each scalar being a constant. 

I I *e> H Z 

Examining Eq. (II ,99a) one finds that, geometrically, the trace is 
conveniently constructed as a four vector sum. One vector for loci on an 
ellipse; two describing line positions; and, one, an archimedian spiral. When 
the full geometry is displayed it is found that the figure looks as "spiral -like"; 
hence, it will be classed as such. A description of this system of vectors is 
developed below: 


(1). The first partial solution (vector) will be described by: 

6 l[ I 2 i( ' 9) ]"[ 2, l +5J 2 ][ T 2« ,+) - I 2] f I- 


(II. 100a) 


Expanded into a matrix format it becomes : 




*2(l-cos<p) 

-2 sin<p 


* ’ 

M 


1 

5 sin9 

5(1 -cos<p) 
. 


. t h. 


(II. 100b) 
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It is an easy exercise to show that this represents a (2:5) ellipse with an origin 
at (4> V) c £ {-2T 3 ,-5T H ). 

An analytic equation for this figure is expressed by the quadric equation: 


/ V + 2T H \2 / \^T 


)*( 


H 




■) 2 = x. 


Necessarily this in-plane trace originates at the coordinate origin; however, 
the figure size is explicitly dependent on the force magnitude (|f | ); and, it 
has a periodicity matching that of the reference orbit. 

(2). One of the ''line" loci found here comes from: 

6 2 k^ £ ?[v 2j 2 I j 2^]v 

As a matrix equation the expression is expanded to: 


6 4 (<p ) 


"o 0 * 


- - 

M 

Otjffi) 

ft 

3(f> 0 


- r H . 


This is a line parallel to the 77-axis; its direction is controlled by sgn (T ); 

H 

and, the rate of displacement along the line is 1 3T I . 

H 

(3). A second "line" trace is developed from the partial solution: 

‘sfc* «“>]• |[ J 1 - 2 J 2 ][ T 2 (V3_) - T 2 (<P+ 0 ; 


or: 


64(<P) 


6 W) 

L J 

3 


sincp 

2 


sinp 0 




L H 


These points (4 , ??) describe a line passing through an origin once each half 
orbit. The amplitude of oscillation here is: 


i6 3 Lax" • 


(II. 100c) 


(II. 101a) 


(II. 101b) 


(II. 102 a) 


(II. 102b) 
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(II -103a) 


(4). The remaining vector from Eq. (II .99a) is: 

6 Ji l 2 £ H" ¥ [ J 1 +2J 2][ B 2 T 2 (<P+) ] ; I • 
As a matrix, this equation can be recast as: 


■8 {#>)' 

= 3£ 

sin<p 

-cos<p 


' 

1 JA 

_8r?«0)_ 

2 

4 

2 cos <p 

2 sin^o 


. T H. 


Geometrically the equation is represented by r = k<p , when a new set of 
coordinates (4, 77/2) is used. That is: 



It is seen that, here, the spirals parameter (k) is |- |r (in-plane) | ; the figure 
completes a circuit once each orbit; and, this solution describes a divergence 
for the full trace geometry. It is evident from Eq. (II .103b) that the secular in- 
fluence of r cannot be removed, without eliminating itself. 

When these partial solutions are added the complete in-plane trace, due 
to Tj» is obtained. Obviously, as the particle moves under this (f action the 
secular influence plays the more dominant role. Therefore, in time, the over- 
all trace becomes T, spiral-like M , and ultimately it will exhibit this characteristic* 
alone. (See Fig. II. 7 for an illustrative sketch). 

(b). In- Plane Displacements; Referred to an Inertial Triad . When the 
above problem is referred to an inertially oriented frame of reference the in- 
plane trace geometry is found to be significantly altered. In order to examine 
this in some detail a solution corresponding to (see Eq. (II . 94)) must be de- 
veloped. Thus, the partial solution to be studied next is; 

A[i 2 ft<*)> T 2 «n[Vi]- 

*A more detailed discussion of these figures can be found in Reference [2] . 
There it is shown that the individual components of Tj leads to spiral-like 
figures also. 


(II. 103b) 


(11.103c) 


(II. 104a) 
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When the coefficient function, T ftp") O ] , is examined it is found that a con- 

Lj f 

venient description of this trace geometry is acquired as a sum of six vectors*. 
These vectors are developed from an expanded form of Eq. (II .104a). The com- 
binations. which are suggested here, are as follows: 


(a). First, taking into account the form of $ , one part of the solution is: 

6 [yM (K +5J 2 ][ T 2 (<P+) - r 2] + Svslw-v*]} V < n - 104b ) 


This equation leads to a three vector sum. Deleting details for expanding to ob- 
tain these, a convenient grouping of terms leads directly to the following three 
vectors : 


6.[l 2 ft(v3)_ 


62 (<p) 



i=l ,2,3 

-2 cos (p 

5 sin tp 

-2 sin tp 

-5 cos 9. 


4 0 
? 


7 • n 

-2 cos 2 <p - — sin 2 (p 

I -2 s in 2<p +7 cos 2 <p 

*■ 4 


^ r 


T_ 


L H 


(11.104c) 


From this expression one can see that 6 [i ft((p)] describes a fixed 

X ^ 

position in the plane of motion. The second component solution represents a 
2 2 2 7 2 

circle; ( 6 H) +(<5 0 H)^ = (2T„) + (— T ) ; and, the last partial result suggests 
2 .2 m 4 H 

a second circle. Here, then, one has a circle whose center is shifted from the 
origin, plus a second circular locus to complete the trace. (This grouping, now, 
actually represents a two-vector sum ). 


Yet to be considered is the remainder of Eq. (II. 104a). That part of 
the solution has the secular coefficient; it is: 


6 [ I 2 B#P) ] E ?[ T 2 (¥5 '>]{[ J 1 +2J 2 ][ J 2 B 2 +B 2 T 2 ( ’ 5+) ]} ? I ‘ 


(II.104d) 


Written in matrix format, this expression expands into: 


*As seen following Eq. (II. 104c) the first three vectors condense to a two - 
vector system; consequently the full solution is described by a five -vector sum. 
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<II.104e) 



6S(ip) 

-?{ 

"0 -3 " 

+ 

-4sin<p 0 


1 

i=4, 5, 6 ^ 

3 0 


4cos<p 0 



“ -sin 2<p 

cos 2 <p ‘ 


"T w " 

+ 



► 

hH 


cos 2<p 

sin 2 <p 

4 


.V 


Examining this result one finds that the first (5^) -vector (which makes 
use of the constant matrix) describes a line-locus, one which moves monotonically 
from an origin. The second and third parts of the solution represent spirals; i.e. : 


IV 


(II.104f) 


and 



+ (6gH) 2 



<p- 


(II . 104g) 


Unfortunately the spiral |§ | is not defined in the coordinate space, per se, 
but the modified space shown by Eq. (II.104f). 


The last trace geometry (6 ) is an archimedian spiral, with a parameter 

o 

which is proportional to the magnitude of (in-plane). 


The figure for the full solution here is found to be spiral-like, in the large; 
however, it does exhibit some unusual characteristics near the origin of motion. 

(A more detailed study of these partial solutions will identify which component 
(vector) is responsible for this)*. As an example of the trace produced here 
Fig. II. 8 shows a sketch for this case, but with the T^scalars chosen as positive 
numbers. Incidentally, it can be shown that if = 0 the resulting trace is a 
spiral -like curve without the 'Viggles" found otherwise present. 


(c). In-Plane Hodograph; Referred to a Rotating Triad . An analytical 
solution for the in-plane relative velocity, referred to a rotating frame of 
reference, has been obtained as: 


y^>= J l ) ] B 2 T 2^'>[ A a] l .v.-f [VJVVI • 


*See Reference [2] for more details on the construction of these traces. 


(11.105) 
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FIG. II. 8. An in -plane, relative motion figure in an inert ially oriented frame of reference. These dis 
placements result from the application of a specific force, ?j. See Eqs. (II .104) for a 
description of this case. 




However, the partial solution pertaining to r is all that is of interest here; 
thus, this investigation will examine only the result: 

wherein, 


T 


= -[2Ji + 5J 2 ]b T 2 (^) + |[j 1+ 2J 2 ][, 


J B + (I -B <P) B T\ 
2 2 ' 2 2 ; 2 2 


«p + >] 


+ i [ J a - 2 J 2 ] B 2 r T 2 -> + T 2 «P + > ] * 

An analysis of this expression finds that the trace, on the hodograph plane, 
is composed of a four vector sum. The individual vectors have been selected so 
that each can be described by a simple geometric figure. The addition of these 
vectors defines the full trace on the (£', t? 1 )-; plane. The selected vectors are 
described below: 

i 

(1). A first vector, extracted from Eq. (II .106), locates a fixed point 
on the hodograph plane; i.e. , 

6 1 [ I 2 A ''^] S f[ J l +2J 2 ][ ! J 2 VlV 


or, as a matrix expression: 




‘o o’ 


*T W ‘ 

6 tj'fto) 

h. * 

1 

3 0 

L J 


. T H, 


It is seen that this fixed locus lies on the 7? '-axis at a distance 3 |r w | from the 
origin. 


(2). A second function locus, obtained here, represents an oscillating 
line point of bounded magnitude which passes through the origin. The partial 
resultant is: 

fl 2 [v ,(v) ] s i[ J r 2J 2>2 [ t 2^-) +t 2^ + )] \ ■ 


(II, 106) 


(II. 107a) 


(II. 107b) 
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Expanded into matrix form, it can be shown the equivalent equation is: 




0 

-COS (p 


cos <p 
2 

0 



• * 


r w 


r 

J 

H 


(3). Next in complexity is the expression: 

^ 'H" f[ 2J l +5J 2] + |[' r i +2J 2]}[ B 2 T 2^ +) ]’ : i ■' 


or, in a matrix form: 
6£'(<p) 
6jj'((p) 


J3 


sin(p 

2 


costp 

2 


2cos(p 2sin^0 


-i 

- 




T 

J 

H 


The point loci, defined here, describe a (4:1) ellipse whose principal axes 

parallel the hodograph axes. Also, it should be noted that if the two vectors, 

6 and 5 , are combined the resulting geometry defines the same ellipse, but 
1 3 

with its center shifted to the point located by 3^. 

(4). The remaining partial solution, from Eq. (11.106), is: 


which can be rewritten as: 

r 


64 » 
6r?'«o) 


- 3<p 


cos (p sinff 
2 2 

-sin (p cos<p 


H 


A quadric equation developed from this is found to be: 

y[2(6 4 4’)] 2 + [6r,-]^ |6 4 | = 3 ; 

it describes an archimedian spiral in the (2£ 77 , )-space. 
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(rr.i07c) 


(n.l07d) 


(II.107e) 





If the four vectors (above) are summed they will describe the overall 
hodograph trace on the (£ r) 1 )- plane. Basically this curve is a spiral type 

(as shown on Fig. II . 9, below). It is seen that the secular influence is apparent 
only in the 6^-vector; however, this divergent character cannot be removed 
without eliminating the in-plane force, per se. Such a condition should not be 
unexpected in view of the fact that the force acts continually, in "fixed" directions; 
and, the divergence is obvious. 


(d). In-Plane Hodograph; Inertial Frame of Reference . A description 
of the (S', H')-hodograph trace is obtained from the analytical expression: 

#[l 2 B-(P)]-T 2 <p-)[ B 2* T + ^]ri, 

wherein and are to be found in expressions following Eqs. (II. 99a) and 
(11.106), respectively. 

After studying the full equation above it has been ascertained that this 
hodograph trace is most conveniently defined by a four vector sum, also. The 
vectors selected for these representations are set down and discussed below: 


(1). The first vector used here defines loci on a circle; one whose 
center is removed from the coordinate origin. The circle is found in the partial 
solution: 


‘6H'(P)‘ 

_ 1 

13 sin 2<p 

-(9 + 11 cos 2 (p)' 


"r w ' 

6H'(p) 

. _ 

4 

1 

(9 -13 cos 2<p) 

-11 sin 29 

-J 


.V 


9 9 

(Note that the circle's center is shifted to a position (H*, H') c - (- “ t , — T w ); 


its radius is f ^ T ) 2 ). 


(2). A second circle, of single orbit frequency, is also acquired from 
the general expanded form of Eq. (II .108). This one can be shown to be: 


(11.108) 


(11.109) 
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ras'W 

6h’(<p) 


*12 


-sin<p Scosp 
cos<p 5sin<p 


-*1 

" t w ' 


r 


L hJ 


7 2 2 

T w +(5T ) . 
tt H 

The secular aspect of the more general equation will produce two spiral 
geometries. The simplest of these is: 


'6E'(<p)' 


"cos<p 0‘ 


‘T w ’ 


= -3<p 





3 

sin^> 0 


, t h. 


which can be reshaped into the archimedian spiral equation: 


1 63 I = 3 T h <P. 

Finally, the remaining terms of the solution equation can be set down as 



S(p 

*cos2<P. sin2<p 

fiH'fto) 

~ 2 

4 

sin2<p -cos2<p 



"T w " 


r 

J 

L HJ 


It should be pointed out that this spiral has double orbit frequency; also, the 
spiral's parameter (k) is determined to be: 

k 1 7*1 (in-plane) | . 

Summing the four vectors (above) it can be shown that the hodograph's 

geometry is 1 ’spiral -like”, generally, but has a frequency of traverse which 

is double orbit frequency. (See Fig. II. 10, below. Since the vector summing 

introduces considerable geometric complexity, the interested readers should 

consult Reference [2] for more details on the construction). An observation 

worth noting is that the trace for t w = 0 is of single frequency and has the 

appearance of a regular spiral. However, the converse situation (T = 0) 

H 

produces a spiraling trace also, but one of double orbit frequency. 


(II. 110) 


(II. Ill) 


( 11 . 112 ) 
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Summary , The material presented in this section has been given to the des- 
cription of how the specific force 7^ affects the in-plane state of a relative 
motion. Primarily, the information set down above has been descriptive, yet 
the graphical representation of these results is limited. In order to expand 
this aspect of the investigation, a compendium of results (equations and figures) 
is presented below. There the reader will find typical traces constructed 
using positive valued coefficients. This collection of data should provide a more 
complete and comprehensive understanding of this aspect of the overall problem. 

In the section to follow, some comments on out-of-plane traces will be 
made and some sample situations illustrated and discussed. Following that 
descriptive material another collection of data will be found; these will pro- 
vide a compendium for that phase of the motion study. 
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VIII. DATA SUMMARY 


A Compendium of Data for In-Plane Trace Geometries. The following illustra- 
tions and information describe the influence of a disturbance 7 (=7 (7 tt , 7. 7 )) 

i i •=• Hz 

on the in-plane trace geometries. These graphs describe typical situations only; 
they are the consequence of a set of arbitrarily chosen positive valued parameters. 
Shown herein are: (1) the general case; and (2) special situations wherein 
each Tj component is applied separately. 

These data appear in the following order: 


(1) . In- Plane Displacements, referred to the Rotating Frame of Reference 

(2) . In-Plane Displacements, referred to the Inertially Oriented Frame 

(3) . In- Plane Hodographs, referred to the Rotating Frame 

(4) . In-Plane Hodographs, referred to the Inertial Frame. 
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In-Plane Displacements 


(a) Referred to the Rotating Frame of Reference 

(b) Referred to the Inertially Oriented Frame of Reference, 
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In-Plane Displacement Diagram; Local, Rotating Frame of Reference . 


The describing equation: 

I 2 1 <p> = [i 2 +3 <j 2 - Vi + [v - 1 b 2 jJ V Vi . 


Partial Solutions 


A L and K are identical to those 
L aJi.v. o 


1. The constructions due to 

obtained and described earlier. These are not reported here. 


2. The Applied Force Solution (with A =*0)* is obtained from: 


a[i 2 *"(<»)> <*YV 


wherein 


i—* 

w 

to 


V [ J l +2J 2][ J 2 B 2 +B 2 T 2 (<0+> ] 

4[V 2 X[ T 2 (,P ">- T 2<V]’ 


and 


T i s W V V- 


(a). Special Cases . 

The situations implied here result in traces 

due to T~ and to r alone. Each component is applied 
~ H 

separately, producing a trace on the displacement plane. 
(See sketch 2a). 


*The motion traces originate at the coordinate origin. 


4 



2. A typical in-plane trace, due to T , 


4 



2a. Typical traces, in the (£, 7?) plane, 

due to T w and T ' respectively. 

e H 
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In-Plane Displacement Diagram; Inertial Frame of Reference. 


The describing equation: 

l/^) = [ T 2 ^') + 3 < J 2 - J 1 )][ A a l ]i.v. +T 2 W '){[ I 2 -| ¥, B 2 J 1 ] K o i + 
Partial Solutions 

1. Constructions due to | A_ . and K n have been 

L a jJi,v. 

obtained previously. These are unaltered, in this problem, and 
do not appear here. 



2. The Applied Force Solution (with & o = &^-0)*. 
This partial solution is described from: 

a[i 2 %)> T 2<<0[Vi]* 

wherein: 

s >? [ J 1 +2J 2][ J 2 B 2 
+ B 2 T 2« 9 ^] + |[ J 1 “ 2J 2 ][ T 2«» _ )' T a<^]- 

and 

(a) Special Cases . 

The cases examined here lead to traces produced 
by the components T_ and r TT , considered separately. 

W I- XI 

(See sketch 2a). 

*The traces originate at the coordinate origin. 


2. Displacement diagram for and r H applied 
simultaneously; this is the combined trace 
from below. 



\ 

I 

/ 




H 


4 


O 


9 

\ 


H 




2a. Examples of the 
traces, due to 

and t h. applied 
separately. 



In-Plane Hodographs 

(a) Referred to the Local, Rotating Frame of Reference 

(b) Referred to the Inertially Oriented Frame of Reference. 
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In-Plane Hodograph Diagram; Local, Rotating Frame of Reference , 

This general hodograph equation is: 

-HvJv^v 


V. l.V. 


Partial Solutions 


1. Hodographs due to A a and K have been described 

i. v. 0 

previously. The contribution due to is the only new information 


to be obtained here. 


2. The Applied Force Solution ( A ~ - 0)* has its trace 


described from: 


wherein 


L 

n S - [ 2j l + 5J 2 ] B 2 T 2 (<P+ > + I [ J 1 + 2J 2 ] { J 2 B 2 

+ [*2 ' + 4 [ J 1 - 2J 2 ] [ B 2 T 2^ _) +B 2 T 2 (<P+) ]- 


VVV V V- 


(a). Special Cases . 

These situations lead to traces due to and r 

~ H 

as applied separate forces „ Typical trace sketches are shown 
on the right. 



*The traces originate at the coordinate origin. 


2a. Hodographs for and T 
applied separately. 



In-Plane Hodograph Diagrams; Inertial Frame of Reference. 


The general hodograph equation is: 


,K'^) = B 2 {[2T 2 (2^[A ai ] LVi + |T 2 «p-)[2J 2 -J 1 -3^B 2 jjK 0i+ T 2 ( V ,-)[$ r ? I ]} + T 2 (^)[^f I ] 


Partial Solutions 


1. Hodographs due to 


v. 


and K q have been des- 


cribed previously. These are not altered in this solution; they 
are not repeated here. 

2. The Applied Force Solution (for ft =ft^ = 0)* ( has a trace 
described from: 


wherein 


A[i 2 ftv<p>]^ 2 «n[B 2 v $ ;]v , . 

S [ 2 J 1 + 5J 2 ] h {tP+) - I 2 ] + f [ J 1 + 2J 2 ] [ J 2 B 2 + B 2 T 2 (<P+) ] 

+ 4[ J l- 2J 2][V>- T 2 <P+) ]' 

*; 5 - [ 2J 1 +5J 2 >2 T 2 (<P+) + 1 [ J 1 +2J 2 ][ J 2 B 2 

+ < V - V > b 2 t 2 (p ■ +) ] + 1 [j x - 2 J 2 ] [b 2 t 2 («p-)+b 2 T 2 «p + ) ] , 


and r I s r i (r H , t r) r z ). 

(a) Special Cases . 

These cases lead to traces produced by the separate 

actions of 7L and T. Sketches are shown at the right. 

«— H 



— t 

A 


2. Hodograph des- 
cribed for Tj 
applied to the 
test particle. 


trr 


2a. Separate effects 
of Tw and 
on hodograph. 


\ T~ 


*The traces originate at the coordinate origin. 



IX. OUT -OF -PLANE TRACE DESCRIPTIONS 


Qut-of- Plane Motion Traces. So far, in these discussions, the relative motion 
have been restricted to the in-plane cases. Now, in this section, the corres- 
ponding out-of -plane representations will be examined and illustrated. Once 
again the motions, as seen in both frames of reference and for the various 
problem types previously noted, will be discussed. 

Earlier in this investigation it was noted that the in-plane coordinates 
were coupled, hence there was a convenient and compact method of notation 
available to describe those relative motions. For the situations to be examined 
here, however, the coordinates and results are not coupled; this leads to an 
inconsistency in notation, one which adds a degree of complexity to the present 
case studies. In this regard it is not so conveient, now, to describe the geo- 
metric traces as simple vector sums; nor is it convenient to represent these in 
terms of elementary curve forms. Thus, only a selected number of cases will 
be detailed below; for a more comprehensive discussion of these geometries the 
interested reader should consult Reference [2] . 

After a review of the various situations described earlier it seems that 
a best grouping, for present purposes, would be according to physical conditions. 
That is, the initial-values -problems will be considered separate from the force 
disturbance cases. In addition, within each of these categories the graphical 
representations in the two reference frames will be described. Finally, in- 
cluded in with these will be discussions of the displacement and hodograph traces 
which develop for each. 

(la). Initial -Values Problem Displacements; Rotating Frame of 
Reference . The analytical solution for this problem type can be acquired 
directly from the general results, Eqs. (II .22) or (n .28). Since here, 

T = 0, necessarily, it follows that the out-of-plane traces retain those gen- 
eral characteristics which are typical of the in-plane cases. That is, these 
traces will exhibit the secular character found earlier; however, this effect 
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can be deleted by imposing the same restrictions as before. In this regard, if 

K = 0 (i.e. ; 2J JL -B I' s 0), by either of the obvious means, then the geo- 
1 ' 1 o 2 o 

metric divergence is removed and closed figures result. Then, if an added 

constraint is imposed, one whereby motion originates at the origin (A = 0), 

the resulting out-of-plane curves are found to be line traces* and ellipses*. 

On the other hand, if the non -secularly restriction is relaxed, here, but the 

condition of X =0 is retained then the line trace becomes an s -shaped on the 
o 

(Vi £) -plane and the divergence reappears. 


The several traces described above are not shown below in graphic 
form. Rather, the interested reader should examine the compendium of re- 
sults, located at the end of this section, or he may consult either Reference 
[l] and/or [2] . These documents are more explicit regarding the geometric 
construction of figures which evolve from the various situations noted above. 


(lb). Initial-Value Problem Displacements; Inertial Frame of Reference . 

This study also exhibits a secular characteristic for the out-of -plane traces. 

As was mentioned in the earlier discussions the geometric divergence can be 

removed if a particular combination of initial values is introduced. In terms 

of the solution format which has been developed, this secular influence vanishes 

for (K ) ^ 0 (i.e. , when J & - J B ft* = 0). To see the consequence of this a 
II 1 o 1 2 o 

bit more clearly, consider the situation of a relative motion originating at the 
origin (6 = 0). For this situation the geometry’s divergence is removed when 
H* o = 0* Physically, this implies that the test particle (Q) is set into motion by 
an impulsive action applied normal to the motion of the main (or reference) 
particle. The subsequent displacement diagram traced out by the test particle 
will be a closed curve, one which returns to the origin periodically. 

Analytical expressions describing these various out-of -plane traces, 
for those conditions indicated above, are acquired from the general solution, 
Eqs. (II .35); this particular non-secular case, for H r = 0, is expressed by: 

*From a broad point of view this motion may be considered as being represented 
by two ellipse-like traces; one a true ellipse, the other (the line) being a limit 
case for the ellipse. 
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<n. 113) 


S (<p) = 2H^ (sin<p- “ sin 2 <p), 

H o = 2H o (\ cos2( P- cos<p ) » 
and Z(<p) = Z^ sin<p. 

(Note that the speed components which appear here satisfy those non-secular 
initial conditions mentioned above. Sketches developed from these expressions 
are found on Fig. II. 11, below). 

(lc) . Initial Values Problem Hodograph; Rotating Frame of Reference . 
The hodographs for this motion bears the same constraints and restrictions as 
did the displacement diagrams in paragraph (la) above. The general expressions 
for this case (with 7 = 0) are given by Eqs. (II. 22) or (11.29); consequently the 
restrictions imposed on this system will modify these equations accordingly. 

The general trace geometries which evolve here are found to be skewed 
ellipses. However, the non-secular cases are found to be ellipses symmetric 
to the coordinate axes. 

Following with the comments noted in paragraph (la), above, curves for 
these traces are not included below. This graphical information is also found in 
References [l] and [2]; however, for reference purposes, there are sketches 
in the data compendium attached to this section. 

(l d) , The Initial Values Problem Hodograph; Inertial Frame of Reference. 
Equation (11.36) is the general expression for the relative speeds; that ex- 
pression includes r as well as the initial state parameters. Of course, in this, 
the present case, 7=0. 

Hodographs for the initial values problem generally exhibit a divergent 
shape. Nevertheless, as indicated in paragraph (lb) this characteristic can be 
removed by setting (K^)j ~ 0. Since these traces are not too well known, a 
sample set is included below, for references purposes*. Interestingly, the 

*Also, see Reference [2] , and the data compendium at the end of this section. 
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FIG. II. lib. Sketch showing the out-of -plane hodographs (H r , Z f ;H f , Z') corres- 
ponding to the non-secular displacements on the preceding figure. 
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graphs exhibit a definite skewed appearance, for the parameter inputs used. 
(It is mentioned that these initial inputs were arbitrarily chosen as positive 
quantities; obviously, varying these terms would introduce a change in the 
geometry of these figures). 


2. Zero Initial -Values Problems. The problem types referred to here 
are a consequence of the forcing quantities, T and f . It should be recalled 
that these vectors represent specific (dimensionless) force systems which are 
aligned with the rotating and inertial triads, respectively. 


(a). In this first situation f = T(T t , T , T r ). The analytical expressions 

s ^7 £ 

illustrating the influence of this parameter on a relative motion are acquired 
from the same general equations noted earlier. As a result, the general state- 
ments set down below follow accordingly. 


First, the effect of T (alone) on the motion traces, that is both displace 
ment and hodograph traces, is one of a secular nature. There are no non- 
trivial force component relationships which can be introduced to remove this 
divergence. Likewise the curve forms which result here are not simple ones. 
Therefore, it is quite difficult to postulate, a priori, the shape of the various 
curves which are to be found*. (Incidentally, the comments above are not 
directed to either frame of reference, in particular; these generalizations 
hold for both). 


In view of the geometric complexities which arise for this class of 
problems, it is not advisable to attempt an in depth analysis of the traces. 
However, as an alternate, typical curves are included in the data section 
which follows this; also, some generalizations maybe gathered from a study 
of those figures. 


(b). 

vvv 


The specific force system aligned with the inertial reference triad, 

T , T ) has much the same influences on the relative motion as does 
-tl z 


*See the data compendium for representative sketches. 
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T. That is, the displacements and speeds are found to exhibit a secular 

character when either or both in-plane components are present „ Therefore, 

the only means available to nullify this influence is to eliminate both and r 

a H 

simultaneously. 

(As noted before, these comments apply to both coordinate reference 
systems; also, the trace geometries*, for both, are affected accordingly). 

From the discussions above it should be evident that these disturbance 
force systems lead to rather complicated relative motion figures. Consequently 
no concerted effort has been made, herein, to codify and classify them. Instead, 
as mentioned elsewhere in these paragraphs, a collection of representative 
sketches has been made; these are found in the data compendium. 

In the event that a reader would like a more in-depth discussion and des- 
cription of these traces he should consult Reference [2] . That document con- 
tains added details on these constructions. 

Summary . In the foregoing paragraphs the various out-of -plane traces have 
been broadly referred to. Generally, it is known that these figures are more 
complicated, in makeup, than the companion in-plane curves. Also, mathe- 
matically, the functional form for these graphs is more complex than that for 
the in-plane cases. However, when one studies these it is found that there is 
a geometric pattern present; and that some estimation of a particular trace 
may be made, a priori. Unfortunately, it would be rather cumbersome to 
attempt to develop a sense, or feeling, for these by means of written descriptions 
Consequently, the reader is advised to study the attached sketches and to ac- 
quaint himself with their overall particulars. 


*See the data compendium for typical traces. 


133 




X. DATA SUMMARY 


A Compendium of Data for the Out-of-Plane Relative Motion Cases . The in- 
formation presented in the following section is a condensation and illustration 
of the trace geometries which could appear as the out-of-plane displacement 
and hodograph diagrams. Since the case studies provided here include the full 
range of problem types considered in this investigation, there is a fairly large 
number of figures, etc. which appear. As an aid toward clarifying the arrange- 
ment of materials presented here, a brief indexing is included below: 

(a). The several subsections are arranged in the following order: 

X. Displacements referred to the Rotating Frame of Reference. 

2. Displacements referred to the Inertially Oriented Frame of 
Reference. 

3. Hodographs for the Rotating Reference Frame. 

4. Hodographs for the Inertial Reference Frame. 

Within each of the above subsections the material is arranged as follows : 

1. Initial Values Problems (T, Tj = 0). 

2. Zero Initial Values Problems (f ^ 0; Tj = 0; Initial relative 
motion state set to zero.) 

3. Cases for Non-Secular Displacements. 

4. The influence of Tj (alone) on the relative motion. (Some 
special situations are noted.) 

For all of the illustrations shown below, the input parameter (initial 
state values, T and Tj components) were arbitrarily selected as positive 
valued constants. In this regard the T{ typical trace geometries” will have 
characteristics related to this selection of quantities. However, it is not 
too difficult a task to examine the pertinent analytical expressions and ascertain 
the effect of (say) a particular negative valued constant. What does present 
difficulties is a determination of how the magnitude of the various inputs would 
Influence the traces. In many cases it would be necessary to examine the out- 
put, in detail, to find an adequate answer to these conditions. Nevertheless 
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one should be aware that the figures found herein are typical of these various 
Influences; and that these general shape characteristics are representative of 
the various cases studied. One final comment regarding the illustrations: in all 
instances the figures represent the relative motion for one orbit of the reference 
particle, only. Also, the various trace geometries are not reproduced here to 
a consistent scale; the primary motivation for this collection of data is to pre- 
sent, in concise format, illustrative information, primarily. 
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Out-of-Plane Displacements, for the 


Rotating Frame of Reference 
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Out-of -Plane Displacement Diagrams ; Rotating Frame of Reference 

The expressions defining these displacement diagrams are: 

(a) the in-plane equation, and (b) the third component relationship. 

(a) I 2 X «P) = [« 2 -2jjT 2 (^-) A a + [l 2 - f PB 2 jjK o + » T f. 

(b) J D X (<P) = J_ [ A cos<p+ JL f &in<p+T (1 -cos<P)l . 

3 3L O o _j 

Partial Solutions : 

1. Construction for Initial State Values (T - 0). 

The scalar displacement equations are: 


£ (<P)“ 2 (-A^ cos<p+A 2 sincp) +K^, 


T?((p) = 4(A l sin9+A 2 cos^) - - <pK t 


C (</>) = £ cosp+^'shup ; 
o o 


wherein 


A, 1; A„ B i'l; 

1 2 L 1 1 oJ 2 2l 2 2 OJ 

K s 2 f 2 J A -J B i'l; K =J A -2JB i ' . 
1 L 1 o 1 2 oJ 2 2o 22o 
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Out-of-Plane Displacement Diagrams; Rotating Frame of Reference (c 


Quadric equations for these planar traces are: 


(a) For the 


\ -plane; 


+ 2 [ (A 1 + C o > " (A 2 + C o 2) ] 008 2 ^ + [Co C 0 - A 1 A 2 ] Sin 2<P 

This trace is a skewed ellipse with its geometric center not at the 
coordinate origin. The principal axes of the figure do not coincide 


with the coordinate axes . 


(b) For the 


(-plane; 


2. Traces for the Zero Initial State {A, A ^ - 0 ; T = t(t^, T y 

are defined from: 


4 (<p)~T t (1 -cos<p) + 2T to-sin<p), 
s V 

3 <p 2 

7?(^) =s 4T ?? (l-cos(p)-2r^ (<p-smp)--=- , 

C^P) = T^(1 -cosg>). 




(la). O is the initial locus of the trace; 

P is a general position on the curve l 



(lb). This trace is not a closed curve due to the 
secular effect of ??. The figure produced is 
a "moving" ellipse, skewed with respect 
to the coordinate axes. Note that the secular 
influence vanishes if K = 0. 
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Out- of- Plane Displacement Diagrams; Rotating Frame of Reference (continued) 


Equations for the Planar Geometries 

(a) Traces for (T ^ , T J) only . (7 components > 0) defined from: 
£(<P) = T£ (1-cosP), 


hence, 

(a,l) 


r}(<p) = - 2T^ (<p- simp), 
C(<p) = T^(l-cos<p); 


Lei _ iM . and izM) «_ 


t 4 


'c 

(a, 2) 


2T 


4 


(<p-sin<p), 


- SS 1 -COS(p. 

T c 


(b) Traces for (r , T w ) only. (T components > 

1 — t? — c — 

£(P) = 2r^(<P-sin<p), 

T?(P) = 4T^[(l-cos<p)-^] , 

Cm-T^Q. -cos<p); 


0) 


hence, 


and 


2*L 


<p- sin<p, 


M s JLtel . 


27. 


; or, 


and 


V 


iM = l-cos 

v 4r 8 r c 



(a, 1). This trace is: 

A line, originating at 
the origin; frequency 
of motion is equal to 
orbit frequency. 

(a, 2). This trace is 



% ) 
O 



A cycloid, originating at 
the origin, incrementing 
in the negative 77 -direction. 


n 

o 


o 


<M> 


(b,2) 


This is a: 

(b,l). Cycloidal trace, beginning 
at origin; its frequency equals 
the orbit frequency. 



(b, 2). A positive sloped line, with 
an accelerated incrementing locus,! 
in the negative ^-direction. 
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Out -of -Plane Displacement Diagrams; Rotating Frame of Reference (continued) 


(c) 


Displacement Geometries for (T ^, T , 


J c 1- 


Note: The resulting traces (here) are combinations of 

the effects exhibited in (a) and (b) above. 


On the (£ , C) plane, the displacement curve is described as: 
(c,l) - — Q (<p- sin<p) = . 

T z T i r ? 

The straight line trace here is modified by the cycloid 

(due to T ); consequently, the £ -symmetry of the cycloid is 

V 

skewed (as shown). 

On the (r} t £) plane, the displacement curve is described by: 
<°- 2 > |^ = -l + 


i 



( 0 , 1 ). 



Here, the cycloid (from 7^ ) is modified (by the line- 
with-accelerating secular effects) to produce the complex 


geometry shown at right. 


Relative motion displacements, traced for one orbit, 
illustrating a typical situation as produced by a fixed 
specific force, T. 




(continued) 
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Out -of- Plane Displacement Diagrams; Rotating Frame of Reference (c 


Quadric Equations for the Planar Traces . 


(a) For the 


on the 


+ {-“ 2 - + C £'} sin2(p. 
t 4 o oi 


(b) For the 


on the 


\ f — \ L-t- mAMrn * 

*{ C o E o'“^} 




(a). The (4, C) trace is a skewed ellipse, as shown. 
Its frequency of motion is the same as the base 
(circular) orbit. If the motion is constrained to 
pass through the origin {A =0), the ellipse re- 
2*i ‘ duces to the line shown on tfiis plane. 


c 


(b). The non-secular (rj, £) trace is a skewed ellipse, 
beginning at O, with a displaced center (K 2 ). If 
the motion is further constrained to pass through 
the origin {Jt 2 0), the skewed ellipse is replaced 
by the symmetric one, beginning at the origin. This 
ellipse has its center located at 
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Out-of- Plane Displacement Diagrams; Rotating Frame of Reference ; Influence. 
The general expressions describing the displacements are: 

(a) the in-plane equation, and (b) the third component expression. 

(a) I 2 iCP) = [l 2 ^(J 2 -J 1 )]T 2 ( V -)[A a ] Lv/ [l 2 -|«pB 2 jjK o + $ T f I 

(b) . J 3 I ^) = J 3 [^ 0 cos<p+A^sin^+T I (1 -cos<P)j . 

Partial Solutions 

1. The only partial solutions which will provide new information 

are those due to the $ -term(s). Only traces involving r are discussed 

T I 

below. 

2. Traces for Zero Initial State (£ = 2 1 =0; 7 T =T r (7L , T„, r )), * 

o o I I a Hz 

are defined from: ^ 

-2r s sin<p- (1 -cosp) J + £ S i n <p-<p C os<p J , 

77 ((p)~ 3 T s Jjp(l+cos<p) -2 sin<p] + 5 t h sin<p- (1 -cos<p)J, 

C(p) = T (l-cos<p). 

Equations for the Planar Geometries 

(a) Traces for (r . T ) only (T components > 0), defined from: 

^(<p) = 2T H j"^psin<p-(l-cos<p)], 

??(<p) = 37L <p(l +cos<p) -2 sin<p I, 

L, —I 

£(<P) = T z (1-cosfp); 

hence, + sin<p, and <p =2 ((p- sincp). 

£ Tm I 4 o | H f 


♦■The motion traces originate at the coordinate origin. 





2a. Example of out-of-plane traces, 
illustrating effects of Tg , T as 
applied force (s). z 
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Out-of~ Plane Displacement Diagrams; Rotating Frame of Reference (continued) 

(b) Traces for (T,„, T ) only (T components >0), defined from: 

H z 

3 T 


£(<p) = — [sin^-<pcos(p] , 
r?(^)=5T H [^psin<p-(l-cosp) j , 


(2b). Traces for r , 

_ ri 


hence, 


and 


C((p) = T (l-cos<p); 
z 


H z 


UM i JLM = 3,3 


5r. 


sin 


H 


(c) Displacement Geometries for T . 

Note: These traces are described by the 

(2c). Sketch of combined effects, 
summing of coordinates in (a), and <j ue to T 

(b). above. 

Sketches depicting the trend, for the present 
situation, are included at right. 





C 



Out - of- Plane Displacements, for the 
Inertially Oriented Frame of Reference 
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Out-of~ Plane Displacement Diagrams; Inertial Frame of Reference 


These general solution equations are: 

V^K 3 ( J 2 - j l> +T 2 ^ "> ]\ +T 2 & ->{ [V - I ' * <*2 J 1 ) ] K 0j + * T T } 

J 3 ^ (<P) = J 3 {^ Q C0S(p+B^ sin<p + f(l -cos<p)j= . 


Partial Solutions, 


1. For the Initial State Values (T=0) Problem ; the displacements are: 
S(<p) +3A = (A cos 2<p-A sin2<p) + (K costp-K sincp) + ^ K sm(p, 

H(p) - 3A 2 = (A x sin2<p+A 2 cos 2 <p) + (K^ sin<p+K 2 cosp) - -p cos<p, 


Z((p) = Z cos<P+Z' sinp; 

^ o o 

wherein 

i— 1 A„ si |k -J ft 1, A +(J BJft r 1 ; 

it 1 2L1 1 oj’ 2 2L2 o l 2 2 ; oJ’ 

CD 

K =2[j ft -J B ft r 1, K ■ - fj ft +2J B rt ft* 1 
1 L 1 o 1 2 oJ 2 L2o 22 oJ 

Quadric equations for the planar traces are : 

(a) For the (E, Z) -plane . 

/ 3K 1 , \2 2 2 2 .2 ,2 2 2 

^E+3A 1 <psm<p J +Z =(A 2 +K 2 +Z^ ) + (A 1 -A 2 )cos 2 <p 

+ (K^-K^+Z 2 -Z ,2 )cos%-A 1 A- sin4<p+(Z Z r -K K )sin2 <p 
1 2 0 o 12 o o 12 

+2^A K cos 2pcosP- A^K 2 cos 2(psin<p- A 2 K^sin2^cos^ 

+ A K sin2<psincp „ 

2 j Li rnJ 



A typical 
Result 


(la). This trace is secular (though appearing to be 
closed). Complex nature of the quadric 
equation does not lend to generalization and 
individual parameter influence descriptions. 


Note: O locates the initial position for each trace. 



Out -of- Plane Displacement Diagrams; Inertial Frame of Reference (c ont inued) 


-Cl 


(b) For the (H, Z) -plane . 

3K 

(h - 3A 2 + <pcos <p ) 2 + z 2 = (A 2 + K 2 + Z’°) + (A® - A 2 ) co S 2 2p 

+ (K - K 2 + Z 2 - z' 2 ) cos V A, A, B ln4* + (K, K + Z Z')sin2<p 
" J- o O 1 C 12 O o 

+ 2 [A K sin2<psinp-f-A K sin2<pcos<p+ A K cos 2<ps in <p 

^ X J- i Z Li 1 


+ A 2 K 2 cos 2(p cos <p j . 

2. The Zero Initial State Solution 

$ o = *'-0; r = r(r^,r^,T^))ls expressed by: 

(H ((p) -f j ) = r^ ^cos9+ ~ cos 2<p J- r^^4sin<p ~sin2(p 

+ 2<p sinp+T^ £cosp+ ^psLn^pJ|, 

(H(<p) +3r ^)“ r ^ [sin<p+ j sin2cp J cos 9 -cos 2<p J 

2 ( p a [r r? [sincp*- ^ cos tp J - cos<p| , 


+ 



(Z(<p) - r^) ** - cos <p, 


(lb), A secular, open trace largely complicated and 
influenced by the Initial Values present. 


Note: O is the initial locus for each trace. 



Out -of- Plane Displacement Traces; Inertial Frame of Reference (continued) 


Equations of the Planar Traces : 

(a) The Influence of (T ^, Tg ) only (with T components > 0), 
(H (<p) + ~ [cosp+ - cos 2<p +2 <p sin<p J 

H(<p) = |jsin<p + ~ sin2<p - 2<pcos<p J 



(Z(p)-T^» - cos (p . 
Hence the describing expressions: 


(2a, b). 


Z - Tv 






and H + 2T^ (<p sin<p)-r^sin<p 


(b) Traces for (T , T^ ) only (T components > 0). 

S flp) = - t £4sinp-sin2p]+2T^£cosp+^slnp J<p 
H(<p) + 3T^“ [jl cos <p - cos 2<p J+2r^|j3in(p - ~~ cos<p 
Z(cp) - T^-~ Tg cos <p , 

(c) Combined effect of (T^ , T , T^) would be represented by 
the additive influences shown above . The general geometries would be 
amplified, in divergence, as is easily seen here. 


The secular nature of the describing equations 
appear to indicate closed curves ; however, the 
divergence is evident (by inspection). 



(2a, b). The divergent nature of these traces is quite 
evident from the descriptive expressions. 



Out -of -Plane Displacement Diagrams, Inertial Frame of Reference (continued) 

3. Traces, with Secular Terms Removed; Initial Value Solution (T = 0), 

Note: To eliminate the divergent (secular) nature of these traces it is 

only necessary to nullify , where = 2 [j^ ] . 

Obviously this can be accomodated by either: (1) letting 

j ft == j B„ ft * ~ 0; or, (2) by setting J ft = J B ft r . 

1 o 1 2 o 1 o 1 2 o 

(a) For illustration the case depicted below merely sets K = 0 

without stipulating how this was accomplished. As a consequence: 

A = - - (J ft); K - 0, but K and A remain as before; and, 

1 2 1 o' 1 2 2 

S (<p) - |s = -|[h o cos2<p+(H o +H^)sin2(p)].+ (H Q + 2Hp sin<p 

H(cp) - — (H +S r )=~-f"-H sin 2<p+ (H + S 1 ) cos 2<p 1- (H + 2E 1 ) cos<p 
v ' 2 ' o o 2 L o o o J o o 



Z(c p) “ Z cos<p+Z T sin<p. 
o o 

In general, these are not easy expressions to discuss, geometrically; 
they do represent closed curves on the coordinate planes; but they are nominally 
complex (of single and double frequency), with an initial point not at the 
coordinate origin. 


Cases 3(a) and 3(b) are shown, sketched 
on the (E, Z)-plane. Case (b) describes 
a non-secular situation; motion is initiated 
at the origin. 

(A typical situation) 
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Out-of-Plane Displacement Diagrams, Inertial Frame of Reference (continued) 

(b) If, In addition to the elimination of the secular influence, the motion 
commences at the coordinate origin (ft = 0), then the trace equations reduce to: 


S ( 0 ) =E ' sin<p [2 - cos tp 3 


H (0 ) - H f = cos<p [cos<p - 2] 


Z(<p) = Z^ sin<p ; 

or, the descriptive equations for the planar traces: 


H 



3. Typical displacement traces for non- 
secular state conditions. 



TST 


Out -of- Plane Displacement Diagrams; Inertial Frame of Reference; Influence, 
General expressions describing the displacements are: 

(a) the in-plane equation, and (b) the third component expression. 


(a) IgRftJ) 


= {T 2 C!(9-) + 3[j 2 -J 1 ]}[A ] 1>Vt+ T 2 fl»-)[l 2 -f?>(B 2 J 1 )]K 0i + T 2 fi9-)[* T r l ] f 


and (b) J 3 B(<p) = J 3 [® Q cos ^ ^ sinfP + Tj (1 -cos<P)J . 

Partial Solutions . 

1. The only partial solutions which provide new information are 

those due to the $ -terms. Only traces involving T are discussed below. 

T _ Jl 

2. Traces for Zero-Initial State (B = B f =0; r T =T t (T w , T tt , T )*) , 
— — o o I 1 * H z 

are defined from: 

H fto ) ~ [2(1 -cos<p)+2(l -cos 2<p) ~3<p(sin<p+ j sin2<p)J 
+ TH^sintp- j sin 2<p- ^ (3 - cos 2<p) j , 

H fP) = Tg Ff (4 cos <p + 3 + cos2<p) -2(sin<p+sin2<p)J 
+ T h[ 5 ^ ~ cos ^) “ 4 (1 -cos 2 <P) + ^ sin 2 <p] , 


and 


z(<p) = r (i-cos<p). 

Z 


Equations for Planar Geometries. 

(a) Traces for (7L , T. ) only (T components > 0), defined from: 

A 'Z 

H(<p) =T,_, [2(1 -cos(P)+2(l-cos2^)-3<p(sin<p+ jsin2(p) j, 


H((p) = T w ^ (3 +4cos<p+cos2<p) ~2(sin<p + sin2<p) J 


1 

1 


v 

V 

Vs. 


o 



and Z((p) -T (1 - cosp); 
z 

hence, the expressions: 

(continued on next page) ^ a * Examples of trace geometry, on the displace 

*The motion traces originate at the coordinate origin. ment planes (H, z) and (H, Z) due to applied 

force (s) TL, T . 

H z 



Out -of- Plane Displacement Diagrams; Inertial Frame of Reference (continued) 


(a) Traces for (T M , T ) only (continued) 

Hz z 

+ 4 ( 1 --^^)]== - 2(sin<p+sin2<p). 

H z z 

(b) Traces for (T„, r ) only (T components > 0), defined from: 

z 

_ r 7 12 1 

*(<p) = T^. [J 5 sln<p - — sin2<p-3^(l - ~ cos 90) J , 

H(9 ) = T h |j 5 (1 - cos 9 ?) - (1 -cos 2 (0) + ^ sin2<p J, 

and Z(p) =T (l-cos(0); 

w z 

ci 

M hence, the expressions: 

H z 


(c) Displacement Geometries for T 
Note: These traces are obtained by summing 
coordinates from (a), and (b). above . 



2b. Trace geometries illustrating effects due to 


H* 



Out-of-Plane Hodographs, for the 
Rotating Frame of Reference 
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Out-of- Plane Hodograph Diagrams; Rotating Frame of Reference 


The general expressions are: 

*2 *■ ' *»> - k + 3 (J 2 - J l> >2 T 2 &>"> A a • I [ B 2 J J K o + V • 

and J 3 i 1 ((p)=J^A^cos(p+ (T-A^)sm(p J . 

Partial Solutions , 

1. Traces for the Initial State Solution (T = 0), may be acquired from: 
£ f ft0) = 2(A sincp+A cosp) 

1 2 Q 

7? r (<p) =4(A 1 cos<p- A 2 sin<p) - - ^ 

C f (<p) = C'cosv>-C sin<p. 
o o 

wherein 



A i 'i[ K r J i A J ; VitWoJ 

K ! " 2 K A o- J 1 B 2 A o) V^V^VVo]' 


(l,a). This trace is an ellipse, symmetrically disposed 
about the origin. Its principal geometric axes do 
not align with the coordinate axes. 


Quadric Equations for the Planar Traces . 

(a) Trace for the (£ r , £ T )-plane . 

* I { [*2 * C / ]' ' [ A l *' <£]} “> W* [ A i A 2 - C„ tj ■' f » . 

(b) Trace for the (r)\ C^-plane . 

[^*?S] 2 *«W4{[ A XK A ^,f]} 


n* 



(l,b). The trace here is an ellipse with its center displaced 
3 

to rj = - -r K . Note, this ellipse is also skewed. 

C ci JL 
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Out-of-Plane Hodograph Diagrams; Rotating Frame of Reference (continued) 


2. Trace for Zero Initial State (A- -A * -0: r = r(T J . 
o o £ 

The scalar expressions are: 

£'(<p) = 2T^(l - cos <p) (sin<p) 

r? , («p) = -2 (1 -cosp) + 4T^(sin<p) -3<pT^ 

Z'(<P) = T^ (sin<P). 

E quations for the Planar Geometries . 

$&) Traces for (T ^ , ) only, (T components > 0) 

are acquired from: 

^((p) = T^sintp t 



(2a). \ 

This trace is a line, symmetric about 
£ T axis, with frequency matching that 
of circular orbit. 



hence 



defined by: 

4 , (^) = 2T^(1-cos<P), 
7^) = 4r^(sin<p)~3<PT^, 
C((p) = T^sin<p; 
hence, the quadrics: 


(2b). This trace is an ellipse, with 
center at £’ c = 2 T^. Motion’s 
frequency is the same as that 
of the base orbit. 


v £ 


= 1; and, 


£M + 3, = £M 

4 r 4 r r„ ' 

n £ 


(2a). The trace here is: 

An ellipse with a displaced center 
(r? c = -2T^). Motion on the ellipse 
has a frequency equal to the orbit 
frequency. 



(2b). This S ^shaped trace is produced 
by a ’line” modified by the secular 
^-influence, (3<pT ). 

V 
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Out-of- Plane Hociograph Diagrams ; Rotating Frame of Reference (continued) 


(c) The Hod ograph Geometries, for ( Tg, T , r» ) combined . 

Note: These traces are due to the combined 

effects shown in (a), and (b). above. 

On the (£', £')-plane the displacement curve 
can be expressed by the quadric: 

(0,1). [4 * fio) - 2r T? ] 2 + [C T (<P ) n 2 = (r| + r|) [4r^- (r^ + ) ] (1 + oos 2(o ) 

-2 T^sin2<p. 

Here, the ellipse (due to r , T^) and the line (from 
the -effect) combine to produce the skewed ellipse shown 
at right. 

.On the (?7 r , ^ T ) -plane, the equation of the curve seen 
there can be written as: 

(c,2). [^ £I + |<p] = ^ £} -^(l-oos(p). 

Consequently, the ellipse (due to , T^) is coupled 
with the s -shaped curve (from r ) leading to the trace shown 
at the right. 

3. Traces for Secular Terms Eliminated ; a Modified 
Initial Value Solution (T s 0) 

Note: It has been noted that the secular (divergent) 

influence can be eliminated by setting = 0 
(see discussion with displacement traces). 



(c.l). 



Typical hodograph traces developed 
for the fixed thrust, T. 



Out -of- Plane Hodograph Diagrams; Rotating Frame of Reference (continued) 

(a) When vanishes it is found that one set of 
hodograph equations is: 

£* (<P ) = ^ cos<p-^ o sin<p, 

Tf(<p) - - 2 [£^sin<p- £ o cos(p ] , 

C cos(p-C 0 sin<p; 

these lead to the following quadric hodograph expressions: \ 

\Both traces (here) suggest ellipses, but these figures have their 
principal axes not parallel to the coordinate axes (i.e. , skewed figures). 

(b) For the added specialization of motion commencing at the 
origin (A s 0), the (above) equations reduce to the following set: 

£ f ftp) £'(0) 

J (a line passing through the hodograph *s origin) 


-1 (an ellipse, aligned with the hodograph’s axes) 


p 



O is the initial locus on these planes. 


7 )’ 



Typical traces depicting a non-secular 
solution. 



tn 
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Out- of -Plane Hodograph Diagrams; Rotating Frame of Reference. 

The general solution equations are; (a) for the in-plane components, and 
(b) for the third component: 

(a) I 2 I •*» - [l 2 + 3(J 2 - J x )] B 2 T 2 (P-)[aJ. v _ - f [b 2 jJvVl , 

and (b) J f .A , (np)= Ji 1 cosp+ (T - -£o) sin<p | . 

3 »• o 1 -i 

Partial Solutions . 

1. Only the partial solutions expressed in terms of the T 
components provide new information* Consequently, only these 
are described below. 

2. Traces for Zero Initial State (A = A 1 =0; T, = TAT* , T_, r ))*, 

O O 1 1 » li z 

are described by: ^ 

a jj^cos^-sin^ J + £<psin<p J , 

r)'(<p) = 3T v [a -cos<p) -psinpj + 2 cos0-sin<p J , 

and C T (<p) = r sinp. 

Equations for the Planar Geometries 

(a) Traces for (T_ , 7 ) only, (T components > 0). 

— — *-E z 

£ '(<P) = T T w [spcosp-sinpj , 

T7'flP) = 3T~ [(1 -eos<p) -<psin(p J ,and C f (^) = T z sin<p; 
hence, the parametric equations: 

+ i!M = 3 p C0S9; 

Tw T 

and ^ ^ + 5^21 p = (l - cos <p ) . 

3 Tb T 

~ z 



2a. Typical hodograph traces due to T v and T 


*The motion traces originate at the coordinate origin. 
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Note: These traces are the composite 
of those diagrams described for 
(a), and (b). above. 


(continued) 



2b. Typical traces for the influence of T . 



Out-of-Plane Hodographs, for the 
Inert ially Oriented Frame of Reference 
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Out-of-Plane Hodograph Traces; Inertial Frame of Reference 
General equations for this solution type are: 

I 2 a 1 «P ) = B 2 { [ 2T 2 (2^ -) ] A &I + T 2 (<P -) t 2 -K -.|( B 2 - J i)^] K o i 

+ t 2 (p _ ) * t t} + t 2 (<p~)Wt 
J R r (<p)~J {ft* cob¥ + (T-S )sin<p} . 

0 O^O O v 

Partial Solutions 

1. Traces for the Initial State Problem (T - Q) , are obtained from: 
H 1 (cp ) = — 2 [a sin2<p+A cos2<p 1+^-K sin<p-K cos<p + “^K cos<p, 

H’fto) =2 [a cos 2p-A sin2p 1- ^ K cos<p-K sinp+ ^ K sintp, 

UX Z -JZX z z x 

Z'(<P) = Z^cos<p- Z Q Sin<p; 
wherein: 

A. =~Tk, - J B 1 ; A ^[k +J bJ'1 

1 2 L 1 loJ 2 2 L 2 2 2 oJ 

K s 2 fj fc &*1; K rt =-fj Q B +2J B ft']. 

1 Llo 1 2 o J 2 L 2 o 22 oJ 

Quadric Equations for the Planar Traces. 

(a) For a trace on the (S r , Z f )-plane . 

2 . 1 t ,2 , _2. , . ..2 .2 


aa). 


(H , -|5k i cos(0) 2 + (Z') 2 = (4A^ + |k 2 +Z 2 )+4(A 2 -A^)cos 2 2<p 

+ (K l - jK 2 +z; 2 -Z 2 )oos 2 V5+4A 1 A 2 sin4^- (1^ + Z Q Zj i )stn2?J 



These traces are necessarily secular, and 
complex, due to the make-up of the scalar 
equations forming the quadrics. This com- 
plexity eliminates the generalizing and 
prediction otherwise expected. 


n 1 

- 4 -A. K sin2<psin<p-A K sin2^cosp+— A K cos 2<Psin<p-A 9 K cos 2<pcos<p . 
L2 11 1Z ZZX -4Z _ 


Note: O is the initial point for each trace. 
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Out-of-Plane Hodograph Trace , Inertial Frame of Reference (continued) 


(b) For the trace on the (H 1 , Z*)-plane 


(H* - K sin<p) 2 + (Z’) 2 * (4A 2 +K 2 + Z 2 )+4(A 2 -A 2 )cos 2 2<p 


2 1 


2 2 o 1 2 


+ |" 7 K 2 - K 2 + Z ,2 - Z 2 lcos 2 ^)+ (— K, K - Z Z r )sin2p 
L41 2 o oJ 2 12 oo 

- 4A A sln4<p~4r“A K cos 2<pcos<p+A.,K cos 2<psin<p 

12 L 2 1 1 12 

- K sin2<pcos<p-A K sin2<psinp 1 . 

U Ca X Z Z J 

2. Trace for Zero Initial State (£^ = B^ = 0; T=T(T^, T , T^)), 

S , (^) = T i .(sin<p-sin2<p) -2T (cos<p - cos 2<p) +<p 2T t cosp 
s V L 9 

+ (sin<p+ -jp cos<p) J . 

H'to) = - (oos^- cos ZP) - 2r n (sin^- sin 2^) + * [s^Bln* 

- r^(cosp- sinp) J . 

Z , (<p) = Tg sin<p . 

Note: O is the initial point for each trace. 


H* 



(lb). This secular curve passes through 
the initial point at = 2tt - (this is 
due to the parametric form of the 
expressions). 
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The secular influence of these expressions 
requires a divergent geometry. (Combined 
effects of 7 >j, 7> would be provided by 
adding displacements on the curves). 


S f ( 9 ) «- r |j2( cos <p- cos 2<p) -<P(sin<p+ cosp) J 
H , {^) = -r^^2(sin ( P*-sin2<p) +(p(cos^)- sin<p) J 


Z 1 ^) = sin<p. 

(c) The combined effects of the full three components 
of the specific force T is obtained by adding ordinates on the traces 
depicted here. Due to the close similarity of the curves it is quite 
apparent that the geometries would not significantly be altered by 
this composition. 

Note: O is the initial position for each trace geometry. 



The additive influence of total thrusting is found 
by combining the ordinates of the curves. 
Secular influences are predominant here. 
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Note: The general requirement for removal of secular effects was 

found to be satisfied by nullifying KL ( = 2 [j ft - J, B_ ft r ] ). 

1 1 o 1 2 o 

It has been noted that this may be carried out in either of two ways. 

(a) Again, for illustration purposes = 0, solely; then A^ ss -“ , J ft^; 
and, consequently: 

H f (cp) = ^H Q stn2 <p- (H q +H^) cos 2<p J + (H q + 2S^) cos<p 
H r (^) = - THq cos 2<p+(H o + H^) S in2<p] + (H o + 2H^)sin<p 



Z'toi^Z’coscp-Z sin <p. 
o o 

As before, these expressions are not readily visualized for their 
geometric interpretation. These equations represent closed curves 
involving functions with single and double frequencies. 


3. Typical traces produced for the 
general and special non-secular 
state conditions. 
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Out -of -Plane Hodograph Diagrams, Inertial Frame of Reference (continued) 


(b) When these non -diver gent relations are further specialized 
to pass through the origin, the equations reduce to: 


H'ftO) = H ^ + 2 cos <p (1 - cos<p)] 

H T (p) -E ^ [2 sintp (1-cosCP)] 


Z'ftfl) = Z^ cos 43 


and, the representative planar trace equations are recast as: 




and 


|^ = 2 ( 1 .^M) sin¥) . 


O 



3. Typical hodograph traces for the non 
secular state conditions noted in (a) 
and (b), at left. 



I - 1 

as 

as 


Out-of-Plane Hodograph Diagrams; Inertial Frame of Reference . 

The general solution equations are: (a) for the in-plane components, 
and (b) for the third component. 

(a) !/'«<» = B 2 {2 [t 2 (2(0 ")][ A a I ]i. v . + i T 2 fi9 ‘, [ 2 J 2 " J 1 <Vl> ] 

and (b) cos<p-^ snKp+TjSin^ j- . 

Partial Solutions . 

X. Only those partial solutions expressed in terms of the f 
components provide new information. Consequently, only these are 
described below, 

2, Traces for Zero-Initial State (ft =ft f =0; s Tj(T g , , T ))*, 
are described by: 

H *(<p) = t_ sin2<p-sin<p-3<p(cos^+ ~ cos 2(p) j 

+ T 1*5 cos (p- 7^* cos 2 p- sin2<p- 7 1 , 

HL 4 2 4 J 

H t (^) = T_ [coscp+ “ (9-13 cos2<p)~3<p(sin<p + ~ sin 2<P) J 
+ T^jjjsinp- — sin2p+^p cos2<P j , 


and Z f (p) = T sinp. 

z 

Equations for the Planar Geometries . 

(a) Traces for (T„, T ) only (T components >0), defined from: 
H 1 (cp) = [sinp(y- cos<p-l)-3<p (cos<p+~ cos 2<p ) J , 

H T (p) = 77 f~ (1 - cos 2^0) - (1 - cos p) -3psin9(l + cos9) 1 , 

H 4 — 




and Z'(<p) = T sin<p; 

z 

♦The traces originate at the coordinate origin. 


2a. Hodographs illustrating the trace geometry due 
to 77 , r . 

a z 
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Out-of- Plane Hodograph Diagrams; Inertial Frame of Reference (continued) 


(a) Traces for (T ,_T ) only (continued) 

hence the parametric equations: 

- . (£ . ! + Z ) (1 - — cos<p) = -3<p(cos<p+ ~ cos 2<p) , 
T_ 7* 2 2 

* z 


and 


and 


_LJ + 3( p ^ + cos(0) = “ (1 - cos 2 <p) - (1 -cosp). 

s z 

(b) Traces for (T „ , T ) only (r components > 0). 

ni H z *i 

H*to) “T ~~ (l-cos2<p)-5(l-cos<p)-3(0(sinpcos<p) , 
H L- 4 -J 

H, ( ( P) =T iif? cos 2(0+5 sin p- ~ sin<pcos<p 1 , 

HL 2 2 J 

Z T (<p) = 7 sin<p; 
z 


hence, the parametric equations 
s » 


^ + 3<p Z ^ cos<p= (1 -cos2<p) -5(1 -cosp), 


H 


and M + Mm coaM ).a col! , i . 

T tt T > 2 / 2 

H z 

(c) Hodograph Geometries for T . 

Note: These traces are defined as the summation 
of coordinates from (a), and (b). above. 



2b. Typical traces illustrating the influence of T^> T 
on the inertially described hodographs. " J 



XI. APPLICATIONS AND SELECTED EXAMPLES 


Introduction. This section of the report will be given to the development, des- 
cription and solution of some selected examples recognized as problems in 
relative motion; and, to a few application situations. In particular the problem 
of intercept for a ’trusting" particle is described, but for the special case of 
a fixed magnitude thrust aligned with the inertial triadic reference system. 
Secondly, some generalizations of the point-to-point transfer problem will be 
included, and descriptive solutions will be found. A third topic included here 
is that of defining appropriate "look-angles '* for the relatively moving particles. 
Since these angles may be inferred for both representative coordinate frames, 
it will be prudent to present a consistent analysis for the two systems. 

Basically, the information found in this section is addressed to the pro- 
blem areas set down above. However, in addition to this there will be found a 
few illustrative cases noted which point to the utility of the results. Unfortu- 
nately time and space limitations do not permit for a more informative and en- 
compassing expose of the analytical (and numerical) capabilities of this work. 

It will be left to the interested reader to describe and develop solutions, using 
the current results, which will evolve from his own applications needs. 

Look Angles . The concept of "look angles", as the name applies here,, is illus- 
trated in the sketch below. On the figure these angles are referred to the local 
(horizon) frame; however, it must be recognized that the same geometry is 
applicable in the inert ially oriented coordinate system also. 

The two angles shown were chosen because: (1) they represent a con- 
sistent nomenclature; and, (2) they are relatively simply to describe mathe- 
matically. Here the angle, a, is referred to as the "cone" angle, while B is 
called the "clock" angle. It should be mentioned that the angular range for 
these two quantities will be: 
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and 


0 ^ a £ 77, 

-77 ^ 0 ^ + 77. 


In the illustration P is the reference (circular orbiting) particle while 
Q is the relatively moving orb iter. Hence the angle pair (a, 0) locates Q re 
lative to P — this is consistent with "location” as defined by the relative posi 

tion coordinates (x, y, z). 

(Note: With the two 



reference frames assumed 
to be coincident, initially, 
there is no initial angular 
separation between them; 
therefore, the descriptions 
for a, 6 will parallel one 
another in both reference 
schemes). 

It is schematically 
apparent that a plays the 
role of a "cone angle", as 
its name infers. However, 


it should be noted that when 


the coordinate y < 0, then 

Fig. 11.12. Sketch depicting the position a> 

angles (a, 0); a = cone angle, 0 = clock 

angle; both angles are positive, as shown. ” c l oc k angle", 

0, will range over a full 
circle (277); also, its defi- 
nition (below) is such that 
0 is an angle consistent 

with the algebraic sign of the coordinates (x, z). In addition, one should 
recognize that 0 is measured in a plane whose normal is the unit vector, e . 
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Taking account of the above information it is seen that when the relative 
coordinates vanish, identically, these angles are undefined. 

Cone and Clock Angles Defined. From the sketch one sees that the relative posi- 
tion vector (r) can be described as: 

r = p + ye , 

where y 

p = xe + z e , 
x z 

(and similarily for the inertially oriented reference system). 

Making use of the appropriate vector and scalar products one sees that 
the circular functions are ascertainable from the following operations: 

r x e ~ sin a, 

y 

pxe ~ sin 0, 
x 

p • e ~ cos 6, 

and x 

p * e ~ cos a ; 

y 

J 

therefore, a description of these quantities can be given by: i 


0 = tan ^ (z/x), 

and 

a -tan (|p|/y); 


where IpI is the magnitude of the vector in Eq, (II. 114b) « 

Expressed in non-dimensional form, and for corresponding frame (s) 
of reference, these angle definitions maybe recast as: 

(a) for the local rotating (horizon) frame: 

tan 0= £/£, 

and 

tan ct=Ji 2 + Z 2 /r); 


(II. 114a) 
(II. 114b) 


(II. 115a) 
(II. 115b) 


(II.U6a) 
(II. 116b) 


170 



(b) for the inertially oriented frame: 

tan 0 - Z/S , 

and 

tan d= J S 2 '+ Z 2 /H. 

The definitions stated above indicate where an observer in ”P” should 
look — relative to his own frame of reference — to ’’locate 1 ' Q. Logically, one 
would (as a first guess) expect that for the converse situation the angles should 
be incremented by 77 radians. Unfortunately, this is not quite the case; a 
’’correction” should be added (or accounted for) before a true definition can be 
had. For instance, true imagery of these angles will exist only when the inertial 
reference frame is used for observational information; and, only then if the 
frames in ”P” and "Q” are aligned. When these conditions exist, the corres- 
ponding angles, from the two particles, are phased at 77 radians. 

A Correction to Look Angles. For the ’’local” rotating frame of reference the 
clock angles are phased as noted above; however, the cone angles cannot be so 
simply related due to an angular rotation needed in locating one local triad re- 
lative to another (see Fig. 11.13, below). 

Assuming that only the in-plane angularity needs the correction, then it 
is easily shown that the reference triads are connected by the transformation: 


e 

X 


cos A <p 

sin A(p 

0 


e 

X 

e 

y 

= 

-sin A <p 

cos A <p 

0 


e 

y 

e 

z 

Q 

0 

0 

1 


e 

z 

% * 


In this expression the subscripts ”Q” and ”P” signify an appropriate origin 
for the triads. Recognizing that the relative position vectors satisfy the relation, 
r = - r, then the positional coordinates, for P relative to Q are given by: 


(11.116c) 
(II. 116d) 


(11.117) 
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y s [ x r y i’ z j > 

then 


X 1 


cos A<p sin A<p 0 


■ « 

X 

y l 

= - 

-sin A<P cos A(p 0 


y 

v 


0 0 1 


z 



► 


• — 


(11.118) 


(An evaluation for A <p can be approximated 
by the relationship: 

A<p — tan 


For most application purposes, consistent 
with the "smallness” of the relative dis- 
placements, this is a valid approximation). 


Fig. 11.13. Sketch, illustrating the 
in-plane angular separation be- 
tween two local triads of refer- 
ence. Note that the relative 
position vector, locating "P" 
from would satisfy the re- 
lationship: r = - r. 


Making use of the above transforma- 
tion, and denoting the new position angles by 
(a , 0 ) # then it ls evident that: 

e i -‘“"'(W. 

and 


(II. 11 9a) 
(II. 119b) 


Here the dimensionless coordinates are described in a manner consistent with 
the definitions for (£, 7 ), £). 


One statement of caution should be made before leaving this section. 

That is, the expressions immediately above apply only to the relative positioning 
of P as measured from Q; and, then only for the local rotating (or horizon 
oriented) triad reference system. Normally, the premise would be that the 
position coordinates are measured from P to Q; hence, the two position angles 
are appropriately defined in Eqs. (11.116). 
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A Powered (r ) Intercept Problem. Earlier , in Reference [l], the problem of a 
two-particle intercept was introduced, solved and illustrated. There, both the 
impulsive action and a fixed magnitude thrust intercept case were documented. 
However, with the introduction of a second frame of reference (the inertially 
oriented one), in this present work, there is reason to go back and reexamine 
the intercept problem. The purpose in doing so is to resolve it for the case of 
a fixed thrust aligned with the inertial reference frame T s triad. 

For clarity and continuity in notation the previously solved thrusting- 
intercept result will be restated (or alluded to) here. The reason for doing so is 
to introduce the notational scheme employed now; but, to do so by means of a 
known resultant. Hopefully this will alleviate any difficulties which could arise 
in notational understanding. 

The meaning of intercept, as inferred here, was described in Article 
III. 6, Reference Cl] ; and, the powered intercept problem was solved in Appendix 
D (see Section D.6.2). Briefly, there one finds that Eq. (D.12b) represents 
(symbolically) the linear solution for a relative motion with thrusting (T). Intro- 
ducing the concept of intercept, then that solution was represented by Eq. (D.13); 
and, the dimensionless resultant, corresponding to this, was written out in Eq. 
(D. 15). 

Now, for the present case (an intercept to occur by an inertially directed 
dimensionless thrust, T ) the solution equations to be manipulated are noted as 
Eqs. (11.30), herein. Symbolically, the general solution can be noted as: 

X = A X +BH' + C T r . 
o o I 

(Here X is a representation for position; X f represents velocity; and A, B, C 
are the appropriate coefficient matrices. For reference purposes, the corres- 
ponding velocity expression, Eq. (II. 31), is expressed as: 

x’=A r X +B'X t +C’T . 
o o I 


(II. 120a) 


(II. 120b) 
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The primes, here, infer differentiation, as before. The "zero” subscript is 
used to signify an initial valued quantity). 


For a solution, the concept of intercept in applied (X = 0), and the value 
of t is determined accordingly. For consistency in notation, the intercept is 
assumed to occur in t* seconds (or in a transfer of 9 * radians); thus, this 
solution for is denoted as ■ T* . Consequently, the symbolic description for 
the intercept is written as: 

T* = - C _1 (A X +BX’), 

I ' o o 

wherein C 1 is the inverse matrix corresponding to C. An expanded repre- 
sentation for this resultant is written as: 


T* 


'4 


\v I 


-1 

0 

-1 

0 

T H 

= -(C A) 

*0 

-(C B) 

’’o 

^ * 

. Z - 




/o_ 


A more useful expression, for present purposes, is set down into the following 


form: 


\ 


r*A= A*x + b*x ! = 
loo 



wherein A is related to the determinant of C \ and A*, B* are the results 
of matrix multiplication noted in Eq. (II. 121b). 


The details of the required matrix operations (above) are not included 
here; however, the scalars for the two matrices (A*, B*), and the value of A 
is noted immediately below. That is: 


a il ” ^ ^ C ° S ^ ~ Sln + 0- “ cos <P) C6(l - cos 9) +23] - 12<p sin 9 , 
2 

a* = 3 [3(9 sin 9 + 2(9 cos 9 - sin 9 ) - (39 + 2 sin 9)(1 - cos 9 )] , 


( 11 . 121 a) 


(II. 121b) 


(II. 122a) 
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a 31 = °* 


a* 2 = ^ 003 V- sin <P)> 

a 22 = “ sin <P + 2(1 - cos <p). 


a* = 0 = a* = a* , . 
32 13 23 


a* - - Acos (p/(l - cos <p); 

oo 


b* =■“(!“ cos <P)(3 <p - 2 sin <p), 

r 3 2 

b* = 4[- sin <p ((P - sin <p) - (1 - cos <p) ] , 


b* = 0, 

31 

b* = 4(4 - cos <p )( 1 - cos <p) - ^ (7 sin <p - 3<P cos <p), 

JL ^ ^ 


b* = 2[~^- sin P “ (3<P + 2 sin <p)(l - cos <p)] , 


and 


b* = 0 = b* = b* , 
32 13 23 


b* 3 = - A/tan | ; 

A . 2 9co 2 

A - 9(<p - sin ip) + (1 - cos <p) + 10(1 - cos <p) , 

Ci 


(II. 122b) 


Equations (11.122) describe the magnitude of T, thrust n needed to produce 
an intercept from an initial state (x^, x^). As before, the Tf time M required for 
the operation is presumed known (a priori). Once the resultant r* is acquired 
a ’time" history of the displacement (X) and the velocity (X 1 ) is obtained from 
Eqs. (11.120)*. In Eqs. (11.120) the proper value for T is (obviously) that 
solved for above ; namely, T * . 

*These results, being analytical, infer a linear solution to this problem. 
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This intercept solution describes the state of motion as it would be 
"seen" in the local rotating frame of reference, even though the "thrust" (Tp 
is aligned with an inertially oriented triad. There is a corresponding, or com- 
panion, solution for this situation — one which would refer the state to an inertial 
frame of reference, per se. In order to acquire this resultant one could proceed 
in the same manner as is indicated by the generalization alluded to in Eqs. (II. 33b), 
and the procedure used to obtain Eq. (11.36). These particular operations are not 
to be carried out here; the interested reader can obtain the desired results, readily, 
but at the expense of some needed algebraic gyrations. 

A General, Point-to-point, Transfer. The general point-to-point transfer in- 
volves the application of an impulse, and/or a specified thrusting action, which 
will drive an orbiting particle from some initial state (X^, Xp* to some final 
state (X^, xp. Normally, the final state is described in terms of either position 
(xp or velocity (xp, as a priori information, but not both. This is the conse- 
quence of preselecting the ff time" (or transfer angle) during which the operation is 
to occur. 

In effect the solution (s) here follow, identically, the pattern set down in 
the foregoing section. For this reason the procedures, here, will be abbreviated; 
however, the cognizant results will be set down in sufficient detail to evaluate 
solutions, if desired. Following the established pattern the present solutions 
evolve, generally, from the state equations written in rotating frame coordinates; 
nevertheless the companion solution (in inertial relative coordinates) can be ac- 
quired as above. There is an exception to this procedure. This is for the case 
of an impulsive transfer which is referred to the inertially oriented frame of 
reference. 

(a) The Impulsive Transfer. The impulsive point-to-point transfer 
is developed from the results in Eqs. (II. 30) and (II. 31). For instance, denoting 
the final state as (X^ } xp, then from these solutions, one can write: 

♦The quantities (>L, x!) are employed as general relative motion state vectors. 
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(II. 123a) 


and 


x_ = A x + B x* , 
f o o 


xl = A'X + B’X r . 
I o o 


Here it must be assumed that A, B, A f , B 1 are known functions of the "time 
required” to complete the intended maneuver. 

In order to ascertain the impulse (X T ) needed to reach x„ from (X , X 1 ) 

o f K o o 

one should solve Eq. (II. 123a). The form of that solution is: 


X' = B _1 X -(B _1 A) X , 
o f o 

wherein B ^ is the inverse of matrix B (see Eq. (II.123a)). Remembering 
that the matrices are known functions of the transfer time (and/or angle), then 
the resultant is acquired as follows: 


Since X. -X(£, 77, £).; X r = X T (£ * , 77 '*, £ r ), the solution is ex- 
1 ' ' 1 o o o o o 


pressed as: 


X 




V 

K 

•ri 

„ .Q l; 

rHj<J 

11 

*f 

X M 

A L 13 j 

% 

K 

w m 


c 

. f . 


c 

0 

• 


for (i, j) = (1, 2, 3). 

Herein the various scalar components of the matrices are: 

4 cos <P -2(1 - cos <p) 0 

2(1 - cos <p) sin<P 0 

0 0 -p- 

sin (p 


b.. = 
ij 


X 


and 


*Subscripts (~) Q and (~)^ denote initial and final values, respectively. 


(II. 123b) 


(II. 124) 


(II. 125)* 


(II. 126a) 
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, wherein 


(11.126b) 


’ll 

a 

12 

a 

13 

21 

a 

a_ 

22 

23 

'31 

a 32 

a 

33 


a il = ^ “ 4 cos + SC 1 - “ cos <P)[(39 -4 cos <P) - 4(<p - sin <p)~\ , . 

a = 2 [3V 3 sin <P - 7(1 - cos V 5 )] , a - 0, 

21 Id 

a 2i = 2 (! - c °s V5), a 22 =-sin<P, a 32 = °» 

a „ = = 0, a “ - A/tan <p , and 

13 23 33 

gcp 2 

A s 4 [sin <P (cos <p - ^ ) + (1 ~ cos P ) ] . 

The consequence here is that the general point-to-point transfer is now 
known; and that the sought for impulse is described by Eq. (11.125). 

Should the desired solution have been one to determine an impulse (x^) 
to reach a final velocity (x^); then, that the result would have appeared as: 

X» = (B f f 1 X' - [(B^'^A 1 ] X . 
o f o 

(Recall that primes in this expression denote differentiation with respect to the 
transfer angle). 

The results above account, in general, for an impulsive transfer from- 
an initial state to a prescribed final state. However, the specific case set down 
here was for a problem referred to the local rotating frame of reference. 

The corresponding results, cast in inertially described relative motion 
coordinates, could be acquired by means of the relations noted in Eqs. (II. 32a); 
or, the more direct approach could be taken using the solution expressions, 

Eqs. (11.39), (11.40). Making use of these latter formulae, and applying to them 
the concept noted in Eq. (11.124) — but introducting the inertially described 
coordinates — then: 


(II.126c) 


(11.127) 
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with K. = X (H, H, Z).; X T =x r (S', H', Z f ) 
I 1 o o o o o 


the solution would be expressed as: 


r~t ■ 

~o 

H' 

0 

= 1 M 

■ ■ 

*f 

H f 

*1 w 

* « 

~0 

H 

0 

Z ! 

0 

* p* 


1 

N 

1 


0 

N 

i 


, for (j,k) =(1, 2, 3). 


Herein the two coefficient matrices are: 


M s 


(2 sin9+ sin 2 9 
-3(p cos <p ) 


[2 cos 9 + (cos 29 -3) 0 

+39 s in 9 ] 


2 cos 9 - “ (3 + cos 2 (p) 2 sin 9 - — sin 29 


and 


w ■ 


2 
0 

39cos9“4 sin9 
3 [9sin9-2(l - cos <p)] 
0 


0 

0 A/sin 9 

3 [9sin9~2(l -cos<p)] 0 

-sin 9 0 

0 -A /tan 9 


with the quantity A being: . 

A = 8(1 - cos 9) - 39 sin 9 . 


This last result describes the same transfer situation as does Eq. 
(11.125); the exception is that here the coordinates are inertially oriented — the 
previous case was for motions referred to the rotating frame of reference. 

These two solutions describe (in a linearized sense) the general impul- 
sive transfer from one prescribed relative position to another. It must be 
recognized that the primary piece of information acquired here is the required 
initial impulse ; the final impulse, and a time history for the state variables - 


(II. 128) 


(II. 129a) 


(II. 129b) 


(11.129 c) 
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during these transfer maneuvers - must be obtained from the general expressions 
(see Eqs. (11.30), (11.31), (11.40)). 

The next problem area to be examined is also concerned with a point-to- 
point transfer, but there the primary consideration will be a determination of 
the thrust required to accomplish the maneuver. In this next sub -sect ion the 
results will be reported only for cases where the coordinates are described in a 
local rotating frame. To determine the transfer in relative inertial coordinates 
the reader should apply the non-familiar transformations for coordinates and 
speeds. 


(b) The Thrusting Transfer. The outline which follows will show re- 
sults obtained for a point-to-point transfer due to thrusting actions. Here, as 
mentioned earlier, the two cases of interest will consider the thrust aligned with 
one or the other of the reference triads. However, both solutions will be for the 
state variables referred to the local rotating (horizon oriented) triad only. 

Since the basic expressions to be manipulated stem from the formal 
solution results, these will be symbolically represented, at first; the specific 
equations will be noted when final statement formats are set down. 

In this regard a general form of the state equations, referred to either 
frame of reference, can be written as: 


x = A x + B x 1 + CT 
o o 

and 

X» = A r X + B »X» + C r T , 
o o ’ 


where the X. are state vectors, T is the dimensionless specific ”thrust ,T , 
primes denote derivative forms and A, B, C, etc. are coefficient matrices. 


Quite simply, the solution, that of defining a thrust to accomplish the 
maneuver, can be represented by: 

r=C“ 1 X.-C" 1 [A X +B x 1 ] . 

I O o 


(II. 130a) 
(II. 130b) 


(11.131) 


180 



This result suggests that (now) a more general solution type is acquired since 
the complete beginning state (X q , xy is present. The fact that only the terminal 
position (xp is preset is a necessary consequence of the general solution. These 
being fixed thrusts (T), then the terminal velocity (x|) will be acquired from a 
time -history description of the motion state, per se. 

To make Eq. (11.131) operationally useful it is necessary, next, to set 
down the results in an appropriate algebraic/matrix format. This is carried 
out below. 

(1) The solution for T - T (T h , T , T*). On the assumption that r 

b *i & 

has the designation indicated, then following the solution described in Eq. (11.131), 
it can be shown that the various coefficient matrices are obtained as indicated 
below: 


First, the solution is cast into the format; 


» - 


r 

m m 


r 



• 

r « 










1 

r i 


r i 



r -j 



T v 

A ' 

L C ij J 


+ a 
L I j - 


V 

o 

+ L b ij] 

% 


k. 


» 

£ 

m ^ 


c 

o 


A 

* 


for (i, j) = (1, 2, 3). 
Here, 


C 11 

- 4(1 - cos 9 ) ■ 

3(p 

2 

'» c 12 = “ 2 ^-sin9), c 13 = ° } 

C 21 

= 2 ( 9 - sin 9 ), 

°22 = 

: (l-cos 9 ), c 23 = 0 , 

C 31 

= °> C 32 = 0 * 

C 33 = 

A / (1 - cos 9); 

a il 

= 4 [ 6 (psin 9 - 7(1 - cos 9 )] , a = 2 ( 9 ~sin 9 ), a = 0 , 

1^ 13 

a 21 

- -2 (<p~ sin 9 ). 

a 

22 

= - (1 - cos 9 ), a 23 = 0 , 


(11.132) 


(II. 133a) 
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(II. 133b) 


a 31 =0, a 32 =0, 


Next, 


3<p 

b = -£- sin(p-4<P (1 - cos <P), 
11 2 


and 


b =2 [2(1 -cos<p) -<psin<p] , 

6 X 


b =0, 

31 * 

b =2 [7<Psinp- 8(1 ~ cos<P)] - (1+cos <P) t 

XCt 2 j 


b 22 = “^( 1 -c°s(p) > 


b =0, 

32 * 

b_ =b oo = 0, 
31 32 ’ 


b -- A sin<P/(l -cos <P). 

OO 


For the above, the quantity A is defined by: 


» _ 5 (p 3 

A = 8(1 - cos(0-<psin<p) + (1 + — cos <p). 

£ o 


(2) The solution for T = r (T^, T , T y. In this solution example the 
motion traces and the transfer arise from the application of a non-dimensional, 
specific, fixed thrusting action; but one where the components are aligned with 
the inertially oriented triad of reference. 

Once more the solution format is that described in Eq. (11.131); however, 
the coefficient matrices are not identical to those set down above, for obvious 
reasons. 

For identification, let the present solution be given in the form: 


(11.133c) 


(II. 133d) 
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(11-134) 



for (j, k) =(t, 2, 3). 

The various scalar components of the coefficient matrices are: 

3 

c = S^sin^-Sfl -cos^), c =— (9cos<p- sin<p), c =0, 

IX X^ & Xo 

°2 1 =3 C 2sin ‘ ifl '^( 1 +oos ^)]> c 22 = 2 sln<lfl ' 2 ( :L ■° os < a )> C 23 =0 ’ 

C 31 =0, C 32 = 0 ’ C 33 =A/(1 ‘° OS<P) - 


Next, 


11 


21 


a 


31 


12 


22 


32 


13 


2 

= 9P(^cos<p-sinp) + 23(l -cos<p) +6(1 -cosp) -1293sin<p, 

2 

= 3 t2 (<P cos <p- sin <p) + 3<p sin 90- (1 - cos 93) (3^ + 2 sin^)], 

= 0 , 

3 

= — (sin 90-93 cos 93), 

3<p . 

= 2(1 - cos <p) - — - sin 90 » 

= 0 , 

= a =0, and a = - A cos 90/(1 - cos 93). 

OO 


Lastly, 

b = (1 - c °s 9 s )(2 s in (p - 3<p ), 

b -4 [“ sin<P(9- sin <p) -(1 -cos<p) 2 ], 
21 2 


(II. 135a) 


(II. 135b) 
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(II.135c) 


3<P 

b =4(1 -cos<p) (4 - cos(P) + - (3(0 cos<P -7 sin (O), 

J. ^ " 

2 

^22 = _ 2(1 -cos<p)(2 sin<p + 3<p), 


b = 0, 
32 


b 13 =b 23 =0> ^ b 33 = - AAan f • 


In the above expressions, 

2 9 2 

A = 9(<p - sin<p) + — — (1 -cos (p) + 10(1 -cos <p) . (II.135d) 

Equations (11.132) and (11.134) describe a "thrust" vector needed to 
accomodate some preselected point-to-point transfer -- one which is to occur 
in a known time (or, over a given transfer angle). To describe the T time -history" 
of such a maneuver, it would be necessary to insert the T vectors into the 
proper state equations (set down) and to evaluate those accordingly. Once again, it j 
must be remembered that these are linear resultants and, correspondingly, they 
have a limited applicability for prediction purposes. 

Also, as mentioned earlier, these results are applicable for local horizon 
state coordinate descriptions only; if the reader desires to have a similar result 
in the inertially defined relative motion coordinates, it will be necessary to trans- 
form these results over to that frame of reference. 

In the next sub-section an operational type example will be described. 

The purpose there will be to illustrate how the information developed here can 
be put to use — a mathematical solution is not set down. 

Example: A Stationkeeping Mode. The idea of "stationkeeping" as it is inferred 
here, has to do with maintaining some position* which is in close proximity to 
(say) an orbiting vehicle. Previously, in Reference [l], ideas were mentioned 
which fit this concept. There the elimination of secular influences allowed for 

*position, here, suggests a region rather than a point . 
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adjacent relative orbits, which fit this idea. The "messenger” transfer, as 
depicted on Fig. III. 7 (Ref. [l] ), could also be classed as a candidate operation - 
even though that maneuver was time restricted in its stationkeeping mode. 

To be more realistic and to present a rather general scheme of operation, 
the maneuvering concept set forth in Reference [ 22 ] will be utilized here. This 
will allow the reader to make ample use of the motion traces, described for a 
general relative motion, and to ascertain how a stationkeeping mode could be 
established. 

For present purposes the reader is referred to (say) Fig, III. 3 in Ref. 

Cl]; and, more specifically, to trace (4) on page 73, and to the hodograph 
(trace (4) on page 81) of this report. These two sketches are reproduced below; 
they will be utilized in describing an in-plane stationkeeping mode - one which 
will "hold” a relative motion in the near vicinity of an orbiting spacecraft. 

On Fig. II. 14a, the origin (P) is assumed to be the locus for the orbit- 
ing spacecraft while Q is an instantaneous position for the relatively moving 
orbiter. Fig. II. 14b is the (relative) hodograph corresponding to the trace on the 

In keeping with the geometry 
shown here, the stationkeeping mode 
will consist of the trace curve (a, b, 
c); the idea is to repeat this trace 
geometry, in a continuous fashion, 
thereby retaining a closed local 
"orbit” for Q, near to P. Obviously, 
if some corrective action is not in- 
troduced then Q will continue to 
fall further and further behind P, 
during its free motion. 


first sketch. 



Fig. II. 14a. Sketch of a general, 
in-plane relative motion trace. P 
is the reference (spacecraft) 
particle: Q is the local ’’orbiter”. 
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Focusing attention, now, on 
Fig. II. 14b, one notes that the hodo- 
graph (for the trace (a, b, c)) is also 
noted by a, b, c (here). That is, at 
position "a", the hodograph locus shows 
particle Q to have a relative velocity 
composed of , -r£) speeds. At 
position T V r the hodograph suggests 
speed components (0, 77 ^); while at 
position n c” the hodograph is similar 
to that at H a ,f except that now the 
speeds are (+4^ , - 77 ^). 

To retrace the position plot 
(a, b, c) it is only necessary to re- 
establish the hodograph locus at ,r a M . 
This can be simply done by adding an impulse at "c M (of size -24 *) — then the 

£L 

initial hodograph is reformed, and the trace (a, b, c) is followed again. 

Consequently, each time particle Q reaches position H c M , the impulse 
is applied again and path (a, b, c) is restored. This obviously establishes (in 
a hueristic sense) the stationkeeping mode for this displacement configuration. 

Some thought on this operation will convince the reader that a like stationkeeping 
path (for Q, relative to P) could be determined for the entire (in-plane) space 
adjacent to P. Of course, in a more restricted sense a stationkeeping mode 
could be developed, adjacent to P, which would eliminate the repeative applica- 
tion of impulses. Such a scheme would be reminiscent of the relative motions with 
secular terms removed. This stationkeeping operation would be repetitive and 
without impulse input, but it would be cyclic in the spacecraft (P) period. The 
operation described above does not utilize the entire hodograph, hence that opera- 
tion is repetitive at less than orbit period. 



Fig. II. 14b. Hodograph for the 
relative motion above. Note that 
the velocity for "Q" is shown, and 
that the speed at a, b, c are des- 
cribed. Also, recognize that 
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Needless to say, there are a variety of operational schemes which can 
be developed from these ideas — the inventive reader could easily envision a 
number of such. 

Here, a word description has been used to develop a stationkeeping mode. 
The mathematics of such a scheme would readily follow from study and imple- 
mentation of the accompanying sketches, coupled with several of the formulae 
provided herein. It is hoped that by now the reader has acquired an appreciation 
and use for the relative motion traces and hodographs. 

Example: A Minimum Time, Relative Transfer . The problem to be solved, next, 
has to do with a minimum tune transfer. The transfer is to begin at an orbiting 
spacecraft (P), and proceed to a lower altitude orbit; but, to a position which is 
close by to P (close enough so that the linearized relative motion equations will 
be applicable). 

Such an operation, as that described here, could depict the ’'rescue" of 
a stranded astronaut from a position below the parent vehicle (see Fig. II 0 15). 

If time would be of the essence, here, then the information gleaned (below) 
could be of great value. 

For mathematical manipu- 
lations this problem will be cast 
into the format set down in Eqs. 
(11.30). That is, all state quantities 
are non -dimens ionalized, and the 
minimum time requirement becomes 
a minimum transfer (A <p) require- 
ment. 

In the parlance of present 
nomenclature, the following quantities 
are employed: 

t- ^ s -i ; T? = C = o ; 

p 





Q 


Fig. II. 15. Sketch of the minimum 
time transfer operation. P is the 
(circular) orb iter; Q is the final 
position to be achieved, Av is the 
impulse, and fi its direction. 
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(11.136) 


W C o = 0; 

£ 1 - sin/? = v sin/3, 

° V P 

17* " - V cos /3; £ f =0. 
o o 

From Eqs. (11.30) this relative motion problem will be expressed (in 
general) as: 

XflP)=A X +BX' +ct ; 
v o o 

and, in particular, for the in -plane format, write: 


» * 
-t 


sinp 

2(1 - cos<p 


m m 

-v sin fi 

0 

■ J 


_ -2(1 -cos<p) 

4 sin(P-3<p 

. . m 


-V cos /3 


From this expression it is found that a general solution statement is: 


• V A 

— - 2 cos p= sin<psin /3 - 2 cos <pcos . 

Examining this expression for an extremal (presuming that~iv and v 
are not variable g) then the condition for an extremal is found to be: 

tan 2 = 2 Ctn i3. 


describing a requirement between the "aiming" angle and the transfer angle. 


When this requirement is incorporated into Eq. (11.138), it is found 
(after some manipulation) that the minimal maneuver must satisfy the con- 
straint condition: 



4 cos /? 


V 2 

4-3 cos /? 


♦No T-terms are apparent since " thrusting" is not assumed. 


(II. 137a) 


(II. 137b)* 


(11.138) 


(11.139) 


(11.140 
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(Note that, according to this result, - ^ 4, in order to satisfy the minimal 
condition). 

With this requirement, above, one sees that; 

(1) If "ft and v are known, a priori, then Eq. (11.140) shows a con- 
straint on for the minimum transfer (time). 

(2) If 'ft and /3 are known, a priori, the equation defines a constraint 
on V for the minimum transfer (time). 

This example concludes the work in this section. Here, the attempt 
was to show how practical results could be extracted from the general solutions ; 
and, to illustrate, by example, a few operational situations which could be 
simulated and studied as by-products of the more formal developments. With 
this insight into relative motion, as might be gained herein, it is obvious that 
there will be many opportunities and situations which can be examined and 
studied each making ample use of the solutions which have been found and set 
down in this report. 

In the following section a few remarks, as concluding statements, will 
be offered regarding this investigation. 



XII. CONCLUDING REMARKS 


The work which has been reported in this document represents the 
results gleaned from an analytical study on a class of relative motion problems. 

The general case here considers two particles orbiting about a central attract- 
ing mass. There is no mutual attraction between particles assumed, and, one 
particle is designed to travel on a circular orbit; this serves as a reference (or 
datum) basis for the system. 

The major portion of effort in this part of the investigation was directed 
to obtaining general solutions for the relative motions. These evolved as a con- 
sequence of an initial state (r , V ) and in response to a system of prescribed 

oo 

forces. Of course the results must be considered as reliable only in the vicinity 
of the reference particle; this is a condition which is imposed in addition to the 
constraint of a simple central attracting center. Nevertheless, the equations 
obtained here are invaluable to understanding the types of motions (variations in 
state, relatively speaking) which do occur. And, in view of the fact that these 
results describe variations in the two primary frames of reference, it is apparent 
that what has been determined here represents a most complete and comprehensive 
data base. 

Aside from acquiring the equations, they have been put to use in describ- 
ing traces of the motion on the principal coordinate planes. This graphical repre- 
sentation adds a new dimension to the results, providing the reader with a pictorial 
as well as analytic description of the influence played by initial values and by the 
applied disturbing forces. It is the investigator's belief that the data summaries, 
which are made up of these sets of information, will materially aid in the study 
and understanding of relative motions. 

Even though the results which have been acquired here are complete, as 
far as they go, they have considered only one particular class of disturbing 
forces. Consequently, it is not unreasonable to expect that the tasks commenced 
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here should he continued; at least to the point of examining other input systems 
to ascertain their influence of the motion. Hopefully this latter task will be 
undertaken and completed. 

In addition, the work described here, and that which is expected to come 
later can be applied to real world problems — or simulations — so that many 
more realistic situations may be examined. It would be a mistake to presume 
that the equations , developed from this investigation, have 'been used"; and, 
that they can be discarded now. Instead, they should be put to use for explicit 
and for implicit data acquisition. Here, in the body of this report, and in the 
references cited, an interested reader can find other suggestions on how the in- 
formation can be utilized. 

Once more the reader is reminded that this investigation has given equal 
emphasis to the relative motion velocity, as compared to the displacements. 

This was done because it is felt that the full state of a motion is essential to a 
complete understanding of these problems and for their manipulation. Generally, 
orbital motions cannot be altered without a change in the total state ; and many 
"changes in state" are due to either a direct or implied change in the velocity. 

Finally, it is remarked that the work described here is complete insofar 
as it goes. Obviously there are other areas and problems which should be 
examined; it is hoped that some of these, at least, can and will be carried for- 
ward to successful conclusions. 
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APPENDIX A 


Definition of a Rotation, Consider an initially oriented Cartesian frame of 
o 

reference, (~) , and a second reference frame, (~ )', which is located by an 
angle 0^ with respect to (~)°. Let this angular displacement be described as 
a T, right-hand rule" rotation; and, let the rotation take place about the z-axis. 
As a consequence of these conditions only the (x, y) plane will undergo a re- 


orientation during the displacement. 



x' 


For convenience, in the description 
and representation of this situation, assume 
that these Cartesian frames are assigned a 
set of orthogonal unit vectors each. To 
identify them, a superscript is added to 
the rotated triad. Therefore, after a dis- 
placement, of (say) 0^, the unit vector 

e , is related to the inertial (e ) vectors 
x 1 i 

by: 


Fig. A.l. Sketch of a displacement. ' 


V " (6 X ' e x' )e X + (e Y ‘ V e Y + (e Z ‘ V 6 Z 


- cos 0 e + sin 0 e + 0 e , 

IX 1 Y Z 

In a corresponding manner the other unit vectors may be described by: 


y 


sin 0 e x + cos 0^ e y + Oe , 


and 


e = e. 


Treated as a matrix operation the expressions above can be written as 


(A.l) 


(A. 2) 
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1 

i CD 

1 


cos 0^ 

sin 6^ 

0 

■ 


~ e x 

e > 

y 

= 

-sin & 

cos 0 
1 

0 




e t 

L z 1 J 


0 

0 

1 

m 


. e z . 


e. t = T(0 + ) (e.) 0 , (i = x, y, z) 


(A. 3a) 


(A. 3b) 


The notation in Eq. (A. 3b) implies that a +-rotation is needed to locate 
the ('v) , -triad from the (^)°-triad. In order to reverse this situation - that is, to 
provide for a reverse rotation - or, to locate the (~)°-triad relative to the (~) f - 
triad - one could write (symbolically), 


(e.)° = T(0 + ) 1 e* f , (i = x, y, z). 

Here T(0 f )^ 1 is the inverse matrix (operator) obtained from T^)! (Note: 


since 

T(0 + ) is an 

orthonormal matrix, then its inverse 

is identical to its trans- 

pose. 

Consequently, 







- 

cos 0 

-sin 0 

0 





1 

1 




in 


sin 0^ 

cos 0 
1 

0 

s T(6~); 



_ 

0 

0 

1 



(A. 4a) 


(A .4b) 


T 

that is, the transpose, T(0 + ) , is identical to the "transform matrix" developed 
for an angular displacement, 0 = - 0^ ). 

There are some interesting and useful relations which can be developed 
for these transforms. Several of these are shown below: 


(a). Demonstrate that T(0”) is indeed the inverse of T(0 + ). 
Proof: If the statement is true then; T(0“) * T(0 + ) = I ; 

O 

that is, 


196 



(A. 5)* 



c 0 

-s0 

0 


’ c e 

s 0 

0 


- 

1 

0 

0 

I = 

60 

c0 

0 


-s6 

c0 

0 

= 

0 

1 

0 


0 

0 

1 


0 

0 

1_ 


_ 0 

0 

0 

-J 


wherein (~) infers the full (3 by 3) unit matrix (Q.E.D.). 
o 

(b). Evaluate the time derivative of the transform operator, T(0"'). 
In this regard: 



C0 

-s 0 

0 


"-S0 

-c 0 

I 

0 

T< 9 "> s ft 

s0 

c0 

0 

• 

= 0 

c0 

-s0 

0 


0 

0 

1 


0 

0 

0 

J 


Also, note that: 


C0 

s 0 

0 


-s 6 

c0 

0 

-s 0 

c 0 

0 

= 0 

-c 0 

-s0 

0 

0 

0 

1_ 


0 

M 

0 

0 _ 


T < 9 > a ft 


Now, having the derivative (s), define a matrix, B, such that T(0”) s BT(9~) 0 
Using the as sinned relationship it is apparent that. 


0BT(0~) T (0"f 1 s B0 = T(0")T(0 - )" 1 , 

^ -i 

at least symbol ically= However, T(0“) = T(0 + ), therefore; 



“-S0 

-C0 

0 


c0 

s 0 

0 

0B = T (0 r )T(0 + )=0 

+c0 

-s0 

0 


i 

01 

cr> 

c0 

0 


0 

0 

0 


0 

0 

1 


which leads directly to, 

*For these operations a shorthand notation is adopted for the circular functions; 
here, s0 = sin 0 and c0 = cos 0. 


(A. 6a) 


(A. 6b) 
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B = 


0 

-1 

0 " 



1 

0 

0 

= 





0 

0 

0 


0 


For reference purposes, it should be remembered that 


T(0") = 0BT(0“). 


For the converse operation define T(0 + ) in a similar manner, and re- 
late this operator to the transform T{0 + ): that is, let 


T(0 + ) s 0 B + T(0 + ); 


then 


or, 


therefore, 



~-s0 

c0 

0 


c 0 -s 0 

0 " 

T(e + )T(e>e 

-c0 

-s0 

0 


s0 c0 

0 


0 

0 

0 


0 

0 

1 











0 

1 

<f 



sb + = e 

X 

-1 

0 

0 

= 6(-B); 




0 

0 

0 




T (0 

h ) = - 

0BT(0 + ) 






(c). Develop some properties for B: 

First note that the product of B with itself is : 


*In view of this definition one can recognize that 0B is a specialization of the 
matrix operator used to replace the vector operator "co x ”, where <o is an 
instantaneous 1 dotation vector". For the present case 0B represents an 
CO - co(0); thus the two elements in this skew- symme trie operator (B). 


(A. 6c)* 


(A. 6d) 


(A. 6e) 
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where I is the two -element idem matrix! 

A 

Making use of the definition for B (Eq. (A, 6c)), then Eqs. (A. 3b) and 
(A, 4b) may be rewritten as: 


T(0 ) = I cos 0 + B sin 0 + (I - I ) a (TJ , 

A o A 


W T(0 rf ) = l 2 oos e-Bsme+(I 3 -I 2 )s(T + ). 

Next, evaluate the operational quantity, T(0~) * B *. T(0 + ): 

<T_)B(T + ) = (cos 0 - 1) + B sin 0+ I 1* B * |~I (cos 0 - 1) - B sin 0+ I ] 

= |I (cos 0 - 1) + B sin 0 + I 1 * f~B (cos 0 - 1) + I sin 0 + B 1 

u Z -J L_ A A -J 

= (cos 0 - 1) + B sin 0 + 1^ J * JjB^ (cos 0' - 1 + 1) + I sin 0 J 


B (cos 0 - cos 0) - I sin 0 cos 0 + B cos 0 + I (cos 0 sin 0 

A A A 


sin 0) + B sin 0 + I sin 0 = B ; 

A A A 


(A, 7a) 

(A. 7b) 
(A. 7c) 
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therefore, 


(TJB 2 (T + )= b 2 . 

Note that past and pre-multiplication of this result leads to the following: 
That is, since (T + ) - (T_) \ and vice versa; 


or 



(T+) [<T_) B 2 (T+)] (T_) = (T + ) B 2 (T_) 
B 2 = (T+) B 2 (TJ. 


(d). The composition for two subsequent rotations , about a common 
axis can be expressed by 

T(6 1 + ) * T(0 2 + ) = T(9+ + e+). 


This is an a priori statement. The proof of the statement follows directly from 
Eq. (A, 7c) wherein; 

T(6+) * T(0p =[l 2 (cos 6 i - 1) - B sin 0 X + I g ]* [l g (cos 63-1) 

- B sin 6 + I 1 ; 

2 3 J 

which, after manipulation, is noted to be, 

= I 2 [cos (0 1 + 0 2 ) - 1] - B 2 sin (0 1 + 0 2 ) + I 3 ; 
or. 


T(0 + ) * T(9p = T(0 + + 0+) 

(Q « E „ D » ) 

(Hueristically, it is apparent that this technique may be used for a sequence of 
displacements: Y 0.). 


(A. 7d) 


(A, 7d2) 


(A, 7e) 
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matrix. 


(e). An application of the above result , applied to the T(8 + ) transform 


According to Eq. (A, 3a) and/or (A. 7c), T(0 + ) when 0-0, gives 

T(0 + ) £ = T(0“)J= I , Correspondingly, for |0 | = 77/2 and 77, respectively; 

T (77/2*) = I - I - B, while T(7r/2") = I ~ I + B; 


3 2 


3 2 


and 

T(77 + ) - I 3 - 2r 2 and T(iT) * - 2I g . 

Now, for instance, T(77 + ) = T(77/2 + ) * T(77/2 + ) - (I - I - B) * (I - I - B) = 

O Li O & 



Correspondingly, the reversed rotations would lead to a like resultant. 


The Consequence of a Rotation. The idea to be illustrated here is that the re- 
lation of a position -vector, projected back onto a reference set of axes, can be 
defined by ineans of the transformation operators, T(0 ± ). For the example 
(here) let the (~)°-triad be an inertial triad (of reference) as before. Also 
suppose that the (~) ! -triad is a rotating frame of reference; one following a 
particle "P". (A graphical description of this situation is been below). 



Fig. A. 2. 


The condition which will be illustrated 
considers a simple rotation of the position 
vector R p ; a displacement through the 
angle, 0. To be described is a representa- 
tion of this displacement in the (~)°-frame 
of reference. 

Since the vector R p (t) can be de- 
fined by 
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V> 5R P e ' x - 

in the moving (or displaced) frame, then its projection onto the (~)°-triad is 
(see Eq. (A.l)) 

H p (t) = R p (cos 0e x + sin 0e y ). 

This is equivalent to transforming from the (~) '-triad back to the (~)°-triad; 
or 

(Rp) 0 =T(e-)(R p )>. (A. 8) 

This is a description of how the transform operators may be employed 
to relate a rotated vector back to an inertial frame of reference. 

Obviously, the operator T(0 + ) should be employed whenever it is desired 
to transform a vector from the inertial -frame to the displaced-frame of reference. 

That is, 

(Pp) 1 s T(0 + )(Rp)°, (A. 9) 
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APPENDIX B 


Some Special Matrices . In this appendix some special matrices, and some 
elementary operations with them, are described. Specifically these matrices 
are introduced here in order to maintain continuity of the main text. 

Ordered Unit Matrices. The general nth order unit matrix is defined as an 
(n x n) unit diagonal matrix. That is 

I 

n 



(B.l) 


In this regard, 3-dimensional problems should lead to a unit matrix, I , defined 

O 

as: 




(B. 2) 


Now, special variations of I lead to the matrices I and I , which are 

3 £ 1 


defined as : 


*2 = 


’ 1 

0 

0 


’ 1 

0 

0 ~ 

0 

1 

0 

, and I = 

0 

0 

0 




1 

- 



0 . 

0 

0 


0 

- 

0 

0 


, respectively. (B.3) 


Next, for convenience, three single unit matrices are defined. These are: 
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" 1 

0 

0 



*0 

0 

0 

J s 

0 

0 

0 

= L; 

J = 

0 

1 

0 

1 

0 

0 

0 

1 

2 

0 

0 

0 


and 


J 3 H 


- VV' 


(B.4) 


Operations with the Matrix, B . Making use of the special matrix, B , {or B), 
from Appendix A, the following operations are performed and identified. Since 


0 

“1 

0 



0 

1 

0 

0 

. 






0 

0 

0 

mt 


0 

. 

0 


then it can be shown that 



‘o 

0 

0 “ 


*0 

-1 

o' 

BJ r 

1 

0 

0 

-J B , and J B= 
2 1 

0 

0 

0 


0 

L 

0 

0 



0 



As a consequence, it is recognized that 

. BJ 1 =BI 1 =(I 2-V B = J 2 B 

and 

J 1 B = I 1 B = B(I 2 -I 1 )=BJ 2 

then, with 


{B. 5)* 


*This result is obvious since B = BI = B (J + J ) = B J 

2 2 v 1 2 7 1 


+ BJ and = BJ +J B, etc. 
2 2 1 1 
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(B.6) 


B 2 =BJ 1 +J 1 B[=J 2 B 


+ BJ 2 ] 5 BJ 1 + BJ 2 =B( Ji+ J 2 )-BI 2 , 


Making an observation of the matrices described in Eq. (B.4), it is apparent 
that 


Vi = J i =J i- 


J 2 J 2 = J 2> J 3 J 3 = J 3- 


while the mixed product of these special matrices vanish. 

Operations with a Combination of the Special Matrices. During the formulation 
and solution to a three dimensional problem, one is likely to have need for the 
inverse of a matrix formed by the special matrices defined above. To illustrate 
the nature of such, suppose there is a matrix combination like 

A H aj + bj + cj + eB 
12 3 


wherein (a, b, c, e) are coefficients describing the problem. Now, from the 
definition of these special matrices, it can be shown that 



by the addition of the matrices, as implied. In forming the inverse one needs 
the determinant of A (Det A = c(ab + e^)), and the adjoint transpose of A (=A^) — 
which can be shown to be: 



be 

ec 

0 


T 

A* = 

-ec 

ac 

0 



0 

0 

(ab+e^) 



(B.7) 


(B. 8a) 


(B. 8b) 


(B. 8c) 
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Now, since A * A^/Det A, then it is immediately apparent that 

be ec 0 


T 


e(ab+e^) c(ab+e ) 


-ec ac 0 


0 0 (ab+e ) 


which can be written as: 


A" 1 - 


bJ l + ^2 ' eB 2 + ^3 

2 c 

(ab + e ) 


In the above operations 


1 

-1 

0 * 


' B 2 

0 ~ 


’l 

0 

0 * 

1 

1 

0 

== 



+ 

0 

1 

0 

0 

0 

1 


0 

1 


0 

0 

0 


B + I, 


(B.8d) 


(B. 8e) 
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